The ability to numerically simulate the behavior of geologic hydrate reservoirs has improved substantially over the past 5 years in terms of both code availability and capabilities [8]. There are currently several numerical models that can simulate the system behavior in hydrate-bearing geologic media (e. g., [30,68,94,130,143,152, 153, 188]). Several of these codes were calibrated against the Mallik 2002 produc­tion test data, and the data from the Mt. Elbert MDT test [2] . A code-comparison study [1, 214] indicated that most of the participating codes appear capable of simu­lating the behavior of hydrates and reservoir fluids during common dissociation scenarios. The current consensus is that the models generally account for the impor­tant physics of the problem, and that validation and calibration (rather than ade­quacy of the numerical code capabilities) will be a constraining factor in the assessment of hydrates as an energy resource [214, 2] . Additionally, while uncer­tainties exist in the description of properties and processes involved in numerical simulators (e. g., thermal properties of composite GH-bearing systems, relative per­meability and capillary pressure, geomechanical properties related to subsidence after the dissociation of the cementing GH from the porous media, etc.), these knowledge gaps are being addressed [125, 127].