页岩气压裂数值模型分析
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
ＤＯＩ ：１０ ．３７８７／ｊ ．ｉｓｓｎ ．１０００‐０９７６ ．２０１１ ．１２ ．０１４
１ 页岩气藏的特点
１ ．１ 特殊的赋存生产机理 页岩既是烃源岩又是储集层 ，就近赋存是页岩气 成藏的特点 。 页岩气的赋存方式多样 ，游离方式 、吸附 状态和溶解状态并存 。 总体上主要以游离气和吸附气 为主 ，吸附状态天然气的含量变化介于 ２０％ ～ ８５％ 。 目前认为页岩气的产出分为 ３ 个阶段 ：① 在压降 的作用下 ，基质系统中的页岩气在基质表面进行解吸 附 ；② 在浓度差的作用下 ，页岩气由基质系统向裂缝系 统进行扩散 ；③ 在流动势的作用下 ，页岩气通过裂缝系 统流向生产井筒 。 由于裂缝空间的有限性 ，因此早期 以游离气为主的天然气产量快速下降并且达到稳定 ， 稳定期的产量主要是基质孔隙里的游离气和解吸气 。 水力压裂可以增大裂缝空间和连通性 ，使更多的吸附 气发生解吸附而向裂缝聚集 。 １ ．２ 纳米级的微观孔隙结构 通过扫描电镜等成像技术和脉冲法等测试技术研 究表明 ，纳米级的有机质孔隙是页岩的主要储集空 间［１］ 和孔隙类型 ，其形成与分布与有机质的丰度密切 相关 ；岩心观察表明天然裂缝较为发育 ，但绝大部分被 矿物充填处于闭合状态 。 孔隙和吼道的尺寸为纳米级
网络出版时间：2011-12-24 21:30 网络出版地址：http://www.cnki.net/kcms/detail/51.1179.TE.20111224.2130.009.html
第 ３１ 卷第 １２ 期 开 发 工 程
·２·
天 然 气 工 业 ２０１１ 年 １２ 月
图 １ 页岩水平井压裂微地震监测结果图
（注 ：１ ｆ ｔ ＝ ０ ．３０４ ８ ｍ ，下同）
纳米级的微孔隙和吼道 ，微米级的大孔隙或天然裂缝 以及厘米级别的人工裂缝等 ，不能用唯一的流动状态 进行描述 。 可采用 Ｋｎｕｄｓｅｎ 系数进行流态的划分 ，在 不同的流动状态下 ，对应以不同的渗流方程进行描述 。 页岩气体的流动状态主要包括 ：在有机质表面的解吸 附 ，可采用 Ｌａｎｇｍｕｉｒ 等温吸附模型进行表征 ；在纳米 级孔隙中的自由分析运动状态 ，近似地认为气体遵循 扩散定律运动 ；在微裂缝中可能出现滑脱现象 ；在人工 裂缝中遵循线性渗流规律 。 其中在纳米级孔隙中的微 观渗流机理目前尚不明确 ，有学者采用流体力学格子 Ｂｏｌｔｚｍａｎｎ 方法结合数字成像技术从介观层次的角度 进行研究 ，提出了一个新的研究方向 。
·１·
页岩气压裂数值模型分析
张士诚１ 牟松茹１ 崔 勇２
１ ．中国石油大学（北京）石油工程教育部重点实验室 ２ ．中国石油海外勘探开发公司
张士诚等 ．页岩气压裂数值模型分析 ．天然气工业 ，２０１１ ，３１（１２） ：８１‐８４ ． 摘 要 水力压裂和水平井开采是页岩气开发的主要技术 ，在我国尚处在工业试验阶段 ，存在很多技术瓶颈 。 在总 结分析了页岩气压裂的特点基础上 ，探讨了网状裂缝形成的主控因素及裂缝扩展模型 、产能预测模型的类型以及优缺 点 。 结果认为 ，特殊的赋存生产机理 、复杂的裂缝形态和多尺度的渗流模式是页岩气压裂的主要特点 ，其目的是形成网 状裂缝 ，扩大储层改造体积 ；网状裂缝的形成主要受天然裂缝与人工裂缝的夹角 、水平主应力差和岩石的脆性等因素的 控制 。 页岩气压裂产能预测模型面临的主要问题是裂缝形态的模拟和气体流态的描述 ，主要有非常规裂缝模型 、离散裂 缝模型和双重介质模型等 ，这些模型和方法在一定程度上表征了页岩气压裂裂缝形态和渗流特点 ，但没有考虑不规则的 裂缝形态等 。 关键词 页岩气 开发 压裂（岩石） 裂缝扩展模型 产能预测模型 渗流 特点
由微地震监测结果确定 ，仅限于本段压裂施工模拟 ，计
算结果不具有普遍适用性 。
线网模型为半解析模型 （图 ２ ） ，泄流区域呈椭圆
形 ，有两组垂直正交的裂缝分别平行于椭圆平面上两
个主应力（σｈ 、σＨ ） ，为避免裂缝间距的主观性 ，ｄ ｘ 、ｄ ｙ
依据实际测试数据确定 ，椭圆的长短轴 ａ 和 ｂ 以及高
２ 网状裂缝的形成及主控因素分析
Ｊｏｎ 等［３］ 采用边界元法对压裂时多裂缝的同时延
伸和它们与天然裂缝之间的相互作用进行了研究 ，认
为在天然裂缝发育的条件下 ，天然裂缝与人工裂缝的
夹角 、拟净压力wenku.baidu.com数是影响网状裂缝形成的主要因素 ，
水平井中人工裂缝和天然裂缝夹角越大 ，值越大则越
容易形成网状裂缝 ；直井条件下不容易形成网状裂缝 。
·３·
２π ｙ ｒ
抄φ 抄ｔ
－
２
抄 抄ｙ
Ｂ（１ ＋ ｒ）ｙｗ ｙ Ｋ ｆ ｙ 抄 ｐ
ｒ２ μｄ ｙ
抄ｙ
＝０
（４ ）
式中 ｑ 为压裂液量 ，ｍ３ ／ｄ ；ｔｉ 为压裂液排量 ，ｍ３ ／ｍｉｎ ；
Ｌｘｉ 、Ｌｙｊ 分别为 ｘ 和 ｙ 方向裂缝段长度 ，ｍ ；Ｗ ｘｉ 、Ｗ ｙｉ 分别
为 ｘ 和 ｙ 方向裂缝段宽度 ，ｍ ；Ｋｆ ｘ 、Ｋｆ ｙ 分别为 ｘ 和 ｙ 方
于更好地了解天然裂缝对于预想的人工裂缝几何形状
和对产能的影响 ，模型的计算结果可以通过微地震监
测进行校正 。 非常规模型的主要问题在于 ：① 天然裂
缝的分布依赖于离散裂缝地质建模的结果 ；② 非常规
复杂裂缝模型对输入参数的精确性要求较高 。
图 ３ 非常规复杂裂缝模型图
４ 产能预测模型
４ ．１ 离散裂缝模型
椭球体的实时扩展 ，考虑了施工参数的影响 ，并计算了
支撑剂在裂缝中的分布情况 。 线网的不足之处在于 ：
① 它必须将油藏改造区域近似为沿井筒对称的椭球
体 ，不能模拟不规则的裂缝形态 ；② 没有建立判断准
则 ，直接地认为天然裂缝与人工裂缝相连接 ；③ 没有考
虑人工裂缝之间的相互干扰 ；④ 裂缝间距和改造体积
别 ，孔隙 、吼道配置关系复杂 ；基质渗透率为纳达西级 别 ，孔隙度一般小于 ７ ％ 。 页岩储层纳米级的微观孔隙结构 ，与相同孔隙度 的微米级孔隙相比提供了更大的比表面积 ，为气体的 吸附提供了条件 。 但是也相应引发如下的问题 ：① 纳 米级储层的物性特征参数难以用常规的方法测量和计 算 ；② 气体在纳米级孔隙中的渗流复杂多变 ，流动规律 目前尚不明确 ；③ 需采用如水平井多级压裂等特殊的 开发方式才能获得经济产量 ，且增产的机理也与常规 压裂不同 。 １ ．３ 水力压裂形成复杂的裂缝形态 常规压裂形成的裂缝一般呈双翼对称裂缝的形 式 。 但页岩气压裂中微地震监测的结果表明 ，裂缝的 形态复杂多变 ，如图 １ 所示［２］ 。 １ 口水平井压裂后微 地震监测结果表明 ，第 １ 、２ 段压裂施工形成了垂直于 水平井段的平面缝 ，第 ３ 、４ 段施工形成了网状裂缝 。 目前认为页岩气压裂目的 ，就是要建立一个独立于传 统意义裂缝半长的更加庞大的裂缝网络系统 ，实现更 大规模的储层改造波及体积 ；生产实践也证明了储层 改造体积越大压后增产效果越好 。 １ ．４ 多尺度的流动状态 页岩储层压裂后形成了多尺度的流动空间 ，包括
模拟网状裂缝面 ，裂缝在空间上呈三维分布互相正交 ；
考虑了稳态和拟稳态流动的情况 ，计算出了产能预测
图 ４ 产能预测数值模型图
图版 。 Ｃｉｐｏｌｌａ 等［９］ 建立了等距分布正交分布的网状 裂缝产能预测模型 ，采用 ＳＲＶ 方法计算改造体积 。 模 拟计算中考虑了支撑剂的分布形态 （支撑裂缝和未支 撑裂缝） ，渗流模式（解吸附 、达西渗流和非达西渗流 ） 以及裂缝参数对产能的影响 。 Ｗｉｌｌｉａｍｓ 等［１０］ 通过裂 缝识别技术确认裂缝面和裂缝数量 ；结合成像测井确 定的裂缝方位 ，建立裂缝分布模型后简化为气藏数值 模拟模型 ，通过历史拟合进行模型校正和参数优化 。 ４ ．２ 双重介质模型 双重介质模型建立的基础是 Ｗａｒｒｅｎ 和 Ｒｏｏｔ 提 出的双重孔隙介质模型 ，该模型强调了裂缝性油藏的 双孔隙的本质 ，简化了裂缝性油藏中的连通性和与模 型规模有关的非均质性问题 。 Ｚｈａｎｇ 等［１１］ 采用了微地震解释获得的潜在裂缝 分布状态结果 ，得到模拟主要横切裂缝的双重介质油 藏模型 。 计算中考虑了裂缝间距 、非达西效应等对产 能的影响 ，表明产能对窜流因子参数最为敏感 。 Ｃｈａｎｇａｎ 等［１２］ 建立了离散裂缝网络模型 ，通过 ＥＳＶ 方法确定压裂有效体积 ，作为进一步约束双重介 质参数的条件 。 该方法利用岩石力学参数和施工参数 估算平均裂缝宽度 ，建立裂缝强度 、裂缝宽度和施工规 模的关系式 ，通过和微地震监测结果对比 ，形成超大规 模的双重介质模拟 。 模拟计算中考虑了等温吸附 、不 稳定扩散和窜流等问题 。 可作为创新点
ｐ － Ｓ ｆｒａｃ
ｈ ，ｍｉｎ
Ｒ ＝ ｎ
Ｓ － Ｓ ｈ ，ｍａｘ
ｈ ，ｍｉｎ
（１ ）
式中 Ｒｎ 为拟净压力系数 ；ｐｆｒａｃ 为净压力 ，Ｍ Ｐａ ；Ｓｈ ，ｍａｘ 为水 平最大主应力 ，ＭＰａ ；Ｓｈ ，ｍｉｎ 为水平最小主应力 ，ＭＰａ 。
Ｇｕ 等［４］ 建立了天然裂缝和人工裂缝相互作用的
离散裂缝模型建立的基础是 Ｆｉｓｈｅｒ 等［７］ 提出的
复杂裂缝形态理论 ，认为页岩气压裂后形成了复杂的
裂缝形态 ，简化为多裂缝或交错分布的形态 ；离散裂缝
模型包括三维的线网模型 、二维的离散裂缝模型以及
随机分布的多裂缝模型（图 ４） 。
重点
Ｂｒｕｃｅ 等［８］ 建立了线网模型 ，采用正交的椭球面
度 ｈ 由微地震结果给出 。
求解方程如下 ：
Ｎｘ
Ｎｙ
∑ ∑ ｑｔｉ － ｈ
Ｌ ｘｉ Ｗ珨 ｘｉ ＋
Ｌ ｙｊ Ｗ珨 ｙｊ ＝ ０
（２ ）
ｉ＝ １
ｊ＝１
２π ｒｘ
抄φ 抄ｔ
－
２
抄 抄ｘ
Ｂ（１ ＋ ｒ）ｘｗ ｘ Ｋ ｆ ｘ 抄 ｐ
μｄ ｘ
抄ｘ
＝０
（３ ）
图 ２ 线网模型图
第 ３１ 卷第 １２ 期 开 发 工 程
判断准则 ，考虑了非正交的模式 ，采用 Ｕ ＦＭ 模型进行
计算 ，结果表明水平面主应力差 、缝内净压力 、天然裂
缝密度以及岩性是影响裂缝形态的主要因素 ；水平主
应力差越小 、天然裂缝发育程度越高 ，易于形成网状裂
缝 ；杨氏模量越大 ，岩石脆性越强 ，易于形成粗糙节理 并保持裂缝开启 ；天然裂缝密度 、基质渗透率等参数也 对裂缝的形态和规模有一定的影响 。 生产实践表明 ，页岩中石英的含量越多岩石的脆 性越大 ，越容易形成网状裂缝 ，获得较好的增产效果 。 提出了脆度指数的概念 ，即岩石中石英的含量占全部 矿物成分的百分比 。 在施工设计中 ，脆度指数越大（大 于 ５０％ ） ，则采用清水压裂 ，大液量 、少支撑剂量的方 法 ，形成网状裂缝获得较好的效果 ；若脆度指数小 （低 于 ３０％ ） ，则采用常规压裂的方法 。
动 、支撑剂的输送和裂缝宽度的弹性变形 ，通过计算应
力阴影考虑相邻人工裂缝间的相互作用 ；采用三层模
型模拟支撑剂输送 ，从底至上依次为支撑剂 ，混砂液和
压裂液 。 非常规模型比线网模型提供了一种更为精确
预测裂缝分布 、几何形态和支撑剂分布的方法 ，充分考
虑了储层岩石力学性质以及不规则的裂缝形态 ，有助
３ 复杂裂缝模型
３ ．１ 线网模型
Ｘｕ 等［５］ 提出了表示复杂裂缝形 态的线网模型
（ｗ ｉｒｅ‐ｍｅｓｈ） ，认为页岩气藏水平井压裂产生的裂缝网
络是沿井筒对称的椭球体 ，通过将该椭球体划分为数
条正交的水平 、垂直均匀截面来描述高渗裂缝 。 线网
模型计算时应用岩石力学方法考虑了压裂过程中裂缝
作者简介 ：张士诚 ，１９６３ 年生 ，教授 ，博士生导师 ，本刊第七届编委会委员 ；长期从事采油工程理论与技术 、油气渗流理论与应用 的教学与研究工作 。 地址 ：（１０２２４９）北京市昌平区府学路 １８ 号 。 电话 ：（０１０） ８９７３３０４７ 。 Ｅ‐ｍａｉｌ ：ｚｈａｎｇｓｃ＠ ｃｕｐ ．ｅｄｕ ．ｃｎ
向裂缝渗透率 ，ｍＤ ；ｒ 为椭圆的宽高比 ；Ｂ 为椭圆积分 。
３ ．２ 非常规模型
Ｗｅｎｇ 等［６］ 提出了非常规裂缝扩展模型 （Ｕ ＦＭ‐
ｕｎｃｏｎｖｅｎｔｉｏｎａｌ ｆｒａｃｔｕｒｅ ｍｏｄｅｌ） 。 该模型为数值模型
（图 ３） ，能够模拟天然裂缝和人工裂缝之间的相互作
用 ，建立了裂缝端部扩展准则 ，考虑了压裂液的一维流
１ 页岩气藏的特点
１ ．１ 特殊的赋存生产机理 页岩既是烃源岩又是储集层 ，就近赋存是页岩气 成藏的特点 。 页岩气的赋存方式多样 ，游离方式 、吸附 状态和溶解状态并存 。 总体上主要以游离气和吸附气 为主 ，吸附状态天然气的含量变化介于 ２０％ ～ ８５％ 。 目前认为页岩气的产出分为 ３ 个阶段 ：① 在压降 的作用下 ，基质系统中的页岩气在基质表面进行解吸 附 ；② 在浓度差的作用下 ，页岩气由基质系统向裂缝系 统进行扩散 ；③ 在流动势的作用下 ，页岩气通过裂缝系 统流向生产井筒 。 由于裂缝空间的有限性 ，因此早期 以游离气为主的天然气产量快速下降并且达到稳定 ， 稳定期的产量主要是基质孔隙里的游离气和解吸气 。 水力压裂可以增大裂缝空间和连通性 ，使更多的吸附 气发生解吸附而向裂缝聚集 。 １ ．２ 纳米级的微观孔隙结构 通过扫描电镜等成像技术和脉冲法等测试技术研 究表明 ，纳米级的有机质孔隙是页岩的主要储集空 间［１］ 和孔隙类型 ，其形成与分布与有机质的丰度密切 相关 ；岩心观察表明天然裂缝较为发育 ，但绝大部分被 矿物充填处于闭合状态 。 孔隙和吼道的尺寸为纳米级
网络出版时间：2011-12-24 21:30 网络出版地址：http://www.cnki.net/kcms/detail/51.1179.TE.20111224.2130.009.html
第 ３１ 卷第 １２ 期 开 发 工 程
·２·
天 然 气 工 业 ２０１１ 年 １２ 月
图 １ 页岩水平井压裂微地震监测结果图
（注 ：１ ｆ ｔ ＝ ０ ．３０４ ８ ｍ ，下同）
纳米级的微孔隙和吼道 ，微米级的大孔隙或天然裂缝 以及厘米级别的人工裂缝等 ，不能用唯一的流动状态 进行描述 。 可采用 Ｋｎｕｄｓｅｎ 系数进行流态的划分 ，在 不同的流动状态下 ，对应以不同的渗流方程进行描述 。 页岩气体的流动状态主要包括 ：在有机质表面的解吸 附 ，可采用 Ｌａｎｇｍｕｉｒ 等温吸附模型进行表征 ；在纳米 级孔隙中的自由分析运动状态 ，近似地认为气体遵循 扩散定律运动 ；在微裂缝中可能出现滑脱现象 ；在人工 裂缝中遵循线性渗流规律 。 其中在纳米级孔隙中的微 观渗流机理目前尚不明确 ，有学者采用流体力学格子 Ｂｏｌｔｚｍａｎｎ 方法结合数字成像技术从介观层次的角度 进行研究 ，提出了一个新的研究方向 。
·１·
页岩气压裂数值模型分析
张士诚１ 牟松茹１ 崔 勇２
１ ．中国石油大学（北京）石油工程教育部重点实验室 ２ ．中国石油海外勘探开发公司
张士诚等 ．页岩气压裂数值模型分析 ．天然气工业 ，２０１１ ，３１（１２） ：８１‐８４ ． 摘 要 水力压裂和水平井开采是页岩气开发的主要技术 ，在我国尚处在工业试验阶段 ，存在很多技术瓶颈 。 在总 结分析了页岩气压裂的特点基础上 ，探讨了网状裂缝形成的主控因素及裂缝扩展模型 、产能预测模型的类型以及优缺 点 。 结果认为 ，特殊的赋存生产机理 、复杂的裂缝形态和多尺度的渗流模式是页岩气压裂的主要特点 ，其目的是形成网 状裂缝 ，扩大储层改造体积 ；网状裂缝的形成主要受天然裂缝与人工裂缝的夹角 、水平主应力差和岩石的脆性等因素的 控制 。 页岩气压裂产能预测模型面临的主要问题是裂缝形态的模拟和气体流态的描述 ，主要有非常规裂缝模型 、离散裂 缝模型和双重介质模型等 ，这些模型和方法在一定程度上表征了页岩气压裂裂缝形态和渗流特点 ，但没有考虑不规则的 裂缝形态等 。 关键词 页岩气 开发 压裂（岩石） 裂缝扩展模型 产能预测模型 渗流 特点
由微地震监测结果确定 ，仅限于本段压裂施工模拟 ，计
算结果不具有普遍适用性 。
线网模型为半解析模型 （图 ２ ） ，泄流区域呈椭圆
形 ，有两组垂直正交的裂缝分别平行于椭圆平面上两
个主应力（σｈ 、σＨ ） ，为避免裂缝间距的主观性 ，ｄ ｘ 、ｄ ｙ
依据实际测试数据确定 ，椭圆的长短轴 ａ 和 ｂ 以及高
２ 网状裂缝的形成及主控因素分析
Ｊｏｎ 等［３］ 采用边界元法对压裂时多裂缝的同时延
伸和它们与天然裂缝之间的相互作用进行了研究 ，认
为在天然裂缝发育的条件下 ，天然裂缝与人工裂缝的
夹角 、拟净压力wenku.baidu.com数是影响网状裂缝形成的主要因素 ，
水平井中人工裂缝和天然裂缝夹角越大 ，值越大则越
容易形成网状裂缝 ；直井条件下不容易形成网状裂缝 。
·３·
２π ｙ ｒ
抄φ 抄ｔ
－
２
抄 抄ｙ
Ｂ（１ ＋ ｒ）ｙｗ ｙ Ｋ ｆ ｙ 抄 ｐ
ｒ２ μｄ ｙ
抄ｙ
＝０
（４ ）
式中 ｑ 为压裂液量 ，ｍ３ ／ｄ ；ｔｉ 为压裂液排量 ，ｍ３ ／ｍｉｎ ；
Ｌｘｉ 、Ｌｙｊ 分别为 ｘ 和 ｙ 方向裂缝段长度 ，ｍ ；Ｗ ｘｉ 、Ｗ ｙｉ 分别
为 ｘ 和 ｙ 方向裂缝段宽度 ，ｍ ；Ｋｆ ｘ 、Ｋｆ ｙ 分别为 ｘ 和 ｙ 方
于更好地了解天然裂缝对于预想的人工裂缝几何形状
和对产能的影响 ，模型的计算结果可以通过微地震监
测进行校正 。 非常规模型的主要问题在于 ：① 天然裂
缝的分布依赖于离散裂缝地质建模的结果 ；② 非常规
复杂裂缝模型对输入参数的精确性要求较高 。
图 ３ 非常规复杂裂缝模型图
４ 产能预测模型
４ ．１ 离散裂缝模型
椭球体的实时扩展 ，考虑了施工参数的影响 ，并计算了
支撑剂在裂缝中的分布情况 。 线网的不足之处在于 ：
① 它必须将油藏改造区域近似为沿井筒对称的椭球
体 ，不能模拟不规则的裂缝形态 ；② 没有建立判断准
则 ，直接地认为天然裂缝与人工裂缝相连接 ；③ 没有考
虑人工裂缝之间的相互干扰 ；④ 裂缝间距和改造体积
别 ，孔隙 、吼道配置关系复杂 ；基质渗透率为纳达西级 别 ，孔隙度一般小于 ７ ％ 。 页岩储层纳米级的微观孔隙结构 ，与相同孔隙度 的微米级孔隙相比提供了更大的比表面积 ，为气体的 吸附提供了条件 。 但是也相应引发如下的问题 ：① 纳 米级储层的物性特征参数难以用常规的方法测量和计 算 ；② 气体在纳米级孔隙中的渗流复杂多变 ，流动规律 目前尚不明确 ；③ 需采用如水平井多级压裂等特殊的 开发方式才能获得经济产量 ，且增产的机理也与常规 压裂不同 。 １ ．３ 水力压裂形成复杂的裂缝形态 常规压裂形成的裂缝一般呈双翼对称裂缝的形 式 。 但页岩气压裂中微地震监测的结果表明 ，裂缝的 形态复杂多变 ，如图 １ 所示［２］ 。 １ 口水平井压裂后微 地震监测结果表明 ，第 １ 、２ 段压裂施工形成了垂直于 水平井段的平面缝 ，第 ３ 、４ 段施工形成了网状裂缝 。 目前认为页岩气压裂目的 ，就是要建立一个独立于传 统意义裂缝半长的更加庞大的裂缝网络系统 ，实现更 大规模的储层改造波及体积 ；生产实践也证明了储层 改造体积越大压后增产效果越好 。 １ ．４ 多尺度的流动状态 页岩储层压裂后形成了多尺度的流动空间 ，包括
模拟网状裂缝面 ，裂缝在空间上呈三维分布互相正交 ；
考虑了稳态和拟稳态流动的情况 ，计算出了产能预测
图 ４ 产能预测数值模型图
图版 。 Ｃｉｐｏｌｌａ 等［９］ 建立了等距分布正交分布的网状 裂缝产能预测模型 ，采用 ＳＲＶ 方法计算改造体积 。 模 拟计算中考虑了支撑剂的分布形态 （支撑裂缝和未支 撑裂缝） ，渗流模式（解吸附 、达西渗流和非达西渗流 ） 以及裂缝参数对产能的影响 。 Ｗｉｌｌｉａｍｓ 等［１０］ 通过裂 缝识别技术确认裂缝面和裂缝数量 ；结合成像测井确 定的裂缝方位 ，建立裂缝分布模型后简化为气藏数值 模拟模型 ，通过历史拟合进行模型校正和参数优化 。 ４ ．２ 双重介质模型 双重介质模型建立的基础是 Ｗａｒｒｅｎ 和 Ｒｏｏｔ 提 出的双重孔隙介质模型 ，该模型强调了裂缝性油藏的 双孔隙的本质 ，简化了裂缝性油藏中的连通性和与模 型规模有关的非均质性问题 。 Ｚｈａｎｇ 等［１１］ 采用了微地震解释获得的潜在裂缝 分布状态结果 ，得到模拟主要横切裂缝的双重介质油 藏模型 。 计算中考虑了裂缝间距 、非达西效应等对产 能的影响 ，表明产能对窜流因子参数最为敏感 。 Ｃｈａｎｇａｎ 等［１２］ 建立了离散裂缝网络模型 ，通过 ＥＳＶ 方法确定压裂有效体积 ，作为进一步约束双重介 质参数的条件 。 该方法利用岩石力学参数和施工参数 估算平均裂缝宽度 ，建立裂缝强度 、裂缝宽度和施工规 模的关系式 ，通过和微地震监测结果对比 ，形成超大规 模的双重介质模拟 。 模拟计算中考虑了等温吸附 、不 稳定扩散和窜流等问题 。 可作为创新点
ｐ － Ｓ ｆｒａｃ
ｈ ，ｍｉｎ
Ｒ ＝ ｎ
Ｓ － Ｓ ｈ ，ｍａｘ
ｈ ，ｍｉｎ
（１ ）
式中 Ｒｎ 为拟净压力系数 ；ｐｆｒａｃ 为净压力 ，Ｍ Ｐａ ；Ｓｈ ，ｍａｘ 为水 平最大主应力 ，ＭＰａ ；Ｓｈ ，ｍｉｎ 为水平最小主应力 ，ＭＰａ 。
Ｇｕ 等［４］ 建立了天然裂缝和人工裂缝相互作用的
离散裂缝模型建立的基础是 Ｆｉｓｈｅｒ 等［７］ 提出的
复杂裂缝形态理论 ，认为页岩气压裂后形成了复杂的
裂缝形态 ，简化为多裂缝或交错分布的形态 ；离散裂缝
模型包括三维的线网模型 、二维的离散裂缝模型以及
随机分布的多裂缝模型（图 ４） 。
重点
Ｂｒｕｃｅ 等［８］ 建立了线网模型 ，采用正交的椭球面
度 ｈ 由微地震结果给出 。
求解方程如下 ：
Ｎｘ
Ｎｙ
∑ ∑ ｑｔｉ － ｈ
Ｌ ｘｉ Ｗ珨 ｘｉ ＋
Ｌ ｙｊ Ｗ珨 ｙｊ ＝ ０
（２ ）
ｉ＝ １
ｊ＝１
２π ｒｘ
抄φ 抄ｔ
－
２
抄 抄ｘ
Ｂ（１ ＋ ｒ）ｘｗ ｘ Ｋ ｆ ｘ 抄 ｐ
μｄ ｘ
抄ｘ
＝０
（３ ）
图 ２ 线网模型图
第 ３１ 卷第 １２ 期 开 发 工 程
判断准则 ，考虑了非正交的模式 ，采用 Ｕ ＦＭ 模型进行
计算 ，结果表明水平面主应力差 、缝内净压力 、天然裂
缝密度以及岩性是影响裂缝形态的主要因素 ；水平主
应力差越小 、天然裂缝发育程度越高 ，易于形成网状裂
缝 ；杨氏模量越大 ，岩石脆性越强 ，易于形成粗糙节理 并保持裂缝开启 ；天然裂缝密度 、基质渗透率等参数也 对裂缝的形态和规模有一定的影响 。 生产实践表明 ，页岩中石英的含量越多岩石的脆 性越大 ，越容易形成网状裂缝 ，获得较好的增产效果 。 提出了脆度指数的概念 ，即岩石中石英的含量占全部 矿物成分的百分比 。 在施工设计中 ，脆度指数越大（大 于 ５０％ ） ，则采用清水压裂 ，大液量 、少支撑剂量的方 法 ，形成网状裂缝获得较好的效果 ；若脆度指数小 （低 于 ３０％ ） ，则采用常规压裂的方法 。
动 、支撑剂的输送和裂缝宽度的弹性变形 ，通过计算应
力阴影考虑相邻人工裂缝间的相互作用 ；采用三层模
型模拟支撑剂输送 ，从底至上依次为支撑剂 ，混砂液和
压裂液 。 非常规模型比线网模型提供了一种更为精确
预测裂缝分布 、几何形态和支撑剂分布的方法 ，充分考
虑了储层岩石力学性质以及不规则的裂缝形态 ，有助
３ 复杂裂缝模型
３ ．１ 线网模型
Ｘｕ 等［５］ 提出了表示复杂裂缝形 态的线网模型
（ｗ ｉｒｅ‐ｍｅｓｈ） ，认为页岩气藏水平井压裂产生的裂缝网
络是沿井筒对称的椭球体 ，通过将该椭球体划分为数
条正交的水平 、垂直均匀截面来描述高渗裂缝 。 线网
模型计算时应用岩石力学方法考虑了压裂过程中裂缝
作者简介 ：张士诚 ，１９６３ 年生 ，教授 ，博士生导师 ，本刊第七届编委会委员 ；长期从事采油工程理论与技术 、油气渗流理论与应用 的教学与研究工作 。 地址 ：（１０２２４９）北京市昌平区府学路 １８ 号 。 电话 ：（０１０） ８９７３３０４７ 。 Ｅ‐ｍａｉｌ ：ｚｈａｎｇｓｃ＠ ｃｕｐ ．ｅｄｕ ．ｃｎ
向裂缝渗透率 ，ｍＤ ；ｒ 为椭圆的宽高比 ；Ｂ 为椭圆积分 。
３ ．２ 非常规模型
Ｗｅｎｇ 等［６］ 提出了非常规裂缝扩展模型 （Ｕ ＦＭ‐
ｕｎｃｏｎｖｅｎｔｉｏｎａｌ ｆｒａｃｔｕｒｅ ｍｏｄｅｌ） 。 该模型为数值模型
（图 ３） ，能够模拟天然裂缝和人工裂缝之间的相互作
用 ，建立了裂缝端部扩展准则 ，考虑了压裂液的一维流