A Predictive Mudstone Porosity Model Which Includes the Effects of Mechanical Compaction, Temperature, Mineralogy, Chemical Diagenesis, and Overpressure

Krushin, James T.¹ (1) Consultant, Houston, TX

Traditional mudrocks compaction curves are simple exponential decay of porosity with depth beginning with a high porosity at the surface. These methods do not accurately model all the parameter controlling compaction of mudrocks. A mudrock compaction model based on desorption isotherms of pure clays and clay mixtures properly incorporates the variables controlling mechanical, thermal, and chemical compaction, as well as variation in mineralogy. Study of sorption isotherms of clays illustrates the defining variables that control the water content in mudrocks. These variables are: CEC (cation exchange capacity), the specific exchangeable cation, effective stress (mechanical compaction), and temperature. Sorption isotherms are simply a measure of the mass of water per gram of dry clay, as relative humidity (p/po) is varied between 0 and 100%. This water includes interlamellar, intraaggregate, surface adsorbed, and interaggregate. Published desorption isotherms of Na-exchanged clays, the most abundant species in the subsurface, show that the amount of water per gram of dry clay increases with increasing CEC values for a given p/po. A plot of p/po vs. mass of water per meq (a unit of CEC measurement) for selected clays and clay mixtures results in a detailed compaction curve for the depths encountered in exploration. Thermodynamic equations enable the conversion of p/po to effective overburden stress, as well as quantify the temperature effects of thermal compaction. The effects of chemical diagenesis (i.e. smectite to illite) and variable mineralogy are incorporated in the compaction model via the bulk CEC parameter, calculated from well logs using a published algorithm. Density log derived porosity and porosity from the model show excellent agreement. Application of this model results in better understanding and quantification of pore pressure prediction, drilling fluid/wellbore interaction, and seismic modeling.