Abstract We present a new improvement to the Alekseev inverse scattering method. This improved inverse scattering method is extended to a double form, followed by the generation of some new solutions of the double-complex Kinnersley equations. As the double-complex function method contains the Kramer-Neugebauer substitution and analytic continuation, a pair of real gravitation soliton solutions of the Einstein’s field equations can be obtained from a double N-soliton solution. In the case of the flat Minkowski space background solution, the general formulas of the new solutions are presented.