En este estudio, se analizó la propagación de ondas solitarias internas en océanos a gran profundidad. Utilizando análisis multiescala y expansión por perturbación, la ecuación básica se simplifica a la ecuación clásica de Benjamin-Ono con coeficientes variables. Para describir mejor las características de propagación de las ondas solitarias, derivamos una ecuación integral diferencial de orden superior con coeficientes variables (Benjamin-Ono). Posteriormente, se derivó la forma bilineal del modelo utilizando el método bilineal de Hirota y se obtuvo una solución multi-solitón. Basándonos en la solución multi-solitón del modelo, estudiamos más a fondo la interacción del solitón, lo que condujo al descubrimiento de la reflexión de Mach. Se obtuvieron algunas conclusiones que son de valor potencial para futuros estudios de ondas solitarias en el océano.