Structural intensity method is used to formulate vibration power flow model in naturally orthotropic plates in the frequency domain for far-field conditions considering bending waves. Dimensionless parameters are used in classical orthotropic plate equations to get this power. Techniques of Fourier transform and finite difference approximation are used in the formulation. Shear force components of vibration power transmission in x-direction and y-direction are obtained separately. Total power is obtained from the idea of far-field conditions. Cross-spectral density functions of field signal are used to facilitate the estimation of power transmission. Structural intensity is formulated which is similar to that of the conventional two-transducer method. A transducer array of two points is required to get an intensity vector in one direction of a point in the plates. A new bending wave number and a modified Laplace operator are also proposed.