For complex geological reservoir modeling, some numerical-simulation-based methods, such as traditional multiple-point statistics (MPS), cannot extract nonstationary patterns of training images (TIs) effectively. CPU-intensive calculations require that variables information only be stored in RAM instead of other storage mediums to avoid huge time consumption if many simulations are performed successively. Generative Adversarial Network (GAN) modeling methods are able to overcome these issues, but when the training datasets are limited, mode collapse easily occurs. Meanwhile, fixed receptive fields might cause the problem that the corresponding reconstruction cannot take local features and global features of the TIs into account simultaneously. Therefore, we propose the Concurrent Multi-Stage U-Net GAN (Con-MSUGAN) based on a single TI. The pyramid structure is used for spatial pattern multi-scale representation in Con-MSUGAN, which can capture the corresponding feature information of TIs under different scales. Meanwhile, concurrent training patterns can make network parameters in adjacent stages correlated with each other to accelerate network convergence speed, thereby realizing stochastic multi-scale reconstructions. Stationary channel and nonstationary delta facies TIs were used to verify the performance of Con-MSUGAN. Results show that simulation of different datasets using this method are more similar in terms of spatial variance, connectivity degree and facies type frequency distribution to the original TIs. Meanwhile, muti-dimensional-scaling (MDS) plots show small discrepancies between simulation results and the corresponding TI. It indicates that Con-MSUGAN can overcome the problem that complex spatial patterns are hard to reproduce. In addition, saved network parameters can be reused. A variety of equiprobability simulation results can be obtained rapidly, thereby realizing stochastic reconstruction of geological reservoir models.