Electromagnetic telemetry systems through complex geological formations have been increasingly investigated in the last decade due to important engineering applications for the oil and gas industry exploration. Many brute-force computational electromagnetic techniques have been used for modeling this scenario. However, they require a tricky treatment of the large conductivity contrasts present in the soil formations. Also, high-cost computational resources are required for the discretization process and the low-frequency instabilities become critical for such large-scale problems. This research presents a semi-analytic formulation for analyzing the electromagnetic field propagation in a biaxially anisotropic and lossy stratified media. The proposed solution employs a novel spectral domain approach where a Hankel-based integral transform is introduced for modeling wave propagation due to azimuthally symmetric current loop sources. The proposed method is employed for analyzing geophysical scenarios analogous to those of the Brazilian Pre-Salt, where high conductivity carbonate rocks are prevalent. Also, the effect of the pre and post-salt formations on the electromagnetic fields and its interaction with the metallic casing of an oil well is then computed for both isotropic and anisotropic environments. It is presented a series of validation results which show that the proposed technique is numerically stable, robust and computationally efficient for modeling several representative problems of wireless oil well telemetry.