4. 地应力场与天然裂缝表征

体积压裂能否形成复杂缝网,很大程度上取决于两个"先天条件":地应力场如何分布,天然裂缝如何展布。地应力决定了水力裂缝的起裂方向与延伸趋势,天然裂缝则提供了缝网转向、分叉与沟通的天然"薄弱面"。本节的任务,是把这两类信息从零散的测量数据,逐步整合为可用于压裂设计的三维耦合模型。

一、地应力大小与方向的确定方法

地应力描述需要四个基本量:最大水平主应力 σH、最小水平主应力 σh、垂向应力 σv,以及最大水平主应力的方位角。工程上通常约定:方位角以正北为 0°、顺时针为正;σH ≥ σh,二者均位于水平面内。全文与图示均采用这一约定。

不同测量手段的"探测尺度"差别很大,必须组合使用:

  • 测井方法:成像测井(FMI/CMI)通过井壁崩落椭圆长轴方向指示 σh 方位,钻井诱导缝走向则平行于 σH;交叉偶极声波测井测量快慢横波速度各向异性,快波偏振方向对应 σH 方位,各向异性幅度可半定量反映应力差。测井覆盖全井段,是方向信息的主力来源。
  • 岩心方法:古地磁定向岩心可做差应变分析(DSA)与波速各向异性测试,直接给出主应力方向;声发射 Kaiser 效应通过加载记忆点估算古应力与现应力大小。岩心法精度高但离散、成本高,适合标定。
  • 微地震与现场测试:压裂微地震事件云的长轴展布方向通常平行于 σH,可反演区域应力方位;小型压裂的瞬时停泵压力与闭合压力可约束 σh 大小,DFIT 测试则给出孔隙压力与闭合应力。

大小计算上,垂向应力由密度测井积分得到:σv = ∫ρ(z)g dz。水平应力常用多孔弹性模型估算,例如:

σh = (ν/(1-ν))·(σv - α·pp) + α·pp + Δσ_构造
σH = σh + Δσ_各向异性(由横波各向异性或DSA标定)

其中 ν 为泊松比,α 为 Biot 系数,pp 为孔隙压力,Δσ_构造 为构造附加应力。示例:某井 2500 m 处 σv=52 MPa、pp=25 MPa、ν=0.23、α=0.9,则弹性部分 σh≈(0.23/0.77)×(52−22.5)+22.5≈31.3 MPa,若横波各向异性标定应力差为 6 MPa,则 σH≈37.3 MPa。该结果需用 DFIT 闭合压力校核,误差通常控制在 5%~10%。

二、天然裂缝分布表征技术

天然裂缝表征的核心是回答三个问题:在哪、多密、什么方向。成像测井能直接识别井壁上的裂缝产状,统计出线密度(条/m)与走向玫瑰图;岩心观察可补充充填程度与开度;地震属性(曲率、相干、蚂蚁追踪)则把井点信息外推到三维空间,刻画裂缝带的平面展布。

需要强调的是,井点裂缝密度不能简单按算术平均外推。更合理的做法是建立裂缝密度与构造曲率、断层距离、岩性组合的关系,用地震属性作为约束面进行协同克里金插值。离散裂缝网络(DFN)建模时,常用"分形维数 + 方位分组"的方式生成裂缝面:先按走向分成 2~3 组,每组给定密度、长度分布(如幂律指数)与开度,再叠加到三维网格中。

示例参数:某页岩段成像测井解释出 3 组裂缝,走向分别为 NE45°、NW310°、近 EW 向,线密度 0.8、0.5、0.3 条/m,充填比例约 60%。DFN 建模时对充填缝降低开度权重,仅保留约 40% 作为潜在渗流通道。

三、地应力场与天然裂缝的耦合建模

把两者"耦合"起来,本质是回答:在给定应力状态下,哪些天然裂缝会张开、哪些会剪切滑动。这直接决定压裂液能否沿天然缝网络扩散。判断依据是莫尔-库仑准则与裂缝面正应力:

裂缝面正应力:σn = σH·cos²θ + σh·sin²θ
(θ 为裂缝走向与 σH 方位夹角)
张开条件:σn < pp(或有效正应力趋零)
剪切滑动条件:τ ≥ C + μ·(σn − pp)

由此可得一条重要工程结论:当天然裂缝走向与 σH 夹角较小(约 0°~30°)时,裂缝面正应力低、易张开;夹角接近 90° 时正应力高,多表现为闭合或仅能剪切。这正是"应力阴影"与缝网转向设计的物理基础。

耦合建模的落地流程通常是:以测井与岩心确定的一维应力剖面为约束,结合地震解释的断层与裂缝带,建立三维非均质应力场(含断层附近的应力扰动);再将 DFN 嵌入该应力场,逐条裂缝计算 σn 与滑动趋势系数,输出"可张开裂缝"与"可滑动裂缝"两套网络。压裂设计时,射孔簇方位应尽量使水力裂缝延伸方向与易张开天然缝组形成有利夹角,从而提升缝网复杂度。

地应力—天然裂缝耦合建模流程与判据 应力数据来源 • 成像测井:崩落/诱导缝 • 偶极声波:快慢波各向异性 • 岩心 DSA / Kaiser 效应 • DFIT / 微地震事件云 • 密度测井积分 σv 裂缝数据来源 • 成像测井裂缝产状 • 岩心观察:充填/开度 • 地震曲率/蚂蚁追踪 • 走向分组与密度统计 三维应力场建模 σv 积分 + 多孔弹性 σh 横波各向异性标定 σH 断层附近应力扰动修正 DFN 离散裂缝网络 走向分组 + 密度 + 长度幂律 充填缝降权,保留有效通道 嵌入三维网格 裂缝稳定性判据 σn = σH·cos²θ + σh·sin²θ 张开:σn < pp 滑动:τ ≥ C + μ(σn − pp) θ 小(0°~30°)→ 易张开 θ 大(≈90°)→ 闭合/剪切 输出:可张开缝网 + 可滑动缝网 → 指导射孔方位与簇间距 → 评估缝网复杂度潜力 压裂设计反馈 射孔簇方位对准易张开缝组 · 簇间距避开高应力遮挡 · 用微地震事件云校核缝网方向 若事件云长轴偏离 σH 预测方向,回查应力方位或 DFN 分组是否需修正

上图把"数据来源—建模—判据—反馈"串成一条闭环。实际项目中,最容易出问题的环节往往不是计算本身,而是坐标系与方位约定不统一:测井解释的裂缝走向、地震属性的展布方向、微地震事件云的长轴,若参照系不一致,耦合结果会系统性偏移。建议在建模前统一建立"正北为 0°、顺时针"的方位基准,并在成果表中同时记录 σH 方位与各裂缝组走向,便于交叉校核。

注意:多孔弹性模型估算的水平应力对泊松比与孔隙压力非常敏感,ν 变化 0.02 即可引起 σh 约 1~2 MPa 的偏差。务必用 DFIT 闭合压力或小型压裂数据反演校正,不可直接套用区域经验值。

完成耦合建模后,应输出一套可直接用于压裂设计的成果:三维应力场网格(含 σH 方位、σh、σv)、分组 DFN、以及按"易张开/易滑动/闭合"分类的裂缝面清单。设计人员据此确定射孔方位、簇间距与泵注规模,再用压后微地震事件云的长轴方向与展布范围做验证——若事件云明显偏离预测的 σH 方向,说明应力方位或天然缝分组需要复核。至此,地应力与天然裂缝从各自独立的测量数据,转化为服务于缝网设计的统一约束条件。