Reciprocal cross-correlation analysis of two-phase seepage processes and reservoir heterogeneities in CO₂ saline aquifer sequestration

When CO2 saline aquifer storage is carried out, the heterogeneity of reservoir rock is an important factor affecting CO2 transport, and the reservoir heterogeneity in numerical simulations is mainly manifested as the heterogeneity of the parameter field. Since the parameter distributions across the reservoir are not available with the existing probes, the stochastic finite element method is combined with a two-phase flow model to establish an unconditional random field of permeability, and computations are performed using the Monte Carlo method. The permeability, CO2 maximum migration distance (M-d) and CO2 sweep area (S-a) were analyzed for mutual correlation. The permeability correlation area affecting M-d and S-a was obtained, and the changes in the correlation area under the coefficient of variation (C-v) and correlation length (lambda(x)) of the permeability field in the different reservoirs were analyzed. The kriging superposition approach (KSA) was subsequently used to estimate both the M-d and S-a of the target reservoir by establishing conditional random fields based on the sampling parameters in regions with different correlations, resulting in errors of 0.66% for M-d and 0.96% for S-a in the high correlation region and 4.86% and 3.12% for M-d and S-a in the low correlation region, which suggested that the sampling results from the high correlation region were less biased in the estimation. Under limited sampling conditions, it is recommended that samples be collected in regions with high correlations to reduce the uncertainty of CO2 transport analysis due to unknown heterogeneity. (c) 2024 Society of Chemical Industry and John Wiley & Sons, Ltd.