Elastic reverse time migration (ERTM) is developed for better characterization of complex structures by imaging multicomponent seismic data. However, conventional ERTM is subject to limitations such as finite recording aperture, limited bandwidth, and imperfect illumination. Elastic least-squares reverse time migration (ELSRTM) can improve imaging accuracy gradually with iterations by minimizing the residuals between the observed and calculated multicomponent data. Conventional ELSRTM suffers from crosstalk artifacts caused by coupled elastic wavefields with different wave modes. Decomposing the coupled elastic wavefields into pure P- and S-waves is an effective method to suppress these crosstalk artifacts. Considering the trade-off between calculation accuracy and efficiency, we have developed a new ELSRTM scheme for isotropic media based on decoupled wave equations to suppress these wave mode-related crosstalk artifacts in the images of conventional ELSRTM. Pure wavefields are obtained by solving the decoupled wave equations using the finite-difference method in our new ELSRTM method. We also derive new decoupled adjoint-state wave equations that are suitable for the elastic velocity-stress equations in isotropic media. Furthermore, we use the gradient equations based on pure wavefields to update the reflectivity images. Synthetic examples demonstrate that our new ELSRTM method can generate images that better represent the subsurface when compared with conventional ERTM and conventional ELSRTM.

Migration is an important step in seismic data processing for oil and gas exploration. The accuracy of migration directly affects the accuracy of subsequent oil and gas reservoir characterization. Reverse-time migration is one of the most accurate migration methods at present. Multi-wave and multicomponent seismic data contain more P- and S-wave information. Making full use of multi-wave and multicomponent seismic data can offer more information about underground structure and lithology, as well as improve the accuracy of seismic exploration. Elastic reverse-time migration (ERTM) has no dip restriction and can be applied to image multi-wave and multicomponent seismic data in complex structural areas and some special lithology structures. However, the surface topography of complex regions has an influence on wavefield and seriously degrades the quality of ERTM’s migration results. We developed a new ERTM method to migrate multi-wave and multicomponent seismic data in the region with complex surface topography. We first fill the layers between the highest and lowest undulating surface with near-surface elastic parameters in a complex topography model to obtain a new model with a horizontal surface. This allows the finite difference (FD) method based on the regular rectangular grid to be used to numerically solve elastic wave equations in the model with complex topography. The decoupled wave equations are used to generate source P- and S-waves and receiver P- and S-waves to reduce crosstalk artefacts in ERTM. A topography-related filter is further used to remove the influence of surface topography on migration results. The scalar imaging condition is also applied to generate PP and PS migration images. Some numerical examples with different complex topographies demonstrate that our proposed ERTM method can remove the influence of complex topography on ERTM’s images and effectively generate high-quality ERTM images.

Microseismic monitoring is an important technology used to evaluate hydraulic fracturing, and denoising is a crucial processing step. Analyses of the characteristics of acquired three-component microseismic data have indicated that the vertical component has a higher signal-to-noise ratio (SNR) than the two horizontal components. Therefore, we propose a new denoising method for three-component microseismic data using re-constrain variational mode decomposition (VMD). In this method, it is assumed that there is a linear relationship between the modes with the same center frequency among the VMD results of the three-component data. Then, the decomposition result of the vertical component is used as a constraint to the whole denoising effect of the three-component data. On the basis of VMD, we add a constraint condition to form the re-constrain VMD, and deduce the corresponding solution process. According to the synthesis data analysis, the proposed method can not only improve the SNR level of three-component records, it also improves the accuracy of polarization analysis. The proposed method also achieved a satisfactory effect for field data.

Migration-based microseismic event location methods based on signal coherency, such as semblance and cross-correlation, have been extensively used for hydro-fracturing monitoring of unconventional resources because they do not depend on first pick accuracy. However, to a great degree, these methods, which are solely based on either semblance imaging function or cross-correlation imaging function, can be significantly influenced by low signal-to-noise ratio (SNR) circumstances. Thus, in order to enhance the location precision, we have presented a migration-based location method with polarity correction by combining waveform semblance and cross-correlation for downhole microseismic events. In this method, the imaging function is constructed by combining the multiplication of P- and S-wave semblance functions with the sum of the cross-correlation functions between P- and S-waves on each trace. This method can efficiently increase the resolution of the imaging function and improve the localization precision of microseismic events compared with those that utilize semblance and cross-correlation imaging functions. According to synthetic data test results, the location method based on cross-correlation-weighted semblance has lower location uncertainty and shows higher resolution in both directions (approximately parallel and perpendicular to the ray path), even under low SNR or incorrect polarity prediction, as compared with the conventional location approaches based on semblance and cross-correlation. The presented method is less affected by velocity errors than the other two methods. Finally, we also tested the above conventional and proposed methods using field microseismic data obtained from one downhole monitoring well, and verified the feasibility of the presented method.