Prediction of Liquid CO₂ Phase-Change Fracture Radius in Coal Seams

Definition and Fundamental Principles

The prediction of liquid CO₂ phase-change fracture radius in coal seams refers to the quantitative determination of the effective fracture zone generated when compressed liquid carbon dioxide is injected into coal seams and undergoes rapid phase transition from liquid to supercritical or gaseous state. This phase change produces a dramatic volume expansion—liquid CO₂ expands approximately 400 to 500 times upon transition to gaseous state at reservoir conditions—generating localized pressures exceeding 100 MPa that exceed the tensile and shear strength of coal rock, inducing radial fractures extending from the injection point.

The fundamental physics governing this process involves three coupled phenomena:

• Thermodynamic phase transition: Liquid CO₂ (injected at pressures typically 8–15 MPa) crosses its critical point (31.1°C, 7.38 MPa) and expands rapidly, creating a pressure pulse within the coal matrix.
• Stress wave propagation: The sudden pressure release generates compressive stress waves that propagate radially through the coal body, transitioning to tensile hoop stresses at the fracture tip.
• Fracture mechanics: When the induced tensile stress exceeds the coal's fracture toughness (K_Ic), crack initiation and propagation occur, defining the fracture radius.

The fracture radius (R_f) is the critical output parameter that determines the effective stimulated reservoir volume (ESRV) and directly governs coalbed methane (CBM) extraction efficiency, coal permeability enhancement, and gas drainage capacity.

Category and Business Positioning

Within the operational framework of Cladding Technology Shanxi Co., Ltd., this research capability occupies a strategic position at the intersection of energy extraction technology and advanced materials engineering. The company's three core technology routes—TIG/MIG weld overlay, hydraulic explosive bonding, and explosion welding—are fundamentally concerned with controlled energy application to achieve metallurgical bonding, structural modification, or material enhancement. The CO₂ phase-change fracturing prediction research extends this philosophy of controlled energy deployment into the geomechanical domain.

This capability is positioned as follows:

• Upstream knowledge transfer: The fracture mechanics and stress wave propagation models developed for CO₂ fracturing prediction share fundamental analytical frameworks with those used in hydraulic explosive bonding and explosion welding process design.
• Cross-sector qualification: Demonstrates the company's capacity for multi-physics modeling and predictive engineering analysis, enhancing credibility when bidding for complex overlay and cladding projects requiring sophisticated process simulation.
• Integrated service offering: Enables the company to provide end-to-end solutions for mining and energy customers—predicting fracture geometry, then applying protective cladding or overlay materials to equipment subjected to CO₂ fracturing operations.

Technical Purpose and Value

Primary Technical Objectives

1. Fracture geometry optimization: Predict the radial extent of fractures to ensure sufficient coal permeability enhancement without creating excessive fracture networks that compromise reservoir integrity or lead to water inrush.
2. Injection parameter calibration: Establish quantitative relationships between injection pressure, CO₂ mass per stage, coal mechanical properties, and resulting fracture radius to enable precise operational control.
3. Equipment design input: Provide fracture radius data to engineers designing injection equipment, wellbore completions, and surface facilities that must withstand or accommodate the predicted fracture geometry.
4. Safety assessment: Determine the maximum fracture extent to ensure it does not intersect adjacent workings, voids, or structural boundaries, preventing gas outbursts or roof collapse.

Engineering Value

Predictive capability transforms CO₂ phase-change fracturing from an empirical operation into an engineered process. Without accurate fracture radius prediction, operators face:

• Under-stimulation: Insufficient fracture radius results in inadequate permeability enhancement and poor gas recovery rates.
• Over-stimulation: Excessive fractures compromise coal seam integrity, increase water ingress risk, and may trigger unexpected gas outbursts.
• Unpredictable wellbore damage: Fractures propagating toward the wellbore can cause casing damage, packer failure, or wellbore instability.

Key Process and Implementation Points

Prediction Methodology Framework

The fracture radius prediction methodology integrates three analytical approaches:

Methodology	Principle	Key Parameters	Accuracy Range
Elastic stress analysis	Hoop stress concentration around injection point using Kirsch solution adapted for coal anisotropy	Injection pressure (Pi), coal Young's modulus (E), Poisson's ratio (ν), in-situ stress (σh, σH, σv)	±15–25%
Fracture mechanics (LEFM)	Energy balance between released strain energy and surface energy; KI = KIc criterion	Fracture toughness (KIc), gas pressure at fracture tip (Pf), fracture surface energy (γ)	±10–20%
Numerical simulation (FEM/CFD)	Coupled thermo-fluid-structural analysis of CO₂ expansion, heat transfer, and stress redistribution	CO₂ injection rate, reservoir temperature, coal permeability, porosity, adsorption isotherm	±5–15%

Key Input Parameters and Data Requirements

Parameter Category	Specific Parameters	Typical Coal Seam Values	Measurement Method
Injection conditions	Injection pressure, CO₂ mass/stage, injection duration	8–15 MPa; 50–200 kg/stage; 30–120 min	Surface pressure monitoring, mass flow measurement
Coal mechanical properties	Compressive strength, tensile strength, Young's modulus, Poisson's ratio, KIc	σc: 10–40 MPa; σt: 1–5 MPa; E: 2–8 GPa; ν: 0.15–0.35; KIc: 0.5–2.0 MPa·m1/2	Uniaxial compression, Brazilian test, SEM/CT, fracture mechanics tests
In-situ stress field	Vertical stress, maximum/minimum horizontal stress, stress ratio	σv: 0.03–0.05 MPa/m; σH/σh: 1.2–2.5	Borehole breakout, hydraulic fracturing, microseismic analysis
Reservoir characteristics	Temperature, pressure, permeability, porosity, gas content, adsorption capacity	T: 25–45°C; k: 0.1–5 mD; φ: 5–15%; Vads: 3–12 m³/t	Well logging, Darcy flow tests, volumetric analysis, Langmuir isotherm

Core Prediction Equation

The simplified analytical expression for fracture radius derived from the energy balance approach:

R_f = √[(P_i − P_f) × V_CO₂ × β / (2π × K_Ic²)]

Where:

• R_f = fracture radius (m)
• P_i = injection pressure (MPa)
• P_f = fracture tip pressure (MPa)
• V_CO₂ = volume of CO₂ injected per stage (m³)
• β = expansion factor (dimensionless, typically 300–500)
• K_Ic = coal fracture toughness (MPa·m^1/2)

Implementation Workflow

1. Phase 1 – Data Acquisition: Collect coal core samples from target seam; conduct mechanical testing per GB/T 23561.1 (compressive strength) and fracture toughness testing per ASTM E399 methodology adapted for coal.
2. Phase 2 – Baseline Modeling: Build analytical model using measured parameters; establish initial fracture radius prediction.
3. Phase 3 – Numerical Validation: Develop finite element model (ANSYS/Abaqus) incorporating CO₂ phase behavior equation of state (Peng-Robinson EOS); simulate injection transient response.
4. Phase 4 – Sensitivity Analysis: Vary key parameters (±20%) to identify dominant factors and establish prediction confidence intervals.
5. Phase 5 – Field Calibration: Compare model predictions with microseismic monitoring data and post-fracturing permeability measurements (flow tests, pressure transient analysis).
6. Phase 6 – Iterative Refinement: Update model parameters based on field data; converge to validated prediction within ±15% accuracy.

Applicable Standards and Acceptance Criteria

Governing Standards

Domain	Standard Number	Relevance
Coal mechanical testing	GB/T 23561.1–2010	Compressive strength determination of coal
Coal mechanical testing	GB/T 23561.2–2010	Shear strength determination of coal
Fracture mechanics	ASTM E399–17	Plane-strain fracture toughness testing methodology
Fracture mechanics	ISO 12118:2012	Metallic materials – fracture toughness – crack propagation resistance
Coalbed methane	SY/T 5487–2009	Coalbed methane reservoir evaluation methods
Coalbed methane	GB/T 26212–2010	Coalbed methane geological exploration and development
CO₂ handling	GB/T 19141–2008	Industrial liquid carbon dioxide specifications
Pressure equipment	TSG 21–2016	Supervision of periodic inspection for fixed pressure vessels
Reservoir simulation	API RP 91–1989	Reservoir simulation guidelines
Geomechanical analysis	ISO 18425:2016	Petroleum and natural gas industries – geomechanical assessment

Acceptance Criteria for Prediction Deliverables

• Fracture radius prediction accuracy: Model prediction within ±15% of field-measured values (validated by microseismic event distribution or post-fracturing permeability tests).
• Parameter sensitivity documentation: Complete sensitivity analysis identifying top-3 governing parameters with quantified influence ranges.
• Uncertainty quantification: Monte Carlo simulation or probabilistic analysis providing P10/P50/P90 fracture radius distributions.
• Reproducibility: Independent verification by second analyst achieving agreement within ±10% on identical input data.
• Safety margin: Predicted maximum fracture radius must not intersect existing voids, workings, or structural boundaries within a minimum safety distance of 30 m (or as specified by mine design).

Common Risks and Controls

Risk Category	Description	Mitigation Measures
Parameter uncertainty	In-situ stress and coal properties vary laterally and vertically; laboratory samples may not represent field conditions	Conduct multi-depth core sampling; use in-situ stress measurements; apply probabilistic analysis with wide input distributions
Model simplification	Analytical models assume homogeneous, isotropic coal; actual coal seams have bedding planes, cleats, and heterogeneity	Employ numerical models with layered/heterogeneous geology; incorporate pre-existing fracture networks from image logs
CO₂ behavior deviation	Actual CO₂ phase behavior may deviate from ideal EOS predictions at reservoir conditions, especially near critical point	Use validated CO₂ property databases (NIST REFPROP); conduct laboratory PVT experiments at reservoir temperature/pressure
Fracture propagation complexity	Fractures may propagate along bedding planes rather than radially; multiple fracture stages may interact	Include bedding plane properties in model; simulate multi-stage injection sequentially; validate with microseismic data
Safety - gas outburst	Excessive fracturing may connect to gas-rich zones, triggering outburst	Limit injection pressure based on predicted fracture radius vs. distance to gas-rich zones; implement real-time pressure monitoring with automatic shut-off
Safety - water inrush	Fractures may connect to aquifers or water-bearing strata	Characterize hydrogeological conditions; limit fracture radius below water-bearing layer; monitor wellbore fluid levels during and after treatment

Application Scenarios Across the Company's Three Technology Routes

Integration with TIG/MIG Weld Overlay

The CO₂ phase-change fracturing prediction research directly supports TIG/MIG weld overlay applications in the following ways:

• Protective overlay design for fracturing equipment: Prediction of fracture radius and associated stress fields informs the design of overlay thickness and composition for injection manifolds, pressure vessels, and wellhead equipment exposed to CO₂ fracturing operations. Equipment must withstand cyclic pressure loads corresponding to predicted fracture pressures.
• Corrosion-resistant cladding for CO₂ service: Liquid CO₂ at elevated pressures becomes mildly corrosive to carbon steel, particularly in the presence of trace water (forming carbonic acid). Fracture radius prediction helps determine the operational pressure envelope, which in turn defines the overlay material selection (e.g., 309L/316L stainless steel overlays per AWS D10.6 or ASME B31.3).
• WPS qualification for high-pressure equipment: The predicted maximum operating pressures from fracture modeling provide the design basis for weld procedure qualification per NB/T 47014 or ASME Section IX, ensuring overlay welds on pressure-containing equipment are qualified for actual service conditions.

Integration with Hydraulic Explosive Bonding

The fundamental physics of hydraulic explosive bonding—where controlled pressure waves drive metallurgical bonding between dissimilar materials—shares analytical frameworks with CO₂ phase-change fracture prediction:

• Pressure wave modeling transfer: The stress wave propagation models developed for CO₂ fracture prediction (solving the wave equation in heterogeneous media) are directly applicable to hydraulic explosive bonding process design, where precise pressure pulse shaping determines bond quality.
• Energy threshold determination: Similar to determining the minimum injection pressure needed to initiate coal fracture, hydraulic explosive bonding requires determination of minimum impact velocity thresholds for bonding specific material pairs. The analytical methodology transfers directly.
• Equipment qualification: Hydraulic explosive bonding systems operate at pressures comparable to CO₂ fracturing injection pressures (10–100 MPa). Fracture radius prediction research validates the pressure containment design of these systems, supporting qualification per TSG 21–2016 and GB 150.

Integration with Explosion Welding

Explosion welding involves detonation-driven collision of material surfaces, generating extreme pressures and velocities. The connection to CO₂ fracture prediction research includes:

• Detonation pressure calibration: Understanding the relationship between energy input and resulting fracture/propagation geometry in coal informs analogous calculations for explosion welding charge design, where detonation pressure determines collision velocity and bonding quality per ASTM A377 or ISO 18754.
• Fracture propagation control: In explosion welding, the fracture mechanics of the collision interface (spall formation, jetting) parallels the fracture mechanics of CO₂-induced coal fracturing. Predictive models for one application enhance understanding of the other.
• Clad plate/pipe specification: Explosion-welded clad plates used in CO₂ fracturing applications (e.g., for high-pressure CO₂ storage tanks, injection cylinders) benefit from fracture radius prediction data that defines the maximum design pressure, enabling proper selection of base/overlay material combinations per ASME SA-240 or NACE MR0175.

Contribution to Qualification Building, Product Delivery, and Customer Value

Qualification Building

1. Multi-physics modeling credential: Demonstrating validated CO₂ fracture prediction capability establishes the company as a provider of advanced engineering analysis, strengthening bids for complex overlay and cladding projects requiring process simulation and predictive design.
2. Cross-disciplinary expertise validation: The research bridges materials science, geomechanics, and process engineering—credentials that support qualification for integrated projects where cladding technology must be designed for specific operational environments (e.g., equipment in CO₂ fracturing service).
3. Standards compliance track record: Engagement with GB/T 23561, ASTM E399, SY/T 5487, and related standards builds a documented track record of standards-based engineering, directly supporting qualification applications under ISO 9001, ISO 3834, and NB/T 47014 frameworks.

Product Delivery Enhancement

• Design-optimized overlay specifications: Fracture radius predictions provide precise operational boundary conditions for overlay design, enabling right-sizing of weld overlay thickness (avoiding over-specification that increases cost, or under-specification that risks failure).
• Performance guarantee basis: Validated prediction models provide the technical basis for performance guarantees on cladded equipment used in fracturing operations, supporting customer confidence and contract compliance.
• Accelerated WPS qualification: Accurate prediction of service conditions reduces the number of qualification trials needed for weld procedures, as design parameters are well-defined rather than conservatively estimated.

Customer Value Creation

• Reduced operational risk: Customers deploying CO₂ fracturing operations gain confidence that equipment (clad/overlaid) is designed for actual predicted conditions rather than generic worst-case assumptions, reducing unnecessary conservatism and cost.
• Integrated service package: The ability to offer fracture prediction alongside equipment cladding/overlay creates a differentiated value proposition—customers receive both process optimization and equipment protection from a single qualified supplier.
• Compliance and safety assurance: Predictive modeling supports regulatory compliance by demonstrating that fracture geometry is controlled within safe boundaries, protecting both customer operations and the company's reputation for quality engineering.
• Lifecycle cost optimization: Accurate fracture predictions enable proper selection of overlay material and thickness, optimizing the balance between initial capital cost and long-term maintenance/replacement intervals.

Summary

The prediction of liquid CO₂ phase-change fracture radius in coal seams represents a sophisticated engineering capability that, while originating in the energy extraction domain, provides substantial value across Cladding Technology Shanxi Co., Ltd.'s full technology portfolio. The predictive modeling frameworks, standards compliance practices, and risk management methodologies developed through this research transfer directly to TIG/MIG weld overlay design, hydraulic explosive bonding process optimization, and explosion welding charge design. This cross-pollination of technical knowledge strengthens the company's qualification position, enhances product delivery precision, and creates integrated value propositions that differentiate the company in competitive markets.