搞水力压裂数值模拟的人很多时候不是在跟最终输出曲线较劲而是在跟“水力裂缝撞上天然裂缝之后到底怎么走”这个岔路口较劲。Abaqus里的cohesive方法之所以长期被主流研究者和工程团队拿来当作主力工具是因为它能在同一套框架里把裂缝的起裂、扩展、张开以及两条裂缝相遇后的穿过、俘获、转向这些行为全部串起来不需要像连续体离散元那样从头折腾接触算法也不用像XFEM那样在交叉处纠结富集函数的配置。这篇文章我就围绕“Abaqus模拟中水力裂缝与天然裂缝相交的cohesive行为”这件事把自己的建模思路、参数设置逻辑、求解器配置以及踩过的坑完整摊开来讲适合正在做页岩储层压裂模拟、研究裂缝网络扩展规律或者被Abaqus收敛性问题折磨得睡不好觉的朋友参考。1. 为什么这个看似“只是加了一条缝”的问题这么难搞1.1 现场真正关心的是裂缝相遇后走哪条路压裂施工中真正影响改造体积的往往不是你造出来的那条主缝本身而是它遇到储层里既有天然裂缝、层理面、节理之后发生的行为。如果水力裂缝直接穿过天然裂缝那么天然裂缝基本不参与导流改造区域就是单一裂缝周围的窄条。如果水力裂缝被天然裂缝截住、打开并沿着弱面滑动那么流体就会沿着天然裂缝方向灌进去把原本不太起眼的旧缝转变成新的导流通道这才会形成复杂的缝网。所以做数值模拟的时候判断标准的重点不是某一条裂缝扩展得多快而是两条裂缝相交那一瞬间应力场和流体压力如何共同决定裂缝走向。在Abaqus里cohesive方法的两条路线——cohesive元素和cohesive surface——都能表达裂缝界面的粘聚强度和断裂能但天然裂缝和水力裂缝在物理上有本质区别。水力裂缝是“从无到有”的它需要起裂条件扩展过程中尖端有流体压力驱动裂缝面还会存在滤失。天然裂缝更像一个“一直在那里”的弱面它不会主动扩展但在水力裂缝逼近时附近应力扰动会导致它提前发生剪切滑移或者拉伸张开。两者在cohesive框架里的角色完全不一样初始化方式、损伤判断和接触响应都得区别对待。很多初学者第一次搭这个模型老老实实建了两条裂缝路径结果发现天然裂缝根本不“听话”水力裂缝到了交叉点直接无视天然裂缝冲过去了或者天然裂缝在远场应力还没怎么变化时就全裂开。问题基本出在“把天然裂缝当成了一条和水力裂缝一样的预置裂缝”而不是把它当作一个应力状态敏感的弱面。1.2 cohesive方法为什么在这个场景里比XFEM更顺手如果只模拟一条水力裂缝的扩展XFEM和cohesive元素都能干互相竞争的点主要是网格依赖性和收敛成本。一旦问题变成“两条裂缝相交”XFEM其实会非常棘手。两条裂缝在几何上相交意味着富集域重叠裂缝交叉处的水平集函数要额外处理Abaqus里对完全相交的富集裂缝支持得并不完善尤其是涉及裂缝面接触和流体连续性问题时收敛极其痛苦。cohesive方法的优势在于它把裂缝路径在网格层面或者接触层面显式表达出来天然裂缝和水力裂缝的几何关系是预先定义好的——要么共用节点要么通过接触面耦合。这样相交处的力学关系就变成了一条简单的本构判定水力裂缝尖端到达天然裂缝表面时外侧主应力、孔隙压力、天然裂缝的残余强度共同决定开裂方向。你不需要操心复杂的富集函数只需要保证两条缝在网格划分时共享交界面然后正确地赋予天然裂缝一个比基质低的强度参数就行。另外从工程实用角度看cohesive方法允许你给天然裂缝的cohesive属性设置一个比水裂缝更低、或者完全不同的初始刚度和断裂能这是XFEM里不太好表达的东西。天然裂缝往往不是单纯的张开型裂缝而是带有较大摩擦角的剪切弱面。cohesive surface天然支持法向-切向耦合破坏准则能区分拉伸张开和剪切滑移这一点对相交行为模拟至关重要。2. 动手建模之前必须先想清楚的物理参数2.1 几何构建天然裂缝用cohesive element还是cohesive surface先说结论对于这个题目我更推荐天然裂缝用cohesive surface基于面的接触行为水力裂缝用cohesive element或者两条缝都用cohesive element但赋予不同材料参数。实际工作中需要根据你对碰撞区域的关注程度来选择。如果你关注的是水力裂缝在天然裂缝表面上的应力重分布以及天然裂缝局部张开和剪切滑移的细节那么cohesive element会给你更丰富的场变量输出SDEG、CSDMG、节点反力、开口量。代价是划分网格时必须把天然裂缝做成一条很薄的实体层通常用零厚度cohesive element来实现网格生成比较复杂尤其当天然裂缝是一条斜线穿越整个模型的时候。用cohesive surface的做法是在天然裂缝两侧定义接触面再在接触属性里绑定一个cohesive behavior的traction-separation本构。好处是省去了插入零厚度单元的麻烦和Abaqus的通用接触框架兼容得更好天然裂缝面的摩擦滑移也会自动承接。问题是cohesive surface里的损伤变量不像element那样可视化得那么直观后处理时需要额外提取相对位移和面法向应力来判断开没开。我在实际工作中偏向于——如果天然裂缝数量少、交汇点唯一用cohesive element如果天然裂缝数量多或者缝网复杂比如几十条随机裂缝用cohesive surface。相交细节的重要性决定了你愿意付出多少网格复杂度这个权衡要靠你自己把控。2.2 材料参数牵引-分离本构里的三组硬骨头cohesive本构在Abaqus里说白了就三件事初始刚度、损伤起始强度、损伤演化断裂能。这三件事在水力裂缝与天然裂缝场景里都有各自的坑。第一初始刚度。对于零厚度cohesive element初始刚度需要取得“足够大但不至于引起数值病态”。一般用基质的弹性模量除以一个特征单元尺寸比如E/H但这个E/H估值经常会让刚度矩阵条件数爆炸。很多收敛问题追来追去最后发现是cohesive初始刚度比相邻实体单元刚度高了四五个数量级。我的经验是取相邻实体单元刚度的10到50倍即可不要盲目追求大刚度。第二损伤起始强度。这里有两个方向都要给法向抗拉强度和切向抗剪强度。对天然裂缝而言切向抗剪强度需要配合摩擦角来给不是随便给一个比基质低的数值就行。你可以用摩尔-库仑准则估算剪切强度粘聚力正应力×摩擦系数。如果天然裂缝面粗糙、有充填物摩擦角可能不小如果是平滑的页岩层理面摩擦角比较低。这个参数直接决定水力裂缝到达时天然裂缝是更容易发生拉伸张开还是剪切滑移。第三断裂能。双线性cohesive本构的核心假设是裂缝从完全张开到失去承载能力的软化段所消耗的能量等于裂缝扩展的临界能量释放率。对水力裂缝这个值可以用岩石I型断裂韧度K_IC换算G_IC K_IC²/E其中E是平面应变模量E/(1-v²)。对天然裂缝断裂能要低得多一般取基质断裂能的20%~50%有时甚至更低因为天然裂缝本来就是岩石中的弱面残余胶结强度很低。2.3 流体与力学的耦合设置在相交问题里至关重要水力裂缝的本质是流体压力驱动裂缝两侧张开Abaqus中要把孔压单元和cohesive单元的孔压自由度结合起来。注意cohesive element必须采用带孔压的单元类型比如COH2D4P平面应变二维这样流体才能在裂缝内流动并驱动缝面张开。天然裂缝是否需要考虑内部流体流动取决于它是否被水力裂缝打开。如果没有预置裂缝内流体天然裂缝一开始是“干”的只有在水力裂缝到达交叉点和它连通后压裂液才会灌进去。Abaqus里实现这个行为需要保证天然裂缝cohesive单元的孔压自由度和基质孔压自由度耦合否则压裂液压力没法传递到天然裂缝内部模拟出的穿行行为会失真。对于相交问题还有一个关键概念水力裂缝尖端到达天然裂缝面前时其实流体的前沿比应力扰动前沿慢——力学扰动先到流体后到。很多模拟结果里水力裂缝提前在天然裂缝另一个位置“穿透”原因就是应力在天然裂缝表面引发局部拉伸而流体还没有到达该处导致力学损伤先行、流体压力后补。所以要不要让cohesive损伤启动跟流体压力同步取决于你是模拟真实的压裂过程压力和流量由注入条件控制还是只想研究静态应力影响。我一般建议真实注入条件下的耦合求解这更贴合现场。3. 实操过程一步一步搭出HF-NF相交模型3.1 从最简二维模型开始不要一上来就建三维我强烈建议先把问题压成二维平面应变模型。一方面水力裂缝在垂直于最大主应力方向扩展二维模型已经能抓住主要机理另一方面二维天然裂缝与水力裂缝的相交关系在几何上就是两条线交叉网格和接触的设置难度低一个量级。模型尺寸建议取50 m × 50 m的方形区域最外侧边界距离水力裂缝端部至少10 m以上减少边界反射对应力场的干扰。中间一条预置水力裂缝路径从注入点一直画到模型边界附近另外一条斜向天然裂缝穿过它。可以设置天然裂缝与水力裂缝成30度、60度、90度三种典型交角通过对比来研究角度影响。网格分区时要把两条裂缝路径放在内部面internal partition上。如果不做partition直接画网格cohesive element没法只沿着裂缝路径产生。Abaqus里可以先用Partition Edge at Vertex或者指定内部面把模型分割成若干区域再在裂缝路径区域单独赋予cohesive元素类型。对于水力裂缝和天然裂缝交汇点网格要加密。交汇点周围至少布置5个以上单元尺寸一致的cohesive element避免单元尺寸跳变引起的虚假应力集中。我试过在交汇点附近粗网格下模拟天然裂缝的开裂模式会被网格尺寸硬生生“逼”出来换一套粗网格结果就变了这种网格无关性问题在相交模拟里尤其严重。3.2 单元类型、材料赋值的坑基体采用CPE4P平面应变带孔压单元天然裂缝和水力裂缝路径上的cohesive单元采用COH2D4P。注意cohesive单元必须有一个方向指向裂缝切口法向否则应力会被错误传递。检查方法是在Abaqus/CAE里显示单元法向确保cohesive单元厚度方向初始牵引方向与裂缝面法向一致。材料赋值上基体给一个简单的线弹性模型杨氏模量E比如30 GPa、泊松比v0.25、渗透率k比如1e-17 m²。cohesive材料要单独定义即便水力裂缝和天然裂缝都是用cohesive element也一定要给两种不同的cohesive材料因为它们的损伤起始强度、断裂能、渗透率完全不同。具体数值我以一个典型页岩场景为例基质抗拉强度5 MPaI型断裂韧度1.2 MPa·m^0.5估算出基质的断裂能约40 N/m。水力裂缝cohesive参数取抗拉强度5 MPa、断裂能40 N/m天然裂缝取抗拉强度1.5 MPa、剪切强度3 MPa、断裂能12 N/m。这样的参数组合下天然裂缝在交叉点附近更易于被激活又不会在远场应力下提前破坏。3.3 地应力平衡这一步做不好后面全是白搭水力压裂模拟最忌讳的是跳过地应力平衡直接加注水压力。如果你直接给模型施加远场应力模型会先发生瞬态变形这个变形产生的位移会在后续孔压计算中造成虚假的裂缝开度得到完全错误的缝宽。正确的做法是在第一个分析步Geostatic里施加远场应力比如水平最小主应力15 MPa、水平最大主应力20 MPa、上覆压力25 MPa关闭孔压变化让模型达到初始平衡。静力平衡后孔隙压力的初始条件也要设置好——初始孔隙压力设为储层压力比如25 MPa取决于深度然后在地应力平衡分析步中用Gravity或者body force让总应力与有效应力达到协调。地应力平衡收敛不好时不要立刻去调网格。先检查初始应力状态是否满足平衡方程即由体力和边界条件产生的应力分布是否等于初始应力场。Abaqus里可以用UNSATURATED或者RESIDUAL选项检查不平衡力。能把地应力平衡做到不平衡力小于0.1%的精度后面的耦合收敛会顺畅很多。3.4 注入与求解设置流量控制比压力控制稳得多注入条件我推荐用流量控制在注入点位置给cohesive元素施加一个体积注入率流比如2e-4 m²/s二维单位面积流而不是直接固定孔压。固定孔压的坏处在于当裂缝展开后孔压在裂缝内高度非均匀你固定一个值就会人为地给裂缝尖端施加一个不切实际的压力分布。流量控制更接近现场施工泵注排量是实际可控的量。Abaqus中实现方式是在注入节点集合上施加集中流量concentrated flux或者给cohesive单元施加分布流。如果是二维cohesive element建议在注入点单元上直接加FLUX边界条件。注入速率的选择要兼顾计算效率和稳定性太快会导致压力急剧上升、convergence失败太慢计算时间长。可以先跑一个10秒注水试算看压力响应是否平滑。分析步类型选SoilsTransient瞬态孔压-应力耦合。时间步长要控制好初始增量步设1e-4秒量级最大增量步不要超过0.5~1秒。显式积分在cohesive严重损伤阶段经常出现负特征值隐式虽然慢但更稳定。收敛问题的终极手段是双管齐下一是打开自动稳定化Automatic stabilization设置一个很小的阻尼系数比如1e-5~1e-6二是增加cohesive本构的粘性正则化在损伤演化里加一个viscous regularization参数取1e-4~1e-5。这不会对结果产生实质性影响但能大幅改善裂纹扩展时的收敛性。3.5 后处理看什么SDEG不够还要看CSDMG和开度结果文件中重点监控三个量SDEG标量损伤刚度退化系数0到1、CSDMGcohesive surface的损伤变量、COPEN接触开度。对水力裂缝SDEG达到1说明该处cohesive单元已经失去承载能力裂缝完全形成对天然裂缝除了SDEG之外还要看它的相对滑移量CSLIP因为天然裂缝在压裂过程中大量以剪切滑移方式激活。判断水力裂缝有没有穿过天然裂缝主要还是看裂缝扩展路径——在交汇点之后水力裂缝是否沿着原来的裂缝路径继续扩展还是拐进天然裂缝方向。看孔隙压力演化云图也能辅助判断如果天然裂缝被激活并导流那么沿着天然裂缝方向会出现一条明显的高孔隙压力带而如果水力裂缝直接穿过天然裂缝另一侧的压力响应滞后明显。这个现象和现场微地震事件的空间分布可以对应起来。4. 相交行为的主要结果模式与判定方法4.1 穿过、俘获、转向——三种结果怎么区分水力裂缝与天然裂缝相交的经典理论框架是Blanton准则和Renshaw-Pollard准则。Blanton准则认为水力裂缝到达天然裂缝面时断裂能否保持在水力裂缝路径上取决于天然裂缝面的摩擦系数、逼近角度以及应力差。近似穿过的条件是水力裂缝尖端的应力强度因子仍然超过基质断裂韧度同时天然裂缝表面剪切应力没有达到摩擦滑动极限。如果天然裂缝表面剪切应力先达到极限裂缝就会被天然裂缝“俘获”发生滑移或张开。在Abaqus模拟里这三种结果都有非常清晰的特征穿过水力裂缝在交汇点不发生偏转继续沿原方向扩展天然裂缝只是出现轻微的剪切滑移但并未大规模张开压力曲线在交汇时刻没有明显拐点。俘获水力裂缝尖端到达天然裂缝后尖端钝化裂缝停止向前扩展流体压力持续升高直到天然裂缝被强行打开。压力曲线上交汇点处有明显压力增升天然裂缝张开度持续增长。转向水力裂缝沿天然裂缝表面滑移/shear扩展一段距离然后在天然裂缝尖端附近重新找到一个应力有利方位产生一条新的拉伸裂缝继续向远处扩展。这种模式在非正交交角比如60度下最为常见。4.2 参数敏感性最影响结果的是交角然后是应力差我用同一套模型跑过大大小小几十个算例最强烈的感受是——交角的影响几乎是一票否决制。当水力裂缝与天然裂缝交角接近90度时天然裂缝表面法向应力很高剪切应力接近零水力裂缝很容易穿过当交角小于45度时天然裂缝面的剪切应力分量显著增大俘获和转向的概率大增。这和室内真三轴实验观察到的结果高度一致。第二个敏感参数是水平主应力差。应力差大比如10 MPa以上时水力裂缝的方向约束力强即使被天然裂缝短暂俘获也会在绕行一段后拉回到最大主应力方向。应力差小比如2~3 MPa时裂缝就容易“横着走”导致天然裂缝全面激活。所以做现场参数研究时这两个参数一定要重点扫一遍。天然裂缝的cohesive强度参数其实排第三位。只有当天然裂缝的抗剪强度极低、或者断裂能极低时它才会从被动的“俘获点”变成主动的“扩展通道”。如果天然裂缝断裂能给得太高它就和基质没有区别模拟结果自然全部变成“穿过”。这里有一个容易被忽略的陷阱有些人为了模拟天然裂缝“承载能力弱”把初始刚度给得很低结果导致裂缝还未建立接触应力就已经发生数值穿透看起来像天然裂缝“提前打开了”其实是计算假象。4.3 与室内实验和微地震观测对照时应该对什么量数值模拟的价值最终要落到与实测数据的对照上。室内真三轴水力压裂实验能观察到裂缝穿过与否的宏观形态数值模拟要对照的变量是裂缝扩展路径角度和注入压力-时间曲线形状。如果你的模拟结果压力曲线在天然裂缝相交时刻出现了一个“平台段”或者“压力突降”而实验曲线也有类似特征说明流体进入了天然裂缝系统这个过程的重现就是模型可靠性的直接证明。现场微地震监测通常没法直接看到单条天然裂缝的开闭但事件点云的空间展布能间接反映缝网复杂度。在模拟中可以把cohesive单元的剪切损伤区域作为“微地震事件候选区”来统计对比事件云的空间扩展方向与模拟的天然裂缝激活范围。这一层面上的对照更偏向趋势验证不必追求逐点吻合。5. 常见问题与调试实录5.1 天然裂缝在远场应力下就提前全裂开了这个问题十有八九是断裂能给低了。天然裂缝虽然是弱面但在储层深处承受着几十兆帕的压实应力不会一碰就碎。如果模拟开始时天然裂缝就大片损伤首先要检查初始应力平衡后天然裂缝面上的正应力和剪切应力是多少然后反推你设置的cohesive参数是否在屈服面以内。另一个原因是初始刚度过低导致数值穿透。我建议在天然裂缝cohesive参数中先做一个弹性试算不加注入只加载地应力观察COPEN和CSLIP是否为零。如果非零说明要么cohesive初始刚度不够要么天然裂缝单元已被错误定义成了“空隙”。5.2 水力裂缝尖端到达天然裂缝后收敛性开始“狂飙”这是相交模拟最经典的一幕。水力裂缝尖端与天然裂缝面接触的瞬间几何拓扑关系发生突变Abaqus的Newton迭代会在几个增量步内反复震荡。此时可以做三件事把最大增量步缩小一个量级给cohesive损伤演化增加粘性正则化系数从1e-5开始试探在天然裂缝的接触属性中打开自动过盈调整adjust only to remove overclosure避免尖端点位落入天然裂缝面网格间隙时产生虚假接触力。如果还是收敛不动就把水力裂缝尖端附近网格进一步加密让cohesive单元的尺寸降到之前的一半。这个位置的收敛问题绝大多数是应力奇异性引起的网格越粗奇异应力越大收敛越难。5.3 注入压力曲线出现锯齿状振荡压力曲线的锯齿一般来自cohesive单元逐个失效时释放的弹性应变能这在隐式求解里表现为节点失衡力的突然升高。如果锯齿振幅不大可以接受如果振幅大到压力回退要考虑减小增量步长或者检查损伤演化断裂能是否低得过于极端。还有另一种锯齿源于网格效应cohesive单元的失效是逐单元触发的失败单元越多锯齿频率越高。这是cohesive方法的固有特性不可能完全消除。工程判断上应该关注压力曲线的包络趋势而不是单个锯齿点。5.4 流体沿着cohesive表面非物理泄漏这个问题很多人会忽视。cohesive element在完全损伤后如果启用了contact裂缝面之间的流体渗流需要单独定义否则会出现流体跨裂缝面“直接穿过去”的非物理现象。Abaqus中要在cohesive element的损伤后行为里设置好渗透系数和失效力学的相互作用或者将cohesive surface与flow properties关联起来。另一种泄漏来自cohesive单元向实体单元的孔压传递。如果cohesive element与相邻孔压单元的孔压自由度没有连接好流体压力就会出现断崖式跌落。排查方法是看裂缝壁面的孔压云图是否有明显的不连续突变如果有多半是孔压自由度不共享导致的。5.5 零厚度cohesive element大面积翻转或负体积零厚度cohesive element在严重剪切变形下容易翻转成负体积这个在相交行为中尤其容易发生因为天然裂缝激活时剪切滑移量非常大。解决方法有三个方向一是采用“几何厚度”方法给cohesive element一个极小但非零的初始厚度比如1e-5 m缓解翻转二是检查单元法向设置确保cohesive单元的厚度方向沿着裂缝面法向三是控制损伤后的最大张开度在cohesive behavior中设置合适的ccdamage后接触刚度。我个人最推荐的还是第三条路线因为翻转本质上是单元在压缩状态下失去抵抗能力导致的设置损伤后的接触约束比单纯给初始厚度更物理。6. 参数标定与验证的整体方法论6.1 先复现单一水力裂缝再叠加天然裂缝我踩过最大的坑就是一上来就搭“基质水力裂缝天然裂缝流体耦合”全套模型然后被无数个同时爆炸的问题淹没。正确路径是分四步走第一步只模拟单条水力裂缝在均质储层中的扩展验证注入压力曲线和裂缝宽度是否符合经典的KGD或PKN解析解第二步加入天然裂缝但不注水只加载地应力观察天然裂缝在应力场中的稳定性第三步水力裂缝和天然裂缝合并但不加流体耦合只给定缝内压力分布研究力学响应第四步打开完整流固耦合。每一步都验证通过了再进入下一步。这么做的好处是你始终知道当前的问题出在哪一层而不是把所有不确定性混在一起反复试错。6.2 Kalman式的参数校准思路cohesive参数虽然有物理意义但在数值模型里多少存在“等效参数”的意味。因此我习惯在模型标定时采用“宏观响应匹配”的思路先用室内实验的注入压力曲线作为目标调cohesive断裂能和强度把压力响应拟合到合理范围再以实验后的裂缝形态为目标校核天然裂缝的激活阈值。这两个目标能同时满足时参数组合的可信度就很高。唯一要提醒的是cohesive参数不要单独某个变量反复微调因为它的响应是高度耦合的。比如降低断裂能可以让压力曲线更早出现拐点但同时也会让裂缝宽度变窄。这时候要同步观察压力曲线和缝宽两个指标才能判断参数调整是否真的指向正确的方向。7. 基于个人经验的收尾建议7.1 预留试错时间开头先用最简单的参数组合跑通全流程如果只能给一条建议那就是第一次跑通全流程之前千万别把参数调到“科学论文级别”的复杂。先用一个非常粗的网格、非常保守的注入速率、非常常规的cohesive参数把模型从地应力平衡到裂缝扩展完整跑通哪怕结果粗糙得和你预期的完全不像只要流程通了后面所有优化都只是替换参数而已。我第一次做相交模型时在参数调优上花了三周结果大量时间浪费在“参数没出现问题、而是材料赋值漏了一项导致行为完全错误”这种低级问题上。7.2 做相交模拟输出频率一定要高相交过程的持续时间往往非常短冲击效应集中在交汇后的零点几秒到几秒内。后处理里的时间步增量如果按照整个注水过程的尺度来设置很可能完全错过交汇时刻的细节。我建议把输出频率设置为每增量步输出一次至少在交汇点附近的时间段内加密输出。这是花钱最少但回报最大的一个设置能让你看到压力降、天然裂缝张开过程和水力裂缝偏转的完整动画。7.3 扩展方向从单一天然裂缝走向缝网这个模型跑通后最自然的扩展方向是把单一天然裂缝变成一组随机裂缝网络。你可以用Python脚本在Abaqus里批量生成随机方向、随机长度、随机间距的天然裂缝然后逐个赋予cohesive surface属性。此时计算成本会明显上升但你能观察到一个比单缝相交丰富得多的现象水力裂缝在缝网中反复分流、转向、被俘获最终形成体积更大的改造区。这才是从“单缝机理研究”走向“缝网形成机理研究”的关键一步。做这类模型两年多我的一个深切体会是Abaqus本身不难难的是你想清楚每一条缝在物理上到底承载什么角色然后把它映射成正确的cohesive参数和接触设置。水力裂缝是主动进攻的一方天然裂缝是被动响应的一方两者在交汇点上博弈而这个博弈的胜负往往在你设定材料参数的那一刻就已经注定了。所以花在参数标定上的时间永远值得别急着跑大规模算例。 Abaqus模拟水力裂缝与天然裂缝相交的cohesive行为要点
搞水力压裂数值模拟的人很多时候不是在跟最终输出曲线较劲而是在跟“水力裂缝撞上天然裂缝之后到底怎么走”这个岔路口较劲。Abaqus里的cohesive方法之所以长期被主流研究者和工程团队拿来当作主力工具是因为它能在同一套框架里把裂缝的起裂、扩展、张开以及两条裂缝相遇后的穿过、俘获、转向这些行为全部串起来不需要像连续体离散元那样从头折腾接触算法也不用像XFEM那样在交叉处纠结富集函数的配置。这篇文章我就围绕“Abaqus模拟中水力裂缝与天然裂缝相交的cohesive行为”这件事把自己的建模思路、参数设置逻辑、求解器配置以及踩过的坑完整摊开来讲适合正在做页岩储层压裂模拟、研究裂缝网络扩展规律或者被Abaqus收敛性问题折磨得睡不好觉的朋友参考。1. 为什么这个看似“只是加了一条缝”的问题这么难搞1.1 现场真正关心的是裂缝相遇后走哪条路压裂施工中真正影响改造体积的往往不是你造出来的那条主缝本身而是它遇到储层里既有天然裂缝、层理面、节理之后发生的行为。如果水力裂缝直接穿过天然裂缝那么天然裂缝基本不参与导流改造区域就是单一裂缝周围的窄条。如果水力裂缝被天然裂缝截住、打开并沿着弱面滑动那么流体就会沿着天然裂缝方向灌进去把原本不太起眼的旧缝转变成新的导流通道这才会形成复杂的缝网。所以做数值模拟的时候判断标准的重点不是某一条裂缝扩展得多快而是两条裂缝相交那一瞬间应力场和流体压力如何共同决定裂缝走向。在Abaqus里cohesive方法的两条路线——cohesive元素和cohesive surface——都能表达裂缝界面的粘聚强度和断裂能但天然裂缝和水力裂缝在物理上有本质区别。水力裂缝是“从无到有”的它需要起裂条件扩展过程中尖端有流体压力驱动裂缝面还会存在滤失。天然裂缝更像一个“一直在那里”的弱面它不会主动扩展但在水力裂缝逼近时附近应力扰动会导致它提前发生剪切滑移或者拉伸张开。两者在cohesive框架里的角色完全不一样初始化方式、损伤判断和接触响应都得区别对待。很多初学者第一次搭这个模型老老实实建了两条裂缝路径结果发现天然裂缝根本不“听话”水力裂缝到了交叉点直接无视天然裂缝冲过去了或者天然裂缝在远场应力还没怎么变化时就全裂开。问题基本出在“把天然裂缝当成了一条和水力裂缝一样的预置裂缝”而不是把它当作一个应力状态敏感的弱面。1.2 cohesive方法为什么在这个场景里比XFEM更顺手如果只模拟一条水力裂缝的扩展XFEM和cohesive元素都能干互相竞争的点主要是网格依赖性和收敛成本。一旦问题变成“两条裂缝相交”XFEM其实会非常棘手。两条裂缝在几何上相交意味着富集域重叠裂缝交叉处的水平集函数要额外处理Abaqus里对完全相交的富集裂缝支持得并不完善尤其是涉及裂缝面接触和流体连续性问题时收敛极其痛苦。cohesive方法的优势在于它把裂缝路径在网格层面或者接触层面显式表达出来天然裂缝和水力裂缝的几何关系是预先定义好的——要么共用节点要么通过接触面耦合。这样相交处的力学关系就变成了一条简单的本构判定水力裂缝尖端到达天然裂缝表面时外侧主应力、孔隙压力、天然裂缝的残余强度共同决定开裂方向。你不需要操心复杂的富集函数只需要保证两条缝在网格划分时共享交界面然后正确地赋予天然裂缝一个比基质低的强度参数就行。另外从工程实用角度看cohesive方法允许你给天然裂缝的cohesive属性设置一个比水裂缝更低、或者完全不同的初始刚度和断裂能这是XFEM里不太好表达的东西。天然裂缝往往不是单纯的张开型裂缝而是带有较大摩擦角的剪切弱面。cohesive surface天然支持法向-切向耦合破坏准则能区分拉伸张开和剪切滑移这一点对相交行为模拟至关重要。2. 动手建模之前必须先想清楚的物理参数2.1 几何构建天然裂缝用cohesive element还是cohesive surface先说结论对于这个题目我更推荐天然裂缝用cohesive surface基于面的接触行为水力裂缝用cohesive element或者两条缝都用cohesive element但赋予不同材料参数。实际工作中需要根据你对碰撞区域的关注程度来选择。如果你关注的是水力裂缝在天然裂缝表面上的应力重分布以及天然裂缝局部张开和剪切滑移的细节那么cohesive element会给你更丰富的场变量输出SDEG、CSDMG、节点反力、开口量。代价是划分网格时必须把天然裂缝做成一条很薄的实体层通常用零厚度cohesive element来实现网格生成比较复杂尤其当天然裂缝是一条斜线穿越整个模型的时候。用cohesive surface的做法是在天然裂缝两侧定义接触面再在接触属性里绑定一个cohesive behavior的traction-separation本构。好处是省去了插入零厚度单元的麻烦和Abaqus的通用接触框架兼容得更好天然裂缝面的摩擦滑移也会自动承接。问题是cohesive surface里的损伤变量不像element那样可视化得那么直观后处理时需要额外提取相对位移和面法向应力来判断开没开。我在实际工作中偏向于——如果天然裂缝数量少、交汇点唯一用cohesive element如果天然裂缝数量多或者缝网复杂比如几十条随机裂缝用cohesive surface。相交细节的重要性决定了你愿意付出多少网格复杂度这个权衡要靠你自己把控。2.2 材料参数牵引-分离本构里的三组硬骨头cohesive本构在Abaqus里说白了就三件事初始刚度、损伤起始强度、损伤演化断裂能。这三件事在水力裂缝与天然裂缝场景里都有各自的坑。第一初始刚度。对于零厚度cohesive element初始刚度需要取得“足够大但不至于引起数值病态”。一般用基质的弹性模量除以一个特征单元尺寸比如E/H但这个E/H估值经常会让刚度矩阵条件数爆炸。很多收敛问题追来追去最后发现是cohesive初始刚度比相邻实体单元刚度高了四五个数量级。我的经验是取相邻实体单元刚度的10到50倍即可不要盲目追求大刚度。第二损伤起始强度。这里有两个方向都要给法向抗拉强度和切向抗剪强度。对天然裂缝而言切向抗剪强度需要配合摩擦角来给不是随便给一个比基质低的数值就行。你可以用摩尔-库仑准则估算剪切强度粘聚力正应力×摩擦系数。如果天然裂缝面粗糙、有充填物摩擦角可能不小如果是平滑的页岩层理面摩擦角比较低。这个参数直接决定水力裂缝到达时天然裂缝是更容易发生拉伸张开还是剪切滑移。第三断裂能。双线性cohesive本构的核心假设是裂缝从完全张开到失去承载能力的软化段所消耗的能量等于裂缝扩展的临界能量释放率。对水力裂缝这个值可以用岩石I型断裂韧度K_IC换算G_IC K_IC²/E其中E是平面应变模量E/(1-v²)。对天然裂缝断裂能要低得多一般取基质断裂能的20%~50%有时甚至更低因为天然裂缝本来就是岩石中的弱面残余胶结强度很低。2.3 流体与力学的耦合设置在相交问题里至关重要水力裂缝的本质是流体压力驱动裂缝两侧张开Abaqus中要把孔压单元和cohesive单元的孔压自由度结合起来。注意cohesive element必须采用带孔压的单元类型比如COH2D4P平面应变二维这样流体才能在裂缝内流动并驱动缝面张开。天然裂缝是否需要考虑内部流体流动取决于它是否被水力裂缝打开。如果没有预置裂缝内流体天然裂缝一开始是“干”的只有在水力裂缝到达交叉点和它连通后压裂液才会灌进去。Abaqus里实现这个行为需要保证天然裂缝cohesive单元的孔压自由度和基质孔压自由度耦合否则压裂液压力没法传递到天然裂缝内部模拟出的穿行行为会失真。对于相交问题还有一个关键概念水力裂缝尖端到达天然裂缝面前时其实流体的前沿比应力扰动前沿慢——力学扰动先到流体后到。很多模拟结果里水力裂缝提前在天然裂缝另一个位置“穿透”原因就是应力在天然裂缝表面引发局部拉伸而流体还没有到达该处导致力学损伤先行、流体压力后补。所以要不要让cohesive损伤启动跟流体压力同步取决于你是模拟真实的压裂过程压力和流量由注入条件控制还是只想研究静态应力影响。我一般建议真实注入条件下的耦合求解这更贴合现场。3. 实操过程一步一步搭出HF-NF相交模型3.1 从最简二维模型开始不要一上来就建三维我强烈建议先把问题压成二维平面应变模型。一方面水力裂缝在垂直于最大主应力方向扩展二维模型已经能抓住主要机理另一方面二维天然裂缝与水力裂缝的相交关系在几何上就是两条线交叉网格和接触的设置难度低一个量级。模型尺寸建议取50 m × 50 m的方形区域最外侧边界距离水力裂缝端部至少10 m以上减少边界反射对应力场的干扰。中间一条预置水力裂缝路径从注入点一直画到模型边界附近另外一条斜向天然裂缝穿过它。可以设置天然裂缝与水力裂缝成30度、60度、90度三种典型交角通过对比来研究角度影响。网格分区时要把两条裂缝路径放在内部面internal partition上。如果不做partition直接画网格cohesive element没法只沿着裂缝路径产生。Abaqus里可以先用Partition Edge at Vertex或者指定内部面把模型分割成若干区域再在裂缝路径区域单独赋予cohesive元素类型。对于水力裂缝和天然裂缝交汇点网格要加密。交汇点周围至少布置5个以上单元尺寸一致的cohesive element避免单元尺寸跳变引起的虚假应力集中。我试过在交汇点附近粗网格下模拟天然裂缝的开裂模式会被网格尺寸硬生生“逼”出来换一套粗网格结果就变了这种网格无关性问题在相交模拟里尤其严重。3.2 单元类型、材料赋值的坑基体采用CPE4P平面应变带孔压单元天然裂缝和水力裂缝路径上的cohesive单元采用COH2D4P。注意cohesive单元必须有一个方向指向裂缝切口法向否则应力会被错误传递。检查方法是在Abaqus/CAE里显示单元法向确保cohesive单元厚度方向初始牵引方向与裂缝面法向一致。材料赋值上基体给一个简单的线弹性模型杨氏模量E比如30 GPa、泊松比v0.25、渗透率k比如1e-17 m²。cohesive材料要单独定义即便水力裂缝和天然裂缝都是用cohesive element也一定要给两种不同的cohesive材料因为它们的损伤起始强度、断裂能、渗透率完全不同。具体数值我以一个典型页岩场景为例基质抗拉强度5 MPaI型断裂韧度1.2 MPa·m^0.5估算出基质的断裂能约40 N/m。水力裂缝cohesive参数取抗拉强度5 MPa、断裂能40 N/m天然裂缝取抗拉强度1.5 MPa、剪切强度3 MPa、断裂能12 N/m。这样的参数组合下天然裂缝在交叉点附近更易于被激活又不会在远场应力下提前破坏。3.3 地应力平衡这一步做不好后面全是白搭水力压裂模拟最忌讳的是跳过地应力平衡直接加注水压力。如果你直接给模型施加远场应力模型会先发生瞬态变形这个变形产生的位移会在后续孔压计算中造成虚假的裂缝开度得到完全错误的缝宽。正确的做法是在第一个分析步Geostatic里施加远场应力比如水平最小主应力15 MPa、水平最大主应力20 MPa、上覆压力25 MPa关闭孔压变化让模型达到初始平衡。静力平衡后孔隙压力的初始条件也要设置好——初始孔隙压力设为储层压力比如25 MPa取决于深度然后在地应力平衡分析步中用Gravity或者body force让总应力与有效应力达到协调。地应力平衡收敛不好时不要立刻去调网格。先检查初始应力状态是否满足平衡方程即由体力和边界条件产生的应力分布是否等于初始应力场。Abaqus里可以用UNSATURATED或者RESIDUAL选项检查不平衡力。能把地应力平衡做到不平衡力小于0.1%的精度后面的耦合收敛会顺畅很多。3.4 注入与求解设置流量控制比压力控制稳得多注入条件我推荐用流量控制在注入点位置给cohesive元素施加一个体积注入率流比如2e-4 m²/s二维单位面积流而不是直接固定孔压。固定孔压的坏处在于当裂缝展开后孔压在裂缝内高度非均匀你固定一个值就会人为地给裂缝尖端施加一个不切实际的压力分布。流量控制更接近现场施工泵注排量是实际可控的量。Abaqus中实现方式是在注入节点集合上施加集中流量concentrated flux或者给cohesive单元施加分布流。如果是二维cohesive element建议在注入点单元上直接加FLUX边界条件。注入速率的选择要兼顾计算效率和稳定性太快会导致压力急剧上升、convergence失败太慢计算时间长。可以先跑一个10秒注水试算看压力响应是否平滑。分析步类型选SoilsTransient瞬态孔压-应力耦合。时间步长要控制好初始增量步设1e-4秒量级最大增量步不要超过0.5~1秒。显式积分在cohesive严重损伤阶段经常出现负特征值隐式虽然慢但更稳定。收敛问题的终极手段是双管齐下一是打开自动稳定化Automatic stabilization设置一个很小的阻尼系数比如1e-5~1e-6二是增加cohesive本构的粘性正则化在损伤演化里加一个viscous regularization参数取1e-4~1e-5。这不会对结果产生实质性影响但能大幅改善裂纹扩展时的收敛性。3.5 后处理看什么SDEG不够还要看CSDMG和开度结果文件中重点监控三个量SDEG标量损伤刚度退化系数0到1、CSDMGcohesive surface的损伤变量、COPEN接触开度。对水力裂缝SDEG达到1说明该处cohesive单元已经失去承载能力裂缝完全形成对天然裂缝除了SDEG之外还要看它的相对滑移量CSLIP因为天然裂缝在压裂过程中大量以剪切滑移方式激活。判断水力裂缝有没有穿过天然裂缝主要还是看裂缝扩展路径——在交汇点之后水力裂缝是否沿着原来的裂缝路径继续扩展还是拐进天然裂缝方向。看孔隙压力演化云图也能辅助判断如果天然裂缝被激活并导流那么沿着天然裂缝方向会出现一条明显的高孔隙压力带而如果水力裂缝直接穿过天然裂缝另一侧的压力响应滞后明显。这个现象和现场微地震事件的空间分布可以对应起来。4. 相交行为的主要结果模式与判定方法4.1 穿过、俘获、转向——三种结果怎么区分水力裂缝与天然裂缝相交的经典理论框架是Blanton准则和Renshaw-Pollard准则。Blanton准则认为水力裂缝到达天然裂缝面时断裂能否保持在水力裂缝路径上取决于天然裂缝面的摩擦系数、逼近角度以及应力差。近似穿过的条件是水力裂缝尖端的应力强度因子仍然超过基质断裂韧度同时天然裂缝表面剪切应力没有达到摩擦滑动极限。如果天然裂缝表面剪切应力先达到极限裂缝就会被天然裂缝“俘获”发生滑移或张开。在Abaqus模拟里这三种结果都有非常清晰的特征穿过水力裂缝在交汇点不发生偏转继续沿原方向扩展天然裂缝只是出现轻微的剪切滑移但并未大规模张开压力曲线在交汇时刻没有明显拐点。俘获水力裂缝尖端到达天然裂缝后尖端钝化裂缝停止向前扩展流体压力持续升高直到天然裂缝被强行打开。压力曲线上交汇点处有明显压力增升天然裂缝张开度持续增长。转向水力裂缝沿天然裂缝表面滑移/shear扩展一段距离然后在天然裂缝尖端附近重新找到一个应力有利方位产生一条新的拉伸裂缝继续向远处扩展。这种模式在非正交交角比如60度下最为常见。4.2 参数敏感性最影响结果的是交角然后是应力差我用同一套模型跑过大大小小几十个算例最强烈的感受是——交角的影响几乎是一票否决制。当水力裂缝与天然裂缝交角接近90度时天然裂缝表面法向应力很高剪切应力接近零水力裂缝很容易穿过当交角小于45度时天然裂缝面的剪切应力分量显著增大俘获和转向的概率大增。这和室内真三轴实验观察到的结果高度一致。第二个敏感参数是水平主应力差。应力差大比如10 MPa以上时水力裂缝的方向约束力强即使被天然裂缝短暂俘获也会在绕行一段后拉回到最大主应力方向。应力差小比如2~3 MPa时裂缝就容易“横着走”导致天然裂缝全面激活。所以做现场参数研究时这两个参数一定要重点扫一遍。天然裂缝的cohesive强度参数其实排第三位。只有当天然裂缝的抗剪强度极低、或者断裂能极低时它才会从被动的“俘获点”变成主动的“扩展通道”。如果天然裂缝断裂能给得太高它就和基质没有区别模拟结果自然全部变成“穿过”。这里有一个容易被忽略的陷阱有些人为了模拟天然裂缝“承载能力弱”把初始刚度给得很低结果导致裂缝还未建立接触应力就已经发生数值穿透看起来像天然裂缝“提前打开了”其实是计算假象。4.3 与室内实验和微地震观测对照时应该对什么量数值模拟的价值最终要落到与实测数据的对照上。室内真三轴水力压裂实验能观察到裂缝穿过与否的宏观形态数值模拟要对照的变量是裂缝扩展路径角度和注入压力-时间曲线形状。如果你的模拟结果压力曲线在天然裂缝相交时刻出现了一个“平台段”或者“压力突降”而实验曲线也有类似特征说明流体进入了天然裂缝系统这个过程的重现就是模型可靠性的直接证明。现场微地震监测通常没法直接看到单条天然裂缝的开闭但事件点云的空间展布能间接反映缝网复杂度。在模拟中可以把cohesive单元的剪切损伤区域作为“微地震事件候选区”来统计对比事件云的空间扩展方向与模拟的天然裂缝激活范围。这一层面上的对照更偏向趋势验证不必追求逐点吻合。5. 常见问题与调试实录5.1 天然裂缝在远场应力下就提前全裂开了这个问题十有八九是断裂能给低了。天然裂缝虽然是弱面但在储层深处承受着几十兆帕的压实应力不会一碰就碎。如果模拟开始时天然裂缝就大片损伤首先要检查初始应力平衡后天然裂缝面上的正应力和剪切应力是多少然后反推你设置的cohesive参数是否在屈服面以内。另一个原因是初始刚度过低导致数值穿透。我建议在天然裂缝cohesive参数中先做一个弹性试算不加注入只加载地应力观察COPEN和CSLIP是否为零。如果非零说明要么cohesive初始刚度不够要么天然裂缝单元已被错误定义成了“空隙”。5.2 水力裂缝尖端到达天然裂缝后收敛性开始“狂飙”这是相交模拟最经典的一幕。水力裂缝尖端与天然裂缝面接触的瞬间几何拓扑关系发生突变Abaqus的Newton迭代会在几个增量步内反复震荡。此时可以做三件事把最大增量步缩小一个量级给cohesive损伤演化增加粘性正则化系数从1e-5开始试探在天然裂缝的接触属性中打开自动过盈调整adjust only to remove overclosure避免尖端点位落入天然裂缝面网格间隙时产生虚假接触力。如果还是收敛不动就把水力裂缝尖端附近网格进一步加密让cohesive单元的尺寸降到之前的一半。这个位置的收敛问题绝大多数是应力奇异性引起的网格越粗奇异应力越大收敛越难。5.3 注入压力曲线出现锯齿状振荡压力曲线的锯齿一般来自cohesive单元逐个失效时释放的弹性应变能这在隐式求解里表现为节点失衡力的突然升高。如果锯齿振幅不大可以接受如果振幅大到压力回退要考虑减小增量步长或者检查损伤演化断裂能是否低得过于极端。还有另一种锯齿源于网格效应cohesive单元的失效是逐单元触发的失败单元越多锯齿频率越高。这是cohesive方法的固有特性不可能完全消除。工程判断上应该关注压力曲线的包络趋势而不是单个锯齿点。5.4 流体沿着cohesive表面非物理泄漏这个问题很多人会忽视。cohesive element在完全损伤后如果启用了contact裂缝面之间的流体渗流需要单独定义否则会出现流体跨裂缝面“直接穿过去”的非物理现象。Abaqus中要在cohesive element的损伤后行为里设置好渗透系数和失效力学的相互作用或者将cohesive surface与flow properties关联起来。另一种泄漏来自cohesive单元向实体单元的孔压传递。如果cohesive element与相邻孔压单元的孔压自由度没有连接好流体压力就会出现断崖式跌落。排查方法是看裂缝壁面的孔压云图是否有明显的不连续突变如果有多半是孔压自由度不共享导致的。5.5 零厚度cohesive element大面积翻转或负体积零厚度cohesive element在严重剪切变形下容易翻转成负体积这个在相交行为中尤其容易发生因为天然裂缝激活时剪切滑移量非常大。解决方法有三个方向一是采用“几何厚度”方法给cohesive element一个极小但非零的初始厚度比如1e-5 m缓解翻转二是检查单元法向设置确保cohesive单元的厚度方向沿着裂缝面法向三是控制损伤后的最大张开度在cohesive behavior中设置合适的ccdamage后接触刚度。我个人最推荐的还是第三条路线因为翻转本质上是单元在压缩状态下失去抵抗能力导致的设置损伤后的接触约束比单纯给初始厚度更物理。6. 参数标定与验证的整体方法论6.1 先复现单一水力裂缝再叠加天然裂缝我踩过最大的坑就是一上来就搭“基质水力裂缝天然裂缝流体耦合”全套模型然后被无数个同时爆炸的问题淹没。正确路径是分四步走第一步只模拟单条水力裂缝在均质储层中的扩展验证注入压力曲线和裂缝宽度是否符合经典的KGD或PKN解析解第二步加入天然裂缝但不注水只加载地应力观察天然裂缝在应力场中的稳定性第三步水力裂缝和天然裂缝合并但不加流体耦合只给定缝内压力分布研究力学响应第四步打开完整流固耦合。每一步都验证通过了再进入下一步。这么做的好处是你始终知道当前的问题出在哪一层而不是把所有不确定性混在一起反复试错。6.2 Kalman式的参数校准思路cohesive参数虽然有物理意义但在数值模型里多少存在“等效参数”的意味。因此我习惯在模型标定时采用“宏观响应匹配”的思路先用室内实验的注入压力曲线作为目标调cohesive断裂能和强度把压力响应拟合到合理范围再以实验后的裂缝形态为目标校核天然裂缝的激活阈值。这两个目标能同时满足时参数组合的可信度就很高。唯一要提醒的是cohesive参数不要单独某个变量反复微调因为它的响应是高度耦合的。比如降低断裂能可以让压力曲线更早出现拐点但同时也会让裂缝宽度变窄。这时候要同步观察压力曲线和缝宽两个指标才能判断参数调整是否真的指向正确的方向。7. 基于个人经验的收尾建议7.1 预留试错时间开头先用最简单的参数组合跑通全流程如果只能给一条建议那就是第一次跑通全流程之前千万别把参数调到“科学论文级别”的复杂。先用一个非常粗的网格、非常保守的注入速率、非常常规的cohesive参数把模型从地应力平衡到裂缝扩展完整跑通哪怕结果粗糙得和你预期的完全不像只要流程通了后面所有优化都只是替换参数而已。我第一次做相交模型时在参数调优上花了三周结果大量时间浪费在“参数没出现问题、而是材料赋值漏了一项导致行为完全错误”这种低级问题上。7.2 做相交模拟输出频率一定要高相交过程的持续时间往往非常短冲击效应集中在交汇后的零点几秒到几秒内。后处理里的时间步增量如果按照整个注水过程的尺度来设置很可能完全错过交汇时刻的细节。我建议把输出频率设置为每增量步输出一次至少在交汇点附近的时间段内加密输出。这是花钱最少但回报最大的一个设置能让你看到压力降、天然裂缝张开过程和水力裂缝偏转的完整动画。7.3 扩展方向从单一天然裂缝走向缝网这个模型跑通后最自然的扩展方向是把单一天然裂缝变成一组随机裂缝网络。你可以用Python脚本在Abaqus里批量生成随机方向、随机长度、随机间距的天然裂缝然后逐个赋予cohesive surface属性。此时计算成本会明显上升但你能观察到一个比单缝相交丰富得多的现象水力裂缝在缝网中反复分流、转向、被俘获最终形成体积更大的改造区。这才是从“单缝机理研究”走向“缝网形成机理研究”的关键一步。做这类模型两年多我的一个深切体会是Abaqus本身不难难的是你想清楚每一条缝在物理上到底承载什么角色然后把它映射成正确的cohesive参数和接触设置。水力裂缝是主动进攻的一方天然裂缝是被动响应的一方两者在交汇点上博弈而这个博弈的胜负往往在你设定材料参数的那一刻就已经注定了。所以花在参数标定上的时间永远值得别急着跑大规模算例。