abaqus水力压裂 1 基于Cohesive单元的二维水力压裂模拟 2基于Cohesive单...

abaqus水力压裂 1 基于Cohesive单元的二维水力压裂模拟 2基于Cohesive单元的三维水力压裂模拟 3水力裂缝与天然裂缝相交模拟-Cohesive单元法 4裂缝发育地层的水力压裂模拟-Cohesive单元法 5基于XFEM的水力裂缝转向模拟 6基于XFEM的水平井多段压裂裂缝的缝间干扰问题研究水力压裂数值模拟总让人头大今天咱们来点硬核实操。用Abaqus搞裂缝扩展核心玩法就两种Cohesive单元和XFEM。先说Cohesive这个万金油二维模型里塞上它就像给地层装了个智能裂缝开关。做二维水力压裂时材料定义脚本得这么玩mdb.models[Model-1].Material(nameCohesive) mdb.models[Model-1].materials[Cohesive].CohesiveZoneProperty(table((0.1, 10.0, 0.5),)) mdb.models[Model-1].materials[Cohesive].elastic.TypeTRACTION mdb.models[Model-1].materials[Cohesive].elastic.table((1e6, 1e6, 1e6),)这里弹性刚度设成1e6是个经验值太小容易数值发散太大又会影响裂缝张开速度。断裂能参数0.5对应页岩的典型值搞砂岩的兄弟得往上调到2.0左右。三维模型就刺激了网格质量直接决定会不会算崩。建议用扫掠网格配合粘性正则化mdb.models[Model-1].steps[Step-1].stabilizationScaleFactor1e-5 mdb.models[Model-1].steps[Step-1].stabilizationMethodDISSIPATED_ENERGY_FRACTION这个粘性系数要是设大了计算结果会失真设小了迭代十几次都收不敛试过凌晨三点调参调到怀疑人生的举个手。abaqus水力压裂 1 基于Cohesive单元的二维水力压裂模拟 2基于Cohesive单元的三维水力压裂模拟 3水力裂缝与天然裂缝相交模拟-Cohesive单元法 4裂缝发育地层的水力压裂模拟-Cohesive单元法 5基于XFEM的水力裂缝转向模拟 6基于XFEM的水平井多段压裂裂缝的缝间干扰问题研究遇到天然裂缝交叉的情况接触定义必须骚操作interaction mdb.models[Model-1].ContactProperty(Fracture-Contact) interaction.tangentialBehavior(formulationPENALTY, directionalityISOTROPIC) interaction.normalBehavior(pressureOverclosureHARD)这里罚函数刚度建议从材料刚度的1%开始试方向性参数用各向同性更稳当。实测用指数衰减接触反而容易在裂缝尖端出现鬼畜抖动。XFEM才是真正的裂缝魔术师搞转向模拟时注意这个参数mdb.models[Model-1].materials[Shale].xfemCrackGrowth True mdb.models[Model-1].materials[Shale].crackCurvatureType KINKING mdb.models[Model-1].materials[Shale].crackKinkingAngle 45.0设置45度转向角时别太耿直实际会根据应力场动态变化。遇到过转向角锁死的兄弟可以在材料属性里把软化定律从线性改成指数试试。水平井多缝干扰问题最能检验模型可靠性重点看这段流量分配代码mdb.models[Model-1].FluidExchangeProperty(nameMultiStage) mdb.models[Model-1].fluidExchangeProperties[MultiStage].dependenceON mdb.models[Model-1].fluidExchangeProperties[MultiStage].table((0.01, 0.0), (0.1, 0.5))流量-压力关系表里第二个参数0.5是缝间干扰系数实际压裂时超过0.3就会明显出现主缝抢流量现象。建议配合现场微地震数据反演这个参数比教科书上的推荐值靠谱得多。最后说个血泪教训千万别在Cohesive单元和XFEM混用的时候开几何非线性算到一半单元扭曲得亲妈都不认识。真要玩高阶的先把单元尺寸砍到1/3再说。