The stress distribution around a circular borehole has been studied extensively.The existing analytical poroelastic solutions,however,often neglect the complex interactions between the solid and fluid that can significantly affect the stress solution.This is important in unconventional gas reservoirs such as coalbeds and shale formations.In order to address this limitation,this paper presents the development of a poroelastic solution that takes into account the effect of gas sorption-induced strain.The solution considers drainage of the reservoir fluid through a vertical wellbore in an isotropic,homogenous,poroelastic rock with non-hydrostatic in situ stress field.The sorption-induced shrinkage of coal is modelled using a Langmuir-type curve which relates the volumetric change to the gas pressure.The redistributed stress field around the wellbore after depletion is found by applying Biot\'s definition of effective stress and Airey\'s stress functions,which leads to a solution of the inhomogeneous Biharmonic equation.Two sets of boundary conditions were considered in order to simulate the unsupported cavity(open-hole)and supported cavity(lined-hole)cases.The implementation was verified against a nu-merical solution for both open-hole and lined-hole cases.A comparative study was then conducted to show the significance of the sorption-induced shrinkage in the stress distribution.Finally,a parametric study analysed the sensitivity of the solution to different poroelastic parameters.The results demonstrate that the developed solution is a useful tool that can be employed alongside complex flow models in order to conduct efficient,field-scale coupled hydro-mechanical simulations,especially when the stress-dependent permeability of the reservoir is of concern.