A three-dimensional adaptive mesh refinement-multiphase lattice Boltzmann flux solver (3D-AMR-MLBFS) is developed for effectively simulating complex multiphase flows with large density ratio and high Reynolds number. In the method, the multiphase lattice Boltzmann flux solver (MLBFS) is used to solve the flow field and the level set method is adopted to capture the interface. In addition, a free energy-based continuum surface tension force (FE-CSF) model, which has a clear physical basis and does not need to introduce a regularization function, is applied to calculate the surface tension force. In order to improve the computational efficiency and at the same time ensure accuracy, the adaptive mesh refinement technology is introduced. This method is used to simulate three-dimensional collision of two equal-sized droplets for a wide range of parametric conditions. We not only verify the developed method by comparing with available literature data but also successfully reproduce four typical outcomes of droplet collision, namely coalescence, reflexive separation, stretching separation, and liquid membrane shattering, for which the evolutions of droplet morphologies and various energies are revealed.