The variation in the cleat permeability of coalbed methane (CBM) reservoirs is attributed primarily to two cardinal processes, with opposing effects. Increase in effective stresses with reduction in pore pressure tends to decrease the cleat permeability, whereas the sorption-induced coal matrix shrinkage actuates reduction in the effective stresses which increases the reservoir permeability. The net effect of the two processes determines the pressure-dependent-permeability and, hence, the overall trend of CBM production with depletion. Several analytical models have been developed and used to predict the dynamic behavior of CBM reservoir permeability during production through pressure depletion, all based on combining the two effects. The purpose of this study was to introduce modifications to two most commonly used permeability models, namely the Palmer and Mansoori, and Shi and Durucan, for permeability variation and evaluate their performance when projecting gas production. The basis for the modification is the linear relationship between the volume of sorbed gas and the associated matrix shrinkage. Hence, the impact of matrix shrinkage is incorporated as a function of the amount of gas produced, or that remaining in coal, at any time during production. Since the exact production from a reservoir is known throughout its life, this significantly simplifies the process of permeability modeling. Furthermore, the modification is also expected to streamline the process of modeling by classifying the shrinkage parameters for coals of different regions, but with similar characteristics. A good analogy is the San Juan basin, where sorption characteristics of coal are so well understood and defined that operators no longer carry out laboratory sorption work. The goal is to achieve the same for incorporation of the matrix shrinkage behavior. Another modification is to incorporate the matrix, or grain, compressibility effect of coal as a correction factor in the Shi and Durucan model so as to assess the permeability variation based on the true shrinkage of coal matrix with reservoir drawdown. Finally, application of the modified models may be carried out for scenarios where the gas content of coal varies with time, either due to injection of a second gas to enhance the recovery of methane, or gas enhancing techniques, such as, bio-stimulation of coal. The original models are currently unable to handle this, particularly when the gas content of the reservoir increases. The research is aimed at simplifying and, in fact, improving the performance of the theoretical models in predicting the variation of coal reservoir permeability.

This dissertation presents contributions to the understanding of the dynamic nature of permeability of Indonesian coal. It is the first-of-its-kind study, first presenting a comparison of experimental results with those obtained using existing analytical permeability models, and then modifying the existing anisotropic model for application to the unique physical structure of Indonesian coal. The first problem addressed in this dissertation was establishing the pressure-dependentpermeability of coal in a laboratory environment replicating in situ conditions for two coal types from the Sanga Sanga basin of Kalimantan, Indonesia. The change in permeability with depletion and the corresponding volumetric strain of coal were measured in the laboratory under uniaxial strain condition (zero lateral strain). Two gases, helium and methane, were used as the flowing fluids during experimental work. The results showed that, decreasing pore pressure resulted in significant decrease in horizontal stress and increased permeability. The permeability increase at low reservoir pressure was significant, a positive finding for Indonesian coals. Using the measured volumetric changes with variations in pressure, the cleat compressibility for the two coal types was estimated. In a separate effort, volumetric strain as a result of desorption of gases was measured using sister samples under unconstrained condition, in absence of the stress effect. Sorptioninduced strain processes were modeled using the Langmuir-type model to acquire the two important shrinkage parameters. All parameters calculated using the experimental data were used for the modeling exercise. The second component of this dissertation is the permeability variation modeling to enable projecting long-term gas production in the Sanga Sanga basin. For this, two commonly used isotropic permeability models were selected. These models, developed primarily for the San Juan coal, were unable to match the measured permeability data. This was believed to be due to the inappropriate geometry used to represent Indonesian coal, where butt cleats are believed to be absent. This was followed by application of the most recent model, incorporating partial anisotropy in coal. This consideration improved the modeling results although there clearly was room for improvement. The final challenge addressed in this dissertation was to consider the coal geometry appropriate for Indonesian coal, stack of sheets as opposed to a bundle of matchsticks. In order to incorporate the structural anisotropy for the stack of sheets geometry, two input parameters were modified, based on geo-mechanical anisotropy. After applying these to the modified model, the permeability modeling results were compared with the experimental data. The matches improved significantly. Finally, the effect of maximum horizontal stress on permeability of coal was estimated by using high and low maximum horizontal stress values and constant vertical and minimum horizontal stresses. The effect of maximum horizontal stress on permeability was found to be significant under uniaxial strain condition for both coals.