High-order optimization methods for large-scale 3D CSEM data inversion