The task of detection statistically significant interaction, its direction and delay between time data series of two oscillatory systems in case of close coupling is investigated with nonlinear modeling approach. Numerical experiments on oscillatory model systems with different coupling function variants are used to study main dependences.