This paper presents a new technique to reduce the order of transmission line circuits simultaneously with respect to multiple parameters. The reduction is based on multidimensional congruence transformation. The proposed algorithm provides efficient means to estimate the response of large distributed circuits simultaneously as a function of frequency and other design parameters.