式，此时可将裂缝作为 源汇项进行处理；
[0044]
(3)
[0045] 其中，Qi为窜流量；ΤΙ;为传导率；λ ;为窜流系统；A为裂缝面积；k为渗透率；<d> 为裂缝到网格的平均距离；Φ为势函数的大小。流程进入到步骤105。
[0046] 在步骤105,重点分析裂缝在基岩中的三类不同嵌入方式下数值计算格式，包括完 全穿透型、部分穿透型和交叉裂缝型，如图3所示；
[0047] 完全穿透时
[0048] 部分穿透E
[0049] 交叉裂缝型knne为裂缝渗透率；cT为两条裂缝间距离；Anne为裂缝开度乘以裂缝 长度；V为网格大小；x n为裂缝到网格边界的距离；kfl，kf2为裂缝渗透率；ωη，为裂缝 开度；L int为网格内裂缝长度；dfl，df2S裂缝中心到网格边界的距离。流程进入到步骤106。
[0050] 在步骤106,基于控制体积有限差分法建立嵌入式离散裂缝网络系统的数值计算 格式；
[0051]
⑷
[0052] 其中，
t裂缝节点处的压力值；ξ为沿 裂缝方向的局部坐标系；Τ为传导率；Vfi为裂缝孔隙体积；qnifi，qffi为基岩与裂缝、裂缝与裂缝 之间的窜流量。
[0053]
[0054] 流程进入到步骤107。
[0055] 在步骤107,将裂缝系统作为基岩系统的源汇项，将两者耦合起来建立嵌入式离散 裂缝流动数学模型。
[0056]
(5)
[0057] 式中，
上述方程的系数矩阵表达式具体如下：[0058]
(6)
[0059] 其中，凡为网格单元总数山=E"。从方程（5)可知：方程（6)中的各系数矩阵仅 与网格单元的几何信息和油藏参数有关，而对网格的几何形状没有特殊要求，求解方便，原 则上适用于任何复杂网格。
【主权项】
1. 基于嵌入式离散裂缝模型的压裂井建模及模拟方法，其特征在于，该基于嵌入式离 散裂缝模型的压裂井建模及模拟方法包括： 步骤1，划分基岩系统正交粗网格并提取粗网格节点及裂缝信息； 步骤2,将系统分为基岩系统和裂缝系统，并对其分别建立流动数学模型； 步骤3,通过流量等效原则对裂缝进行降维处理，建立相应的嵌入式离散裂缝网络模 型； 步骤4,基于模拟有限差分法建立基岩系统数值计算格式，将裂缝作为源汇项进行处 理； 步骤5,分析裂缝在基岩中的H类不同嵌入方式下数值计算格式； 步骤6,基于控制体积有限差分法建立嵌入式离散裂缝网络系统的数值计算格式；W 及 步骤7,将裂缝系统作为基岩系统的源汇项，将两者禪合起来建立嵌入式离散裂缝流动 数学模型。2. 根据权利要求1所述的基于嵌入式离散裂缝模型的压裂井建模及模拟方法，其特征 在于，在步骤2中，建立的基岩系统流动模型为：(1) 其中，Vb为网格块体积；4为孔隙度；！为绝对渗透率；kf,为油水相的相对渗透率，j =1为油化j = 2为水相；ii j为油水相的粘度；Pj为油水相压力值；Y j为油水相密度；D 为高度。3. 根据权利要求2所述的基于嵌入式离散裂缝模型的压裂井建模及模拟方法，其特征 在于，在步骤2中，建立的裂缝系统流动模型为：(2) 其中，为连接点的质量交换量。4. 根据权利要求1所述的基于嵌入式离散裂缝模型的压裂井建模及模拟方法，其特征 在于，在步骤4中，建立的基岩系统数值计算格式为：(3) 其中，Qi为窜流量；TIi为传导率；A 1为窜流系统;A为裂缝面积；k为渗透率；<d>为裂 缝到网格的平均距离；4为势函数的大小。5. 根据权利要求1所述的基于嵌入式离散裂缝模型的压裂井建模及模拟方法，其特征 在于，在步骤4中，H类不同嵌入方式包括完全穿透型、部分穿透型和交叉裂缝型，裂缝在 基岩中的H类不同嵌入方式下数值计算格式为： 完全穿透时 部分穿透时交叉裂缝型k"M为裂缝渗透率；CfM为两条裂缝间距离；a"m为裂缝开度乘W裂缝长度. V为网格大小；X。为裂缝到网格边界的距离；kfi，kf2为裂缝渗透率；《 fi，《 f2为裂缝开度； Lmt为网格内裂缝长度；dfi，df2为裂缝中必到网格边界的距离。6. 根据权利要求1所述的基于嵌入式离散裂缝模型的压裂井建模及模拟方法，其特征 在于，在步骤6中，基于控制体积有限差分法建立嵌入式离散裂缝网络系统的数值计算格 式为：(4) 其中：；ffi二VfiQfi ;Pfi为裂缝节点处的压力值；S为沿裂缝 方向的局部坐标系；T为传导率；Vfi为裂缝孔隙体积；Qmfi ;qffi为基岩与裂缝、裂缝与裂缝之间 一 C (1 ，裂缝单元与另一裂缝单元相交 W1〇，裂缝单元不与其它裂缝单元相交。7. 根据权利要求1所述的基于嵌入式离散裂缝模型的压裂井建模及模拟方法，其特征 在于，在步骤7中，建立的嵌入式离散裂缝流动数学模型为：巧、 式中，= [乂i- ] :; = [/，L ];巧》, = [/J?。； fm =[。，其中乂，,,? = I。,恥,'孤；上述方 程的系数矩阵表达式具体如下：其中，化为网格单元总数；Ii = Em,从方程妨可知：方程（6)中的各系数矩阵仅与网 格单元的几何信息和油藏参数有关，而对网格的几何形状没有特殊要求，适用于任何复杂 网格。
【专利摘要】本发明提供一种基于嵌入式离散裂缝模型的压裂井建模及模拟方法，该方法包括：划分基岩系统正交粗网格并提取粗网格节点及裂缝信息；将系统分为基岩系统和裂缝系统，并对其分别建立流动数学模型；建立相应的嵌入式离散裂缝网络模型；基于模拟有限差分法建立基岩系统数值计算格式；分析裂缝在基岩中的三类不同嵌入方式下数值计算格式；基于控制体积有限差分法建立嵌入式离散裂缝网络系统的数值计算格式；以及将裂缝系统作为基岩系统的源汇项，将两者耦合起来建立嵌入式离散裂缝流动数学模型。该基于嵌入式离散裂缝模型的压裂井建模及模拟方法避免了复杂的非结构化网格剖分过程，提高了压裂井裂缝描述和模拟方法的准确性和高效性。
【IPC分类】G06F17/50
【公开号】CN105653746
【申请号】
【发明人】杨勇, 李友全, 孙业恒, 张奎祥, 于伟杰, 高安邦, 阎燕, 王杰, 李弘博, 韩凤蕊, 郭建福, 黄成江
【申请人】中国石油化工股份有限公司, 中国石油化工股份有限公司胜利油田分公司地质科学研究院
【公开日】2016年6月8日
【申请日】2014年11月10日