We have implemented two-component spin-orbit relativistic effective core potential (SOREP) methods in an all-electron relativistic program DIRAC. This extends the capacity of the two-component SOREP method to many ground and excited state calculations in a single program. As the test cases, geometries and energies of the small halogen molecules were studied. Several two-component methods are compared by using spin-orbit and scalar relativistic effective core potentials. For the $I_2$ molecule, excitation energies of low-lying excited states agree well with those from corresponding all-electron methods. Efficiencies in SOREP calculations enhanced by using symmetries are also discussed briefly.