We investigate the saturation effect, which describes the violation of the linearity between the measured photoacoustic amplitude and the object's optical absorption coefficient in functional photoacoustic imaging when the optical absorption in the object increases. We model the optical energy deposition and photoacoustic signal generation and detection in a semi-infinite optical absorbing obje...