Abstract: In this paper,NREL 5 MW wind turbines are used as basic model and a large-eddy simulation with an actuator line technique was integrated as a numerical method to simulate a wind farm with multiple wind turbines in nature atmospheric boundary layer.In order to simulate the complex inflow conditions of wind farms,this paper first simulats an atmospheric boundary layer with a volume of 3000 m（length）x3000 m（width） ×1000 m（height）,and verifies the simulation results.The results show that the potential temperature at different heights is unchanged below the capping inversion layer,the average wind speed satisfies the shearing characteristics,and the pulsating wind speed satisfies the turbulence spectrum characteristics;Secondly,the value law of the number of grid nodes（N） on the diameter of the wind rotor and Gaussian distribution factor（ε） in the actuating line model are studied.It is found that when s takes the grid scale（η） as the independent variable,the value of N increase with increase of the coefficient of η.When N is63,the coefficient of η can be 2 or 3.However,η can only be 1.2 when N is 25;Finally,8 wind turbines are arranged in the atmospheric boundary layer using the actuation line model.The wind farm is simulated and the mutual interference between the wind turbines was studied.It was found that the power of the first row of wind turbines is significantly higher than other wind turbines,accounting for 40.3% of the total power output of the wind farm.