The future power system, with the unprecedented penetration of renewable energy resources, will be faced with many uncertainties which may lead to problems in grid operation. Consequently, the uncertainty assessment of system performance is essential. Hence, the conventional power flow methods, given the constraints of the electricity and gas networks, may not be suitable for distribution networks such as the multi-carrier microgrids. This paper presents an effective method for optimizing the simultaneous utilization of different energy infrastructures in an environment with different uncertainties considering the constraints of the network. The aim is to study the effect of uncertainties on the optimal operation management of an interconnected microgrids. The fluctuation behavior of loads, renewable resources, and electricity price has been investigated in the proposed model based on probabilistic power flow technique. The results of the system are extracted as random variables, which are depicted in probabilistic and cumulative distribution forms. In this study, the multi-carrier microgrids not only exchange energy with the main grid, but also, energy interchange among the microgrids is possible. Simulations are presented, applying the procedure to an illustrative system of three interconnected MCMGs, and the results justify the effectiveness of the proposed technique.