Kirkwood-Buff (KB) integrals provide a connection between microscopic properties and thermodynamic properties of multicomponent fluids. The estimation of KB integrals using molecular simulations of finite systems requires accounting for finite size effects. In the small system method, properties of finite subvolumes with different sizes embedded in a larger volume can be used to extrapolate to macroscopic thermodynamic properties. KB integrals computed from small subvolumes scale with the inverse size of the system. This scaling was used to find KB integrals in the thermodynamic limit. To reduce numerical inaccuracies that arise from this extrapolation, alternative approaches were considered in this work. Three methods for computing KB integrals in the thermodynamic limit from information of radial distribution functions (RDFs) of finite systems were compared. These methods allowed for the computation of surface effects. KB integrals and surface terms in the thermodynamic limit were computed for Lennard-Jones (LJ) and Weeks-Chandler-Andersen (WCA) fluids. It was found that all three methods converge to the same value. The main differentiating factor was the speed of convergence with system size . The method that required the smallest size was the one which exploited the scaling of the finite volume KB integral multiplied by . The relationship between KB integrals and surface effects was studied for a range of densities.