Being capable of serving multiple users with the same radio resource, non-orthogonal multiple access (NOMA) can provide desirable performance enhancements in a fair and spectral efficient manner. In this paper, we investigate the resource allocation (RA) and cluster formation (CF) aspects of NOMA for downlink (DL) uplink (UL) decoupled (DUDe) heterogeneous networks (HetNets). A non-ideal NOMA scheme is considered with power disparity and sensitivity constraints (PDSCs), delay tolerance, and residual interference after cancellation. Taking the PDSCs into account, we analytically show that using the DL decoding order limits UL-NOMA performance by that of OMA, while employing an inverse order result in a performance gain that is mainly determined by the channel gain disparity of users. Thereafter, a generic CF method is proposed for any type of user graph, which iteratively forms clusters using Blossom algorithm. Finally, highly non-convex RA problem is converted into a convex form by employing geometric programming (GP) where power and bandwidth are optimized to maximize network sumrate and max-min fairness objectives.