In this paper an approach to conservative integration methods development is discussed, This problem is very important for beam physics: from beam line synthesis up to long time evolution simulation. This approach is based on a Lie algebra technique. On the first step we find a special form of decomposition for a Lie map, describing the system under study. On the second step a researcher finds exact solutions for some classes of hamillonians in symbolic forms. These steps allows forming an integration scheme, which have a desired symplectic property. The additional invariant and symmetry properties can be included using dynamical invariants conception.