To achieve efficient recovery of subsurface energy resources, a suitable trajectory needs to be identified for the production well. In this study, a new approach is presented for automated identification of optimum well trajectories in heterogeneous oil/gas reservoirs. The optimisation procedures are as follows. First, a productivity potential map is generated based on the site characterisation data of a reservoir (when available). Second, based on the fast-marching method, well paths are generated from a number of entrance positions to a number of exit points at opposite sides of the reservoir. The well trajectory is also locally constrained by a prescribed maximum curvature to ensure that the well trajectory is drillable. Finally, the optimum well trajectory is selected from all the candidate paths based on the calculation of a benefit-to-cost ratio. If required, a straight directional well path, may also be derived through a linear approximation to the optimised non-linear trajectory by least squares analysis. Model performance has been demonstrated in both 2D and 3D. In the 2D example, the benefit-to-cost ratio of the optimised well is much higher than that of a straight well; in the 3D example, laterals of various curvatures are generated. The applicability of the method is tested by exploring different reservoir heterogeneities and curvature constraints. This approach can be applied to determine the entrance/exit positions and the well path for subsurface energy system development, which is useful for field applications.