The aim of the study is to identify specifically temporal and spatial patterns of the intermediate seismic activity in Vrancea seismic zone using new approaches. We have investigated the influence of the principal lunar semidiurnal tidal component M2 on intermediate seismic activity in Vrancea (Romania) sub-crustal region from 1934 to 2009 with a special regard for the time series of events from 1980 to 2009. The constituent is assigned by HiCum stacking method according to the earthquake occurrence. “Schuster” and “Permutation” independent tests are applied to distributions found by stacking. Null hypothesis between seismic activities and selected tidal periodicities is rejected when the statistical p-values obtained by the two tests are less than 5% level of confidence in term of statistics. The stacking function is applied to time series of events belonging to windows shifted in time and space, respectively, to evaluate the variability of correlations in both cases. In the case of 3D shifting domain, a specific algorithm, called “statistical tidal tomography”, is described. The results reveal important issues: a) There is a specific temporal footprint of the p-values around the larger earthquakes; b) A Fast Fourier Transform on the n-order polynomial least squares fit (LSF) of the p values variations emphasizes a long-term period about 17 – 18 years; c) Following the 3-D distribution of p<5% values in different sliding time windows we observe a certain pattern confirmed by the CN algorithm for the earthquake prediction and the future strong Vrancea events monitoring; d) the statistical tidal tomography of M2 component has similar patterns with the analysis of seismicity patterns introduced by others for the Vrancea seismic region.