Riverbank filtration (RBF) is widely used in drinking water production all over the world. The horizontal and radial collector wells are among the most important water production facilities in case of the RBF aquifers. In this study several types of modeling methods are used to analyze the hydraulic conditions of horizontal and radial collector wells. The unsteady and steady-state flow in homogeneous, infinite, semi-infinite, and square-shaped bounded aquifers are investigated. The effect of delayed gravity response is modeled. Three-dimensional analytical and numerical finite-difference (FD) solution techniques are applied during the simulations to provide more reliable calculation results. Comparative analyses are carried out to optimize the wellbore drawdown in FD simulations. This important value is controlled by the local head loss due to the radial flow in vertical plane normal to the axes of screens. The multimesh FD simulation software FLOW is used for infinite and semi-infinite aquifers, while the MODFLOW Revised Multi-Node Well (MNW2) package is applied to the square-shaped aquifers. The late time partitioning (LTP) method is introduced to estimate the late time drawdown in infinite aquifers using semianalytical steadystate drawdown equations considering circular recharge boundaries. The latter value is derived by the semianalytical CW software. The use of results of various mostly analytically based solutions confirmed or improved the accuracy of the approximate numerical solutions. Such analytical support is vital for proper evaluation of the operation of horizontal and radial collector wells.