Data-driven estimation of three-phase saturation during gas hydrate depressurization using CT images