The key aspects in the hydrodynamic modelling of wind-generated water waves including irregularity and non-linearity of the wave motion are presented in the report. On the first stage of the investigations the intensity and frequency structure of irregular surface wave field ware determined by using generalized 6th parametric double-peak Hassellman frequency spectrum. Irregular wave surface ware calculated in the frame of classical spectral theory on the so-called quasi-stationary temporal intervals, but irregular wave pressure field and fluid particle velocity field have been derived as a spectral solution of the specially formulated quasi linear boundary value problem for surface waves. On the second stage local group characteristic of wave surface realizations such as amplitude envelope, phase perturbations and corresponding local values of the wave numbers and frequencies ware derived by using Hilbert transform procedure. All these values are considered as input data for the non-linear generalization of the initial quasi linear wave model on the next stages. In the procedure of the non-linear generalization of the model first of all irregular wave surface must be “saturated” by high order bounded harmonics (up to 20th order) with special phase shift to obtain asymmetric surface profiles to be specific for wind generated waves. On the next stage fluid particle velocities on the wave surface ware calculated by using specially developed half-inverse solution of the fully non-linear boundary value problem. For the non-linear calculation of wave particle velocity and pressure fields in the fluid domain Cauchy integral formula (in 2D model) has been used. Testing and systematic numerical calculations and animated pictures displayed realistic results for irregular wave surface and hydrodynamic fields in the fluid domain.