Positive numerical integration methods for chemical kinetic systems