In this paper, we propose and analyze by means of simulations the use of surrogate data algorithms for blind detection of nonlinearities in multiple-echo ultrasonic signals. We assume a blind scheme so that no information about the input (emitted ultrasonic pulse) can be used. The metrics and equations that model some nonlinear situations are carefully reviewed. Also, closed form equations of the third-order metrics from a simplified second-order Volterra kernel are derived. Computer simulations show that the surrogate data technique is a potentially powerful tool for blind detection of nonlinearities in multiple-echo ultrasonic signals if adequate metrics are chosen. They also reveal interesting trade-offs among parameters that model ultrasonic systems and detection percentages.