For the first time, we report the effect of interference between different optical channels on the two-photon absorption (TPA) process in three dimensions. We have employed response theory as well as a sum-over-states (SOS) approach involving few intermediate states to calculate the TPA parameters like transition probabilities (δ(TP)) and TPA tensor elements. In order to use the limited SOS approach, we have derived a new formula for a generalized few-state-model (GFSM) in three dimensions. Due to the presence of additional terms related to the angle between different transition moment vectors, the channel interference associated with the TPA process in 3D is significantly different and much more complicated than that in 1D and 2D cases. The entire study has been carried out on the two simplest Reichardt's dyes, namely 2- and 4-(pyridinium-1-yl)-phenolate (ortho- and para-betain) in gas phase, THF, CH(3)CN and water solvents. We have meticulously inspected the effect of the additional angle related terms on the overall TPA transition probabilities of the two 3D isomeric molecules studied and found that the interfering terms involved in the δ(TP) expression contribute both constructively and destructively as well to the overall δ(TP) value. Moreover, the interfering term has a more conspicuous role in determining the net δ(TP) associated with charge transfer transition in comparison to that of π-π* transition of the studied systems. Interestingly, our model calculations suggest that, for o- and p-betain, the quenching of destructive interference associated with a particular two-photon process can be done with high polarity solvents while the enhancement of constructive interference will be achieved in solvents having relatively small polarity. All the one- and two-photon parameters are evaluated using a range separated CAMB3LYP functional.