A full probability method is presented to statistically evaluate the variation of geomechanics numerical simulation results subject to multi-random parameters. Taking the interested results as stochastic implicit response variables

firstly the expectations of some special functions of response variables are determined by the high performance multidimensional numerical integration. Then the full probability distribution function of the interested numerical simulation results is obtained to quantitatively analyze the variation effect using two proposed method including numerical probability analysis method based on fast Fourier transform (FFT) and the maximum entropy distribution solver software developed by the authors. The numerical examples show that this method has higher accuracy and the computational effort is reduced.