Hydraulic fracturing is the primary method of stimulation in unconventional reservoirs, playing a significant role in oil and gas production enhancement. A key issue for the analysis of hydraulic fracture initiation is to accurately determine the stress distributions in the vicinity of the borehole caused by the injection of pressurized fluids. This paper develops an exact, three-dimensional, poroelastic coupled analytical solution for such stress analysis of an arbitrarily inclined borehole subjected concurrently to a finite-length fluid discharge and
in situ
stresses, using Fourier expansion theorem and the Laplace–Fourier integral transform technique. The complicated boundary conditions, which involve the mixed boundary values at the borehole surface and the coupling between the total radial stress and injection-induced pore pressure over the sectioned borehole interval, as well as the fully three-dimensional far field
in situ
stresses, are addressed in a novel way and deliberately/elegantly decomposed into five fundamental, easier to handle modes. The rigour and definitive nature of the proposed analytical methodology facilitates fundamental understanding of the mechanism underlying the stress responses of the borehole and porous medium. It can be and is used here as a benchmark for the numerical solutions obtained from the finite-element analysis commercial program (ABAQUS).