4. 净压力拟合的数学模型框架

净压力拟合的本质，是用井口或井底压力反演裂缝的几何形态与延伸行为。要做到这一点，必须先把"压力—缝宽—缝长—滤失"之间的关系写成可计算的方程，这就是裂缝扩展模型。工程上常用的四类模型——PKN、KGD、径向（Penny-shaped）与拟三维（P3D）——并不是互相替代的"精度等级"，而是针对不同裂缝几何假设建立的近似解。选错模型，拟合曲线即使贴合，反演出的缝高、缝长也可能是错的。因此本节的重点不是背公式，而是搞清楚每个模型的假设边界，以及这些假设在净压力曲线上留下什么"指纹"。

一、四类模型的基本假设与几何图像

所有模型共享同一套力学骨架：裂缝内净压力 p_net 与缝宽 w 通过弹性方程联系，物质平衡方程约束注入体积等于裂缝体积加滤失体积，缝端满足断裂准则。差别在于对裂缝形状和缝高行为的简化方式。

PKN 模型假设裂缝高度恒定且受上下隔层限制，缝长远大于缝高，垂直于缝长方向的每个横截面独立变形，水平面内呈椭圆断面，缝内流体沿缝长一维流动。其关键推论是：缝宽与缝高成正比、与缝长关系较弱，净压力随缝长增长而缓慢上升，典型净压力—时间曲线呈平缓上翘。

KGD 模型同样假设缝高恒定，但认为裂缝在水平面内近似矩形、缝宽在缝长方向变化，且缝宽由垂直于裂缝面的平面应变控制。它适用于缝长与缝高相当、甚至缝高大于缝长的情形，净压力随时间上升更快，对滤失系数更敏感。

径向模型假设均质无限大地层中形成圆形（饼状）裂缝，缝高不受隔层约束、自由向上向下对称扩展，缝内为径向流。它常用于小型测试压裂、无遮挡薄层或校准孔眼摩阻的阶段，净压力曲线在早期表现为较陡的上升段。

拟三维模型放松了"缝高恒定"这一最强假设，把裂缝沿缝长和缝高两个方向离散，用二维弹性方程（或等效的缝高—缝宽耦合关系）计算缝宽，同时引入上下隔层的地应力差与断裂韧性，允许缝高在延伸中动态变化。它比前三者更接近真实，但输入参数更多，对地应力剖面和隔层性质的准确性高度依赖。

二、假设差异如何映射到净压力曲线

净压力拟合的判读逻辑，正是把实测曲线形态与上述"指纹"对照。若实测净压力在停泵前呈长时间平缓上翘、斜率稳定，且压后微地震或缝高监测显示缝高基本被隔层锁住，PKN 是合理起点；若净压力上升明显偏快、对前置液量敏感，同时施工压力绝对值不高，则更接近 KGD 描述的短宽缝。测试压裂阶段若净压力在几分钟内快速抬升后趋稳，径向模型常被用来解释近井筒的饼状缝或射孔段附近的初始扩展。

最需要警惕的是缝高失控。当实测净压力曲线出现明显的斜率转折甚至台阶式下降（净压力降低往往意味着裂缝进入低应力层、缝高突然扩张、缝内压力被"摊薄"），此时继续用恒定缝高的 PKN 或 KGD 拟合，会系统性地高估缝长、低估缝高。拟三维模型的价值就在这里：它把上下隔层的地应力差 Δσ 作为输入，当缝内净压力接近或超过 Δσ 时，缝高开始穿层，缝宽重新分配，净压力响应自然出现转折。因此，P3D 不是"更高级所以更好"，而是当存在缝高穿层证据时，它是唯一能自洽解释曲线转折的框架。

三、模型选择与拟合的实操路径

实际拟合时，建议按以下顺序推进，而不是一上来就套最复杂的模型：

1. 先做诊断性判断。用测试压裂或小型注入段确定闭合压力、瞬时停泵压力（ISIP）与孔眼摩阻，得到较干净的净压力序列；观察净压力随时间的变化形态，初步判断属于平缓型、快速上升型还是带转折型。
2. 用简单模型定基线。先以 PKN 或径向模型拟合早期段，确定滤失系数与缝宽量级；若拟合残差集中在后期，说明几何假设在后期失效，而不是滤失参数错了。
3. 引入 P3D 处理转折。把隔层地应力差、层厚、断裂韧性作为输入，允许缝高变化，重新拟合全时段曲线；重点检查转折点出现的时间是否与净压力台阶对应。
4. 用独立信息校验。拟合得到的缝长、缝高应与微地震事件分布、压后井温测井或示踪剂结果交叉验证，避免"曲线拟合得好、几何反演得错"。

四、参数敏感性：哪些输入最影响结论

四类模型对输入的敏感程度并不相同。PKN 与 KGD 对滤失系数和流体黏度较敏感，因为缝长主要由物质平衡决定；径向模型对注入排量和地层渗透率更敏感；P3D 则对隔层地应力差、层间弹性模量差异和断裂韧性最敏感——Δσ 估计偏差 1～2 MPa，就可能导致缝高穿透与否的判断完全反转。因此做 P3D 拟合时，应优先用测井或小型压裂数据约束地应力剖面，而不是把 Δσ 当作自由拟合参数随意调整。下表给出一个示例性的参数敏感性对照，数值为教学假设，仅用于说明相对影响方向。

需要强调的是，模型只是把压力翻译成几何的"词典"，词典选错，翻译再流畅也无意义。工程实践中更稳妥的做法是：先用简单模型快速给出量级判断，再用 P3D 处理复杂段，最后用独立监测数据闭环校验。净压力拟合的可靠性，最终取决于假设与地层的匹配程度，而不是方程的复杂程度。