We discuss the variational formulation of the Symmetric Autoencoder (SymAE) and its role in achieving disentanglement within the latent space to extract coherent information from a collection of seismic waveforms. Disentanglement involves separating the latent space into components for coherent information shared by all waveforms and components for waveform-specific nuisance information. SymAE employs a generative model that independently generates waveforms based on coherent and nuisance components, and an inference model that estimates these components from observed wavefield. By assuming the independence of waveforms conditioned on coherent information, the model effectively accumulates this information across multiple waveforms. After training, a metric based on Kullback-Leibler divergence is used to evaluate the informativeness of individual waveforms, enabling latent-space optimization and the generation of synthetic seismograms with enhanced signal-to-noise ratios. To demonstrate the efficacy of our proposed method, we applied it to a data set of teleseismic displacement waveforms of the P wave from deep-focus earthquakes. By training the SymAE model on high-magnitude events, we successfully identified seismograms that contained robust source information. Furthermore, we generated high-resolution virtual seismograms enriched with relevant coherent source information and less influenced by scattering noise, allowing a deeper understanding of the characteristics of the earthquake source. Importantly, our method extracts coherent source information without relying on deconvolution, which is often used in traditional source imaging. This enables the analysis of complex earthquakes with multiple rupture episodes, a capability that is not easily achievable with conventional approaches.