1. Geomechanical In-Situ Closure Stress on Proppant Packs
The effective closure stress (\(\sigma_{\text{closure}}\)) acting on proppant grains in the fracture plane is the net compressive stress driving proppant crushing and embedment:
\[ \sigma_{\text{closure}} = \sigma_{h,\text{min}} - P_{\text{pore}} \]
During hydrocarbon production, as the well is drawn down from initial reservoir pressure (\(P_{\text{res}}\)) to flowing bottomhole pressure (\(P_{\text{wf}}\)), the effective stress on the proppant pack increases dramatically to its peak flowing value:
\[ \sigma_{\text{closure, flowing}} = \sigma_{h,\text{min}} - P_{\text{wf}} \]
where minimum horizontal stress is evaluated from poroelastic theory: \(\sigma_{h,\text{min}} = \frac{\nu}{1 - \nu}(\sigma_v - \alpha P_p) + \alpha P_p\).
2. API RP 19D Long-Term Permeability Degradation
Under severe closure stress, proppant grains experience crushing into fine particles that plug the pore throats. Per API RP 19D baseline testing, retained permeability is modeled by exponential stress degradation:
\[ k_f(\sigma) = k_0 \cdot \exp(-C_{\text{crush}} \cdot \sigma_{\text{closure}}) \cdot \eta_{\text{gel}} \]
where \(k_0\) is reference permeability at zero stress, \(C_{\text{crush}}\) is the material-specific crushing coefficient, and \(\eta_{\text{gel}}\) is the fractional cleanup efficiency accounting for residual fracturing fluid polymers.
3. Non-Darcy Inertial Choking (Cooke Beta Factor)
Near the wellbore, convergent high-velocity gas flow generates severe turbulent inertial resistance governed by the Forchheimer equation:
\[ -\frac{dP}{dx} = \frac{\mu v}{k_f} + \beta \rho v^2 \]
Cooke (1973) experimentally correlated the \(\beta\)-factor for proppant packs as \(\beta = \frac{b}{(k_f)^a}\). The effective apparent permeability is reduced by the Reynolds number term:
\[ k_{f,\text{eff}} = \frac{k_f}{1 + \frac{\beta k_f \rho v}{\mu}} \]
4. Cinco-Ley & Samaniego Dimensionless Fracture Conductivity (\(F_{CD}\))
The dimensionless ratio comparing fracture deliverability to reservoir deliverability is:
\[ F_{CD} = \frac{k_{f,\text{eff}} \cdot w}{k \cdot x_f} \]
The equivalent pseudo-skin factor (\(s_f\)) is derived from Cinco-Ley's polynomial solution:
\[ s_f = \frac{1.65 - 0.328 \ln F_{CD} + 0.116 (\ln F_{CD})^2}{1 + 0.18 \ln F_{CD} + 0.034 (\ln F_{CD})^2} - \ln\left(\frac{x_f}{r_w}\right) \]
delivering an equivalent wellbore radius \(r_w' = r_w e^{-s_f}\) and production fold of increase \(\text{FOI} = \frac{\ln(r_e / r_w)}{\ln(r_e / r_w) + s_f}\).