CO₂相变致裂岩石损伤范围计算及影响因素研究技术解析
一、技术定义与基本原理
CO₂相变致裂技术（CO₂ Phase-Change Fracturing）是一种以液态二氧化碳（Liquid CO₂）或超临界二氧化碳（Supercritical CO₂）作为致裂工作介质，利用其相变膨胀产生的巨大压力（通常可达数十兆帕至数百兆帕量级）对岩石进行水力压裂替代性致裂的非常规增产改造技术。与常规水力压裂不同，该技术不依赖大量水基压裂液，而是利用CO₂从液态向气态转变时体积膨胀约数百倍的物理特性，在储层岩石中产生并扩展裂缝网络，从而实现储层改造目的。
《CO₂相变致裂岩石的损伤范围的计算及影响因素研究》这一课题的核心在于：建立CO₂相变致裂过程中岩石损伤场（Damage Field）的定量计算模型，系统研究影响损伤范围（Damage Range）的关键参数及其耦合机制。所谓"损伤范围"，指的是在CO₂相变膨胀压力作用下，岩石产生微裂隙萌生、扩展、连通并形成宏观裂缝网络所覆盖的空间区域。该研究属于岩石力学与断裂力学的交叉前沿领域，涉及多物理场耦合（热-力-化）的复杂本构关系。
CO₂相变致裂的基本物理过程可归纳为以下几个阶段：
- 注入阶段：液态CO₂在高压下通过井筒注入储层，此时CO₂处于高压液态或超临界态（温度高于31.1°C，压力高于7.38MPa时为超临界态），密度高、流动性好。
- 相变膨胀阶段：CO₂进入储层后，由于储层温度通常高于CO₂临界温度，液态CO₂迅速吸收地层热量发生相变（沸腾），体积急剧膨胀，产生巨大的膨胀压力作用于周围岩石。
- 裂缝萌生与扩展阶段：当相变膨胀应力超过岩石的抗张强度（Tensile Strength）时，岩石中萌生微裂隙；随着压力持续作用，微裂隙扩展、连通，形成宏观裂缝。
- 损伤场形成阶段：裂缝扩展带动周围岩石产生应力集中和微损伤区，形成以裂缝尖端为中心的损伤范围。
二、所属大类与业务定位
该技术条目属于非常规油气/地热储层改造工艺技术范畴，在科雷丁技术（山西）有限公司的业务体系中，其定位具有双重属性：
2.1 技术研发支撑定位
CO₂相变致裂损伤范围计算研究属于公司基础理论研究与工程技术储备板块。该研究虽不直接涉及金属复合材料制造，但为公司在以下领域提供了关键的技术输入和工程依据：
- 明确CO₂致裂作业中井下管柱、井口装置、地面管线所承受的载荷特征与服役工况，为复合管/复合板的设计选材提供边界条件；
- 揭示CO₂（尤其是超临界CO₂）与金属材料的腐蚀-应力耦合损伤机制，指导耐蚀堆焊过渡层与功能层的选材与工艺设计；
- 为公司在CO₂驱油（CO₂-EOR）、CO₂封存（CCUS）、地热开发等新兴领域的产品配套提供技术依据。
2.2 产品配套与工程服务定位
CO₂相变致裂作业对井下设备提出了极为苛刻的要求。超临界CO₂对碳钢具有显著的腐蚀-应力开裂（CO₂ Corrosion / Stress Corrosion Cracking）风险，尤其在高压、含微量水及硫化物的条件下，碳钢及低合金钢的腐蚀速率可急剧升高。这直接催生了对以下产品的刚性需求：
- 耐蚀复合油管/套管（如碳钢基体+不锈钢堆焊内衬管）；
- 耐蚀复合井口装置用复合板（如碳钢+316L/2205双相不锈钢复合板）；
- 耐蚀法兰、阀门内衬等关键承压部件的堆焊修复与制造。
三、技术目的与核心价值
3.1 技术目的
本研究旨在实现以下技术目标：
- 建立定量计算模型：基于断裂力学与损伤力学理论，建立CO₂相变致裂过程中岩石损伤范围的解析或数值计算模型，实现损伤范围的可预测性；
- 识别关键影响因素：系统辨识并量化分析影响损伤范围的主要参数，包括注入压力、注入速率、储层温度、储层压力、岩石物理力学参数（杨氏模量、泊松比、抗拉强度、断裂韧性）、裂缝初始形态等；
- 优化致裂参数：基于损伤范围计算结果，为现场CO₂相变致裂作业提供最优参数窗口建议，实现"可控致裂"；
- 保障作业安全：通过精确预测损伤范围，避免因过度致裂导致井壁失稳、套管损坏或储层破坏等安全事故。
3.2 对公司核心价值
该研究成果对公司具有多维度的战略价值：
- 资质建设价值：参与CCUS（Carbon Capture, Utilization and Storage）、CO₂-EOR等国家重点领域的技术研发，有助于公司积累非常规油气开发、碳封存等新兴领域的技术资质，拓宽业务范围；
- 产品交付价值：损伤范围计算结果直接关联井下设备的设计载荷与选材要求，为公司提供精准的耐蚀复合产品技术规范；
- 客户价值：为油气田客户和地热开发客户提供从储层改造到井下设备配套的一体化技术方案，提升公司综合服务能力与客户粘性。
四、关键计算方法与实施要点
4.1 损伤范围计算模型
CO₂相变致裂岩石损伤范围的计算通常采用以下理论框架：
4.1.1 应力场分析
在CO₂相变膨胀压力下，储层岩石中的应力状态可用弹性力学理论描述。对于径向对称的裂缝扩展问题，可采用Kirsch解或Bunger解进行应力分布计算：
在远场应力（σ_h, σ_H, σ_v）和近场相变膨胀压力（P_CO₂）共同作用下，裂缝尖端附近的应力集中系数（Stress Concentration Factor, SCF）可用以下公式估算：
K_t = 1 + 2(a/r)，其中a为裂缝半长，r为距裂缝尖端的距离。
当局部应力超过岩石抗拉强度时，该区域即进入损伤状态。
4.1.2 损伤力学模型
基于连续损伤力学（Continuum Damage Mechanics, CDM）理论，岩石损伤变量D（Damage Variable）定义为：
D = 1 - σ_eff/σ_nom，其中σ_eff为有效应力，σ_nom为名义应力。
损伤变量D的演化方程通常采用Kachanov-Rabotnov型演化律：
dD/dt = (1/Ω) · (σ_eff/σ_0)^n · (1/ψ(D))
其中Ω、n、ψ(D)为材料损伤参数，σ_0为参考应力。
4.1.3 断裂力学判据
裂缝扩展的判据采用应力强度因子准则（Stress Intensity Factor Criterion）：
K_I ≥ K_IC 时裂缝扩展
其中K_I为I型应力强度因子，K_IC为岩石的断裂韧性（Mode I Fracture Toughness）。对于CO₂相变致裂产生的径向裂缝，K_I的计算需考虑相变膨胀压力的分布特征。
4.1.4 数值模拟方法
对于复杂地质条件下的损伤范围计算，通常采用有限元法（FEM）或离散元法（DEM）进行数值模拟：
- 有限元法（FEM）：适用于连续介质假设下的损伤场分析，可耦合热-力-化多物理场；
- 离散元法（DEM）：适用于模拟岩石的非连续变形与破碎过程，能更真实地反映裂缝的萌生与扩展；
- XFEM（扩展有限元法）：无需在网格中追踪裂缝路径，适用于动态裂缝扩展模拟。
4.2 关键影响因素分析
影响CO₂相变致裂岩石损伤范围的因素可分为以下四大类：
| 因素类别 | 具体参数 | 影响机制 | 对损伤范围的影响趋势 |
|---|---|---|---|
| 注入参数 | 注入压力（P_inj） | 决定相变膨胀驱动力大小 | 压力↑ → 损伤范围↑ |
| 注入速率（q_inj） | 影响压力传播波速与压力场分布 | 速率↑ → 压力集中↑ → 损伤范围可能↓（压力来不及扩散） | |
| 注入量（V_inj） | 决定总膨胀功与裂缝总体规模 | 注入量↑ → 损伤范围↑ | |
| 储层条件 | 储层温度（T_res） | 影响CO₂相变速率与膨胀比 | 温度↑ → 相变速率↑ → 损伤范围↑（但过度升温可能导致CO₂超临界态流动性增强，压力分散） |
| 储层压力（P_res） | 与注入压力形成压差，驱动裂缝扩展 | 压差↑ → 损伤范围↑ | |
| 地应力（σ_h, σ_H, σ_v） | 控制裂缝走向与扩展方向 | 最小主应力方向控制裂缝走向；应力差↑ → 裂缝更定向 → 损伤范围在垂直方向可能↓ | |
| 岩石力学参数 | 杨氏模量（E） | 反映岩石刚度 | E↑ → 应力集中↑ → 局部损伤↑但损伤范围可能↓ |
| 泊松比（ν） | 影响应力场分布 | ν↑ → 横向约束↑ → 损伤范围可能↓ | |
| 抗拉强度（σ_t） | 决定裂缝萌生的临界条件 | σ_t↑ → 萌生门槛↑ → 损伤范围↓ | |
| 断裂韧性（K_IC） | 控制裂缝扩展的临界条件 | K_IC↑ → 扩展门槛↑ → 损伤范围↓ | |
| CO₂物性 | 相变膨胀比（β） | 液态→气态体积膨胀倍数 | β↑ → 膨胀功↑ → 损伤范围↑ |
| CO₂纯度与含水率 | 影响相变温度与腐蚀行为 | 含水率↑ → 相变温度↑ → 损伤范围↓；但腐蚀风险↑ |
4.3 典型计算案例参数
| 参数 | 典型取值范围 | 单位 | 备注 |
|---|---|---|---|
| 注入压力 | 20-60 | MPa | 取决于储层深度与破裂压力 |
| 注入速率 | 0.5-5.0 | m³/min | 受泵送能力与管柱承压限制 |
| 储层温度 | 60-180 | °C | 地热储层可达200°C以上 |
| 储层压力 | 10-50 | MPa | 与埋深相关 |
| 岩石杨氏模量 | 5-40 | GPa | 砂岩10-30 GPa，碳酸盐岩20-40 GPa |
| 岩石抗拉强度 | 5-30 | MPa | — |
| 断裂韧性K_IC | 0.5-3.0 | MPa·m^(1/2) | — |
| 损伤范围（预测） | 1-50 | m | 取决于上述所有参数的耦合 |
五、执行标准与验收依据
5.1 相关技术标准
- GB/T 23561-2009《岩石力学参数试验方法》——岩石物理力学参数测试依据；
- GB/T 22374-2008《岩石断裂韧性试验方法》——断裂韧性K_IC测试依据；
- SY/T 6610-2016《压裂改造设计方法》——压裂/致裂设计参考标准；
- API RP 90《Production Logging Interpretation》——生产测井解释参考；
- ISO 17154-1:2016《Petroleum and natural gas industries — Determination of rock fracture toughness》——岩石断裂韧性测定国际标准；
- ASTM E399-21《Standard Test Method for Linear-Elastic Plane-Strain Fracture Toughness K_IC》——材料断裂韧性测试方法（适用于人工岩样）；
- NACE MR0175/ISO 15156系列标准——含CO₂环境中抗蚀材料选材依据（与井下设备配套直接相关）；
- GB/T 13296-2013《无缝不锈钢管》——耐蚀复合管内衬管选材依据；
- NB/T 20318-2018《承压设备用复合钢板》——复合板制造依据；
- ASME B31.3《Process Piping》——工艺管道设计（含CO₂管线）；
- ASME BPV Code Section VIII Div.1——压力容器设计（含CO₂致裂地面设备）。
5.2 验收依据
该研究课题的验收通常包括以下维度：
- 理论模型验证：计算模型需通过与室内物理模拟实验（如CO₂相变致裂岩石模拟实验）或现场先导试验数据进行对比验证，损伤范围预测值与实测值的偏差应控制在±20%以内；
- 敏感性分析：需完成对各影响因素的敏感性排序（Sensitivity Ranking），明确主控因素；
- 工程应用报告：形成可指导现场作业的参数优化建议与操作窗口。
六、常见风险与控制措施
| 风险类别 | 具体风险描述 | 风险等级 | 控制措施 |
|---|---|---|---|
| 井控风险 | CO₂相变膨胀压力失控导致井喷或井壁失稳 | 高 | 精确计算破裂压力与损伤范围，设置安全注入压力上限；配备防喷器（BOP）系统；实施实时压力监测 |
| 套管损坏 | 过度致裂导致套管受挤变形或破裂 | 高 | 损伤范围计算中纳入套管位置与承压能力；控制裂缝扩展范围不超过套管外壁安全距离 |
| CO₂腐蚀 | 超临界CO₂或含CO₂地层水对井下管柱、井口装置造成严重腐蚀 | 高 | 选用符合NACE MR0175/ISO 15156标准的耐蚀材料；采用不锈钢堆焊内衬复合管；实施缓蚀剂注入 |
| 泄漏风险 | CO₂高压泄漏造成人员窒息或环境污染 | 中 | 地面设备按ASME BPV Code设计制造；设置CO₂浓度检测报警系统；制定应急预案 |
| 计算偏差 | 损伤范围计算模型简化导致预测偏差过大 | 中 | 采用多模型交叉验证（解析解+数值模拟+物理实验）；开展不确定性分析（Uncertainty Analysis）；设置安全裕度 |
| 储层破坏 | 裂缝网络过度扩展导致储层连通性破坏 | 中 | 基于损伤范围计算优化注入参数；实施分段控制注入；后期进行生产测井评估裂缝形态 |
七、在公司三大技术路线中的应用场景
7.1 TIG/MIG堆焊技术路线
CO₂相变致裂损伤范围研究为TIG/MIG堆焊技术提供了明确的选材依据与工艺参数设计输入：
- 耐蚀内衬管堆焊：CO₂相变致裂作业中，井下油管/套管长期承受高压CO₂环境。研究明确了CO₂的相态、温度、压力、含水率等工况参数，为堆焊过渡层（如309L）与功能层（如316L、321、2205双相不锈钢）的选材提供边界条件。例如，当储层温度超过150°C且CO₂含水率较高时，需选用耐氯离子应力腐蚀开裂（Cl-SCC）性能更优的合金体系。
- 堆焊层厚度设计：基于CO₂腐蚀速率预测（可依据NORSOK M-506或API RP 945），结合损伤范围计算中确定的服役周期与工况波动，确定堆焊功能层的最小厚度与堆焊道数。典型设计参数如下：
| 设计参数 | 典型取值 | 依据 |
|---|---|---|
| 过渡层合金 | 309L (UNS S30908) | 匹配碳钢基体与奥氏体功能层的膨胀系数与碳含量 |
| 功能层合金 | 316L (UNS S31603) / 2205 (UNS S32205) | 满足NACE MR0175/ISO 15156对CO₂环境的耐蚀要求 |
| 堆焊层总厚度 | 1.5-3.0 mm | 基于腐蚀速率预测+安全裕度（通常取10年服役周期） |
| 焊接方法 | TIG打底 + MIG盖面 | 保证熔合质量与堆焊效率 |
| 层间温度控制 | ≤150°C | 防止热影响区（HAZ）敏化与晶间腐蚀 |
| 焊后热处理 | 300-350°C × 2h 去应力退火 | 消除残余应力，降低CO₂应力腐蚀开裂风险 |
7.2 水压复合技术路线
水压复合（Hydrostatic Expansion / Hydraulic Expansion Bonding）技术在CO₂相变致裂领域的应用主要体现在地面高压容器与管路的制造：
- CO₂高压储液罐：CO₂相变致裂前需要将CO₂液化并高压储存。CO₂液化储罐工作压力通常达15-20MPa，且CO₂对碳钢有腐蚀风险。采用水压复合技术制造碳钢基体+不锈钢内衬的复合储罐，既保证承压强度，又满足内表面耐蚀要求。复合板的结合强度需满足NB/T 20318-2018要求，水压复合压力通常为设计压力的1.5-2.0倍。
- 高压注CO₂管线：从CO₂储液罐到井口的输送管线承受高压CO₂，采用水压复合管（如碳钢+316L复合管）可有效解决内表面腐蚀问题。水压复合工艺参数（复合压力、保压时间、加热温度）需根据管径、壁厚、材料组合进行专项评定。
- 井口装置壳体：井口装置（Christmas Tree）的壳体、阀门体等承压部件，可采用碳钢+不锈钢复合板制造，利用水压复合工艺实现冶金结合，满足ASME B31.3与API 6A的相关要求。
7.3 爆炸焊复合技术路线
爆炸焊复合（Explosive Cladding / Explosive Welding）技术在CO₂相变致裂领域的应用主要集中在大尺寸厚壁耐蚀复合构件的制造：
- 大型CO₂储罐/反应器复合板：对于大型CO₂液化储罐或CO₂-EOR地面处理装置，爆炸焊复合可高效制造大面积、厚壁（基体板厚可达50-100mm）的碳钢+不锈钢复合板。爆炸焊复合的界面结合强度高（通常超过母材强度），且不受材料厚度比的限制，适合制造厚壁耐蚀构件。
- 井下封隔器/配水器外壳：CO₂相变致裂作业中使用的井下工具（如封隔器、配水器）外壳需承受高压CO₂环境。采用爆炸焊复合技术制造碳钢+哈氏合金（Hastelloy C-276）或碳钢+2205双相不锈钢复合管，可满足极端腐蚀工况下的长期服役要求。
- 复合管接头/异径管：井下管柱中的异径接头、弯头等异形件，可采用爆炸焊复合技术制造，确保接口处的耐蚀性与机械强度。
7.4 三大技术路线对比与应用选择
| 技术路线 | 适用场景（CO₂致裂领域） | 优势 | 局限 | 典型产品 |
|---|---|---|---|---|
| TIG/MIG堆焊 | 井下油管/套管内表面耐蚀堆焊；法兰密封面堆焊；阀门内衬堆焊 | 灵活性强，可现场施工；适用于修复与制造；堆焊层厚度可控 | 堆焊层厚度有限（通常≤5mm）；效率较低；对操作者技能要求高 | 耐蚀复合油管、堆焊内衬套管、堆焊法兰 |
| 水压复合 | 地面高压储液罐、高压注CO₂管线、井口装置壳体 | 结合强度高且均匀；适合中等厚度（基体≤30mm）；工艺可控性好 | 对管径有上限（通常≤600mm）；对材料组合有选择性；设备投资较高 | 碳钢+不锈钢复合管、复合板容器壳体 |
| 爆炸焊复合 | 大型储罐/反应器复合板、厚壁井下工具外壳、异形复合件 | 适合厚壁（基体≥30mm）大面积复合；结合强度高；材料组合范围宽 | 安全性管理要求高；设备投资大；异形件工艺难度大；批次一致性控制要求高 | 厚壁复合板、大型复合容器、复合异径管 |
八、对公司资质建设、产品交付与客户价值的综合作用
8.1 资质建设维度
《CO₂相变致裂岩石的损伤范围的计算及影响因素研究》这一课题的开展，对公司资质建设具有显著推动作用：
- CCUS/CCS领域技术资质：参与CO₂封存与利用领域的核心技术研究，有助于公司取得CCUS相关项目的设计、制造与服务资质；
- 非常规油气开发配套资质：CO₂相变致裂是页岩气、致密油气、页岩油等非常规油气储层改造的重要技术之一，相关研究积累有助于公司进入非常规油气开发供应链；
- 地热开发配套资质：CO₂相变致裂在地热储层增强型（EGS, Enhanced Geothermal System）开发中具有重要应用前景，相关研究有助于公司拓展地热领域业务；
- 科研合作与平台资质：该课题通常需要与高校（如中国石油大学、中国矿业大学等）或科研院所合作，有助于公司建立产学研合作平台，提升企业技术形象。
8.2 产品交付维度
研究成果直接赋能公司产品的技术规格确定与质量控制：
- 精准选材：损伤范围计算中确定的CO₂相态、温度、压力、含水率等参数，直接决定了耐蚀复合产品的合金体系选择。例如，在超临界CO₂环境下（T>31.1°C, P>7.38MPa），需选用对CO₂腐蚀有优异抵抗能力的合金（如316L、2205、C-276等）；
- 工艺参数优化：基于服役工况的精确描述，可优化堆焊工艺参数（如热输入、层间温度、焊后热处理制度），确保堆焊层的组织性能满足长期服役要求；
- 质量检验标准制定：研究成果有助于公司制定针对CO₂环境服役产品的专项检验标准（如堆焊层耐CO₂腐蚀试验方法、复合界面结合强度测试方法等），提升产品交付质量的可追溯性与可靠性。
8.3 客户价值维度
该研究为客户提供了超越单一产品供应的整体技术解决方案：
- 储层改造+设备配套一体化方案：公司可基于损伤范围计算结果，为客户提供从储层改造参数优化到井下耐蚀设备选型的一体化技术方案，降低客户的多供应商协调成本；
- 风险预警与寿命预测：基于CO₂腐蚀与应力环境分析，公司可为客户提供井下设备的腐蚀寿命预测与风险预警服务，帮助客户制定预防性维护计划，降低非计划停机风险；
- 降本增效：通过精准的选材与工艺设计，避免过度设计（Over-design），在保证安全性的前提下降低材料成本与制造成本，为客户创造经济价值。
九、总结与展望
《CO₂相变致裂岩石的损伤范围的计算及影响因素研究》虽然属于岩石力学与储层工程领域的理论研究课题，但其研究成果对科雷丁技术（山西）有限公司的核心业务——双金属及多层金属复合材料设计与堆焊制造——具有重要的工程指导意义。该研究为公司在TIG/MIG堆焊、水压复合、爆炸焊复合三大技术路线中，提供了明确的工况边界条件、选材依据与工艺参数设计输入。
随着全球"双碳"目标的推进，CO₂-EOR、CCUS、地热开发等领域将迎来快速发展期，对耐蚀复合材料的需求将显著增长。公司通过参与此类前沿技术研究，不仅能够巩固现有技术优势，更能够在新兴市场中抢占技术制高点，实现从"产品制造商"向"技术解决方案提供商"的战略转型。
核心启示：CO₂相变致裂损伤范围研究不仅是储层工程的技术课题，更是耐蚀复合材料产品设计的"前置输入"。精准的损伤范围预测意味着精准的服役工况描述，而精准的服役工况描述意味着精准的选材与工艺设计。这一技术链条的贯通，是公司从"制造驱动"向"技术驱动"升级的关键路径。