《Applied Sciences》：Inverse-Problem Approach for 3MA Electromagnetic NDT on Laser-Hardened Materials
编辑推荐：
非常规储层的强非均质性导致裂缝扩展复杂，给定量压裂评估带来挑战。本研究整合真三轴压裂实验、三维CT重建、多场耦合数值模拟和多元线性回归，研究了古龙页岩（Q1、Q9）和致密砂岩在超临界二氧化碳（SC-CO?）压裂下的裂缝行为。研究人员提出了一个包含分形维数、空间
非常规储层的强非均质性导致裂缝扩展复杂，给定量压裂评估带来挑战。本研究整合真三轴压裂实验、三维CT重建、多场耦合数值模拟和多元线性回归，研究了古龙页岩（Q1、Q9）和致密砂岩在超临界二氧化碳（SC-CO?）压裂下的裂缝行为。研究人员提出了一个包含分形维数、空间均匀性和开度分布的裂缝复杂性指数（FCI）作为定量指标。结果表明，与水相比，SC-CO?显著降低了破裂压力并增加了裂缝复杂性。对于Q9页岩，SC-CO?的破裂压力为32.91 MPa（比水低10.46%），分形维数为2.41，FCI为8.92。在致密砂岩中，SC-CO?破裂压力为34.12 MPa，而水将其升至44.50 MPa；分形维数和FCI仅为2.05和3.40，均低于水压裂的页岩。多元线性回归量化了贡献权重：岩性弱面发育主导裂缝复杂性（41.6%），远超脆性指数。注入速率主要控制压裂规模：当速率从10 mL/min增加到40 mL/min时，损伤面积比从1.79%升至2.90%。水平应力差是复杂性提升的关键：当应力差从0 MPa升至4 MPa时，分形维数从1.9230增至1.9901。数值模拟进一步揭示了热-流-固多场耦合效应。所提出的基于FCI的评估和回归模型为优化非均质非常规储层SC-CO?压裂设计提供了定量框架。
**论文解读文章**
**研究背景与问题**
随着全球能源需求演变，页岩和致密砂岩等低渗透储层的高效开发成为能源安全的关键挑战。水力压裂作为核心增产技术，暴露出水基流体的局限性：高耗水量、黏滞损伤和水锁效应。超临界二氧化碳（SC-CO?）压裂凭借低黏度、高扩散性、低破裂压力、复杂裂缝网络以及节水、减碳、降低储层损伤等额外优势，成为替代方案。然而，现有研究多集中于水基流体，对SC-CO?与不同岩性（尤其是高度非均质的陆相页岩如古龙页岩）相互作用的系统认识不足。不同岩性（页岩与致密砂岩）在SC-CO?下的裂缝行为及网络复杂性缺乏系统对比，多参数（流体相、注入参数、应力条件）的耦合影响与敏感性仍需探究。因此，本研究旨在通过物理实验与数值模拟结合，定量揭示岩性、流体和工程参数对裂缝扩展的控制机制，为SC-CO?压裂设计提供理论依据。
**研究方法与技术**
研究人员采用真三轴多流体压裂实验系统，对古龙页岩Q1层、Q9层和致密砂岩三种岩心进行压裂实验，涵盖SC-CO?、液态CO?和水三种流体，并设置不同应力状态和注入速率。实验后通过CT扫描进行三维裂缝重建，提取分形维数、裂缝密度、裂缝高度和开度变异系数，构建裂缝复杂性指数（FCI）。同时，基于Q9层数据，在COMSOL Multiphysics平台上建立考虑热-流-固-损伤多场耦合的数值模型，采用Weibull分布描述岩石非均质性，Span-Wagner状态方程和Vesovic输运模型实时更新CO?物性，模拟不同注入速率、流体温度和水平应力差下的裂缝扩展。最后，利用多元线性回归量化岩性、应力差、注入速率和温度对破裂压力和FCI的贡献权重。样本来源于古龙页岩Q1和Q9层及致密砂岩，实验数据包括矿物组成、力学参数和真三轴压裂结果。
**研究结果**
**2.1 岩性与力学参数**
通过XRD分析和三轴压缩实验，Q1页岩含40–50%黏土和30–40%石英，Q9页岩含50–60%黏土和20–30%石英，致密砂岩含45–55%石英和10–15%黏土。脆性指数梯度：Q9（30–40）< Q1（35–45）< 砂岩（55–65）。力学参数显示Q9页岩弹性模量最低（10.5–13.2 GPa）、泊松比最高（0.30–0.35），致密砂岩弹性模量最高（28.5–36.4 GPa）、泊松比最低（0.18–0.24）。尽管致密砂岩脆性指数高，但缺乏弱面，导致“高脆性、低复杂性”悖论，后续CT和数值分析证实弱面系统比脆性指数更决定性地控制裂缝网络复杂性。
**2.2 实验准备与流程**
多流体真三轴系统实现了注入压力、围压和温度的实时监测。试样加工为Φ50 mm×100 mm圆柱，钻孔并安装PEEK管，经泄漏测试、加载、加热、流体注入至破裂，再卸压冷却后进行CT扫描。实验设计涵盖三种岩性、三种流体、三种应力状态和不同注入速率，共完成19组有效测试。
**2.3.1 不同流体与岩性分析**
所有岩性下破裂压力顺序为：SC-CO? < 液态CO? < 水。Q9页岩对应值分别为32.91、36.65和36.83 MPa。致密砂岩流体敏感性更强：SC-CO?破裂压力34.12 MPa，液态CO?和水分别升至43.59和44.50 MPa（增幅28%和30%）。原因在于SC-CO?超低黏度和近零界面张力使其能渗透微纳米孔喉，形成均匀孔隙压力场；而水黏度高，压力集中井筒附近。CT三维重建显示，Q9页岩形成多分支曲折网络，Q1页岩以沿最大主应力方向主裂缝为主，致密砂岩形成单一平直裂缝。
**2.3.2 不同注入速率分析**
注入速率对破裂压力影响因岩性而异。Q1页岩和致密砂岩中破裂压力随速率单调增加；致密砂岩速率从2升至6 mL/min时破裂压力从33.27升至38.42 MPa，且裂缝形态简化。Q9页岩呈V型趋势，中速率（4 mL/min）时裂缝网络最复杂，表明弱面存在允许“扩散-积累”平衡。
**2.3.3 不同应力条件分析**
当轴向应力小于围压时，破裂压力升高；相等时最低。Q1页岩在非对称应力下分形维数达2.37（比对称应力高约3.9%），裂缝更频繁穿越层理，形成阶梯状三维连通网络。致密砂岩分形维数增幅<2%，复杂性改善有限。
**2.4 基于CT的压裂效果评价**
采用FCI综合分形维数、裂缝浓度、裂缝高度和开度变异系数。流体顺序：SC-CO? > 液态CO? > 水。Q9页岩SC-CO?压裂分形维数2.41、裂缝密度12.3 fractures/cm3、FCI 8.92；水压裂分别为2.12、7.2 fractures/cm3和3.96。岩性方面，页岩FCI显著高于致密砂岩：SC-CO?下Q9页岩FCI 8.92，Q1页岩4.43，致密砂岩仅3.40，甚至低于水压裂页岩的3.96。致密砂岩裂缝孤立、短直，缺乏分支。
**2.5 岩性对裂缝复杂性的定量对比**
补充指标：分支角、曲折度和连通率。Q9页岩分支角达±75°，曲折度1.48，连通率0.86；Q1页岩分别±15–±40°、1.31、0.74；致密砂岩<±10°、1.12、0.61。三者一致表明弱面赋予页岩更高裂缝扩展自由度，岩性结构是控制裂缝网络复杂性的根本因素。
**3.2 不同相态流体分析**
数值模拟与实验一致：SC-CO?破裂压力最低，液态CO?次之，水最高。Q9页岩模拟值SC-CO?比液态CO?低2.74 MPa（7.69%），比水低3.92 MPa（10.64%）。形态上，水压裂形成对称双翼主裂缝，SC-CO?压裂诱导多条不规则裂缝，多分支且不规则扩展。
**3.3.1 注入速率影响**
速率从10升至40 mL/min时，损伤面积比从1.79%增至2.90%，但分形维数仅从1.9224变至1.9613（变化1.98%），表明注入速率主要扩大压裂规模，对几何复杂性提升有限。
**3.3.2 注入温度影响**
温度从40°C升至70°C，损伤面积比从1.95%增至2.21%，分形维数从1.9431增至2.0120，增幅分别为11.76%和3.42%，温度通过优化物性（密度、黏度降低）产生正向但中等的作用。
**3.3.3 水平应力差影响**
应力差从0 MPa升至4 MPa，损伤面积比从1.75%增至2.31%，分形维数从1.9230增至1.9901。非对称应力场下裂缝路径更曲折，激活更多起裂点，促进多向扩展，显著提升三维连通性和复杂性。
**讨论与总结**
讨论部分指出，SC-CO?压裂过程受流体性质、渗流和流-岩相互作用的多尺度耦合控制。低黏度SC-CO?改善低渗储层渗流效率，均匀化孔隙压力场，降低局部应力集中，增加偏转和分支概率。此外，SC-CO?对有机质和碳酸盐矿物的溶解弱化作用进一步促进裂缝扩展。岩性差异通过弱面发育、脆-延性响应决定起裂阈值和扩展模式，FCI定量证实Q9页岩的弱面系统使其裂缝网络最复杂。多元线性回归量化了主控因素贡献：对裂缝复杂性，岩性（41.6%）主导，其次为水平应力差（29.5%，负向约束）、注入速率（17.1%）和温度（11.8%）；对破裂压力，水平应力差（38.2%）为主控，脆性指数（25.5%）、注入速率（24.5%）和温度（11.8%，负向）为次要。协同控制机制概括为“岩性控复杂性，应力控方向，工程参数控规模与扩展动力学”。研究结论强调：岩性结构（尤其是弱面发育）对裂缝网络复杂性起决定性作用，超过脆性指数；SC-CO?显著降低破裂压力并增强复杂性；操作参数呈岩性依赖效应；多元线性回归定量揭示了各因素贡献权重；对于强非均质页岩，推荐中等注入速率与SC-CO?结合以最大化网络复杂性；对于均质致密砂岩，需要更高注入速率补偿弱面缺失。
