Data reconstruction is a crucial step in seismic data preprocessing.To improve reconstruction speed and save memory,the commonly used three-dimensional(3D)seismic data reconstruction method divides the missing data in...Data reconstruction is a crucial step in seismic data preprocessing.To improve reconstruction speed and save memory,the commonly used three-dimensional(3D)seismic data reconstruction method divides the missing data into a series of time slices and independently reconstructs each time slice.However,when this strategy is employed,the potential correlations between two adjacent time slices are ignored,which degrades reconstruction performance.Therefore,this study proposes the use of a two-dimensional curvelet transform and the fast iterative shrinkage thresholding algorithm for data reconstruction.Based on the significant overlapping characteristics between the curvelet coefficient support sets of two adjacent time slices,a weighted operator is constructed in the curvelet domain using the prior support set provided by the previous reconstructed time slice to delineate the main energy distribution range,eff ectively providing prior information for reconstructing adjacent slices.Consequently,the resulting weighted fast iterative shrinkage thresholding algorithm can be used to reconstruct 3D seismic data.The processing of synthetic and field data shows that the proposed method has higher reconstruction accuracy and faster computational speed than the conventional fast iterative shrinkage thresholding algorithm for handling missing 3D seismic data.展开更多
Seismic time-frequency(TF)transforms are essential tools in reservoir interpretation and signal processing,particularly for characterizing frequency variations in non-stationary seismic data.Recently,sparse TF trans-f...Seismic time-frequency(TF)transforms are essential tools in reservoir interpretation and signal processing,particularly for characterizing frequency variations in non-stationary seismic data.Recently,sparse TF trans-forms,which leverage sparse coding(SC),have gained significant attention in the geosciences due to their ability to achieve high TF resolution.However,the iterative approaches typically employed in sparse TF transforms are computationally intensive,making them impractical for real seismic data analysis.To address this issue,we propose an interpretable convolutional sparse coding(CSC)network to achieve high TF resolution.The proposed model is generated based on the traditional short-time Fourier transform(STFT)transform and a modified UNet,named ULISTANet.In this design,we replace the conventional convolutional layers of the UNet with learnable iterative shrinkage thresholding algorithm(LISTA)blocks,a specialized form of CSC.The LISTA block,which evolves from the traditional iterative shrinkage thresholding algorithm(ISTA),is optimized for extracting sparse features more effectively.Furthermore,we create a synthetic dataset featuring complex frequency-modulated signals to train ULISTANet.Finally,the proposed method’s performance is subsequently validated using both synthetic and field data,demonstrating its potential for enhanced seismic data analysis.展开更多
基金National Natural Science Foundation of China under Grant 42304145Jiangxi Provincial Natural Science Foundation under Grant 20242BAB26051,20242BAB25191 and 20232BAB213077+1 种基金Foundation of National Key Laboratory of Uranium Resources Exploration-Mining and Nuclear Remote Sensing under Grant 2024QZ-TD-13Open Fund(FW0399-0002)of SINOPEC Key Laboratory of Geophysics。
文摘Data reconstruction is a crucial step in seismic data preprocessing.To improve reconstruction speed and save memory,the commonly used three-dimensional(3D)seismic data reconstruction method divides the missing data into a series of time slices and independently reconstructs each time slice.However,when this strategy is employed,the potential correlations between two adjacent time slices are ignored,which degrades reconstruction performance.Therefore,this study proposes the use of a two-dimensional curvelet transform and the fast iterative shrinkage thresholding algorithm for data reconstruction.Based on the significant overlapping characteristics between the curvelet coefficient support sets of two adjacent time slices,a weighted operator is constructed in the curvelet domain using the prior support set provided by the previous reconstructed time slice to delineate the main energy distribution range,eff ectively providing prior information for reconstructing adjacent slices.Consequently,the resulting weighted fast iterative shrinkage thresholding algorithm can be used to reconstruct 3D seismic data.The processing of synthetic and field data shows that the proposed method has higher reconstruction accuracy and faster computational speed than the conventional fast iterative shrinkage thresholding algorithm for handling missing 3D seismic data.
基金supported by the National Natural Science Foundation of China under Grant 42474139the Key Research and Development Program of Shaanxi under Grant 2024GX-YBXM-067.
文摘Seismic time-frequency(TF)transforms are essential tools in reservoir interpretation and signal processing,particularly for characterizing frequency variations in non-stationary seismic data.Recently,sparse TF trans-forms,which leverage sparse coding(SC),have gained significant attention in the geosciences due to their ability to achieve high TF resolution.However,the iterative approaches typically employed in sparse TF transforms are computationally intensive,making them impractical for real seismic data analysis.To address this issue,we propose an interpretable convolutional sparse coding(CSC)network to achieve high TF resolution.The proposed model is generated based on the traditional short-time Fourier transform(STFT)transform and a modified UNet,named ULISTANet.In this design,we replace the conventional convolutional layers of the UNet with learnable iterative shrinkage thresholding algorithm(LISTA)blocks,a specialized form of CSC.The LISTA block,which evolves from the traditional iterative shrinkage thresholding algorithm(ISTA),is optimized for extracting sparse features more effectively.Furthermore,we create a synthetic dataset featuring complex frequency-modulated signals to train ULISTANet.Finally,the proposed method’s performance is subsequently validated using both synthetic and field data,demonstrating its potential for enhanced seismic data analysis.