Low-Frequency Expansion Approach for Seismic Data Based on Compressed Sensing in Low SNR

Featured Application: The proposed method can improve the anti-noise-interference ability in the inversion procedure of compressed sensing, improve the low-frequency expansion effect for seismic data with a low SNR and further enhance the accuracy of the subsequent inversion and interpretation work. Abstract: Low-frequency information can reﬂect the basic trend of a formation, enhance the accuracy of velocity analysis and improve the imaging accuracy of deep structures in seismic exploration. However, the low-frequency information obtained by the conventional seismic acquisition method is seriously polluted by noise, which will be further lost in processing. Compressed sensing (CS) theory is used to exploit the sparsity of the reﬂection coefﬁcient in the frequency domain to expand the low-frequency components reasonably, thus improving the data quality. However, the conventional CS method is greatly affected by noise, and the effective expansion of low-frequency information can only be realized in the case of a high signal-to-noise ratio (SNR). In this paper, well information is introduced into the objective function to constrain the inversion process of the estimated reﬂection coefﬁcient, and then, the low-frequency component of the original data is expanded by extracting the low-frequency information of the reﬂection coefﬁcient. It has been proved by model tests and actual data processing results that the objective function of estimating the reﬂection coefﬁcient constrained by well logging data based on CS theory can improve the anti-noise interference ability of the inversion process and expand the low-frequency information well in the case of a low SNR.


Introduction
The propagation distance of the low-frequency components in seismic data is much longer than that of the high-frequency components, and the content it carries is more abundant. Due to the slow attenuation of low-frequency information and strong penetrating ability, it can be used to improve the imaging quality of complex structures [1,2]. Seismic exploration on land faces challenges from complex surface conditions, such as severe lateral velocity variation and complex terrain [3]. Deep seismic exploration in this complex situation has higher requirements for low-frequency information. In the full waveform inversion, the frequency information below 10 Hz plays an important role [4]. If the seismic data contain high-fidelity low-frequency information, the accuracy requirements for the initial model will be greatly reduced. Zhang et al. analyzed in detail the impact of missing low-frequency information on the resolution of the profile and waveform inversion [5]. At present, low-frequency information below 10 Hz in seismic data obtained by conventional acquisition methods is severely contaminated by noise, and the effective low-frequency information will be further lost during processing such as surface wave suppression. Therefore, the recovery of low-frequency information has very practical significance for deep seismic exploration.
Deconvolution is an indispensable part of seismic data processing. It compresses the wavelets by eliminating the filtering effect of the earth to obtain the seismic reflection coefficient and ultimately improve the resolution of the data [6,7]. In the 1950s, Silvia and Robinson first proposed a time series convolution model [8], which regarded the seismic trace as a convolution of a white noise reflection sequence and a minimum phase seismic wavelet, and proposed a predictive deconvolution method. Then, Wiggins proposed the minimum entropy deconvolution algorithm by adjusting the inverse wavelet to maximize the output timing variance mode, finally obtaining the reflection coefficient. The algorithm does not need to make any assumptions about the wavelet phase [9]. Claerbout, Gray, Ooe, Godfrey and Walden have successively conducted in-depth research on this method [10][11][12][13][14]. Larue et al. and Bann et al. proposed a deconvolution algorithm based on minimizing the mutual information rate and implemented the deconvolution algorithm in the frequency domain [15,16]. Rosa et al. introduced the spectral simulation algorithm into the deconvolution process, which further improved the resolution of seismic data [17]. In addition, Guan et al. [18] used predictive filtering methods in the frequency domain to extrapolate the frequency spectrum of seismic data and successfully expanded the low-frequency components of the seismic signal. Masoomzadeh et al. [19] and Guo Shuxiang et al. [20], respectively, proposed spatial spectrum whitening and frequency division denoising processing methods, which can effectively protect the low-frequency energy in seismic data. After that, Woodbum et al. [2] used the zero-phase matching operator to perform frequency expansion on the low-frequency part of the zero-phase wavelet to improve the amplitude spectrum of a low frequency signal and realize the expansion for the low-frequency component. Han et al. [21] and Zhang et al. [5] successively proposed a low-frequency expansion method based on the compressive sensing (CS) theory and achieved certain results. However, this method is greatly affected by noise. When the noise interference is weak, the reflection coefficient obtained is more accurate, but when the noise is strong, the accuracy of the obtained reflection coefficient will be greatly reduced.
The collected data from well logs have a very high vertical resolution and are rich in information of different frequency bands. Well logging is a significant tool that assists geophysicists with making constraints, guidance and quantity analysis regarding low resolution seismic data [22]. Zhang et al. proposed an inversion method using logging data to improve the resolution and the signal-to-noise ratio (SNR) of seismic records [23]. Firstly, the initial reflection coefficient model is constructed using logging, seismic and geological information. Then, the reflection coefficient is convoluted with the estimated seismic wavelet to obtain a synthetic seismic record. Finally, the correlation between synthetic records and actual seismic records is analyzed, and the reflection coefficient model is modified repeatedly to maximize the correlation coefficient between synthetic records and actual seismic records, thereby obtaining the optimal reflection coefficient. This method has a relatively large dependence on the accuracy of the initial reflection coefficient model during inversion. If the initial model has a wider frequency band, the final reflection coefficient will be more accurate. Pedersen-Tatalovic et al. proposed a multivariate interpolation method based on logging data to improve the accuracy of lowfrequency prediction by estimating the layer velocity, layer depth, and layer thickness of seismic data [24]. This method can significantly improve the wave impedance information above 2-8 Hz. Due to considerations of the geological background, the modeling process of this method is more complicated.
In order to improve the accuracy of low-frequency compensation and avoid the dependence of the well-constrained seismic trace inversion method on the accuracy of the initial model, this paper develops low-frequency protection and frequency expansion processing technology of well-constrained seismic data under the framework of CS and uses the well data to constrain the sparse inversion process so as to realize low-frequency compensation for noisy seismic data.

CS Theory
CS theory was proposed by David Donoho and Emmanuel Candès [25][26][27][28]. It is a technique that reconstructs sparse, compressible signals from under-determined random linear measurements [29]. Generally, CS theory consists of a sparse representation, a measurement matrix and a sparse promotion algorithm. The two principles of the theory are the sparsity of the signal and the uncorrelation between the measurement matrix and the sparse transformation matrix.
Assuming that transform domain set Ψ can sparsely represent discrete data d∈R N , the matrix form is expressed as follows: where Ψ = [Ψ 1 , . . . , Ψ N ]∈R N is the basis function vector in the transform domain, Ψ T is the inverse transformation of Ψ, N is the total number of basis functions, vector e = [e 1 , e 2 , . . . ,e N ], e j = <d, Ψ j > (j = 1, . . . , N) is the coefficients in the sparse domain.
Measurement signal d obs ∈R n can be regarded as complete signal d obtained by mapping the sampling matrix, Φ∈R n×N : where n is the number of the collected data points, n << N, Φ: = RM,M is the measurement matrix and R is the constraint matrix. Substitute Formula (1) into Formula (2): where Θ = [Θ 1 , . . . , Θ P ] is the model matrix, and Θ i : = RMΨ i T ∈R n×N . According to the theory of CS, the premise of stably reconstructing the K i maxima of the ith feature of an N dimensional signal is that Θ i needs to satisfy the restricted isometry property (RIP) [30,31]. For a K sparse vector (length is N), if Θ∈R M×N satisfies the following conditions: where δ is the RIP constant. Generally, if 0 < δ < 1 is satisfied, it is considered that all submatrices of matrix Θ have at most K columns that satisfy the condition, and Θ satisfies the K-order RIP property with δ. Any matrix can be used as a measurement matrix for data reconstruction as long as it conforms to the RIP properties described in Equation (4) [32,33], which theoretically guarantees that the K sparse signal d can be completely recovered from the measurement signal d obs .
In CS theory, the number of observation signals is much smaller than the length of the signal itself, and the recovery of K sparse signal d by Formula (3) requires solving the underdetermined equations. If the signal d is sparse, the problem of solving underdetermined equations can be transformed into a minimum zero norm problem: where A CS is the model matrix, which consists of the sampling matrix and sparse transformation matrix. Combinatorial optimization is an NP problem. When N is very large, it cannot be effectively realized numerically, and the anti-noise ability is very poor. When measurement matrix Φ satisfies RIP, Formula (5) can be transformed into a convex optimization problem that is numerically easy to deal with the l1 constraint: The method of solving Equation (6) can be roughly divided into three categories: a greedy algorithm, convex optimization algorithm and other algorithms (such as Bayesian class, iterative threshold, minimum total variation and related improved algorithms). The greedy algorithm [34][35][36] selects a local optimal solution in each iteration of the solution process, and finally obtains the global optimal solution. The most important feature of this algorithm is its fast solving speed. The disadvantage of the algorithm is that it requires a great deal of input data, and the accuracy of the optimal solution is low. A convex optimization algorithm [37,38] is a more global optimization to improve the inefficiency in the greedy algorithm. It converts the non-convex optimization problem into a convex optimization problem to calculate the approximate value of the signal and minimizes it under the global criterion. This algorithm has a high calculation accuracy and requires less measurement data. The disadvantage of this algorithm is its slow calculation speed. Iterative threshold technology is often used in sparse optimization algorithms [39][40][41]. For example, when processing high-dimensional data, this technology can reduce the computational complexity and improve computational speed in each iteration.

Low Frequency Expansion
The basic idea of the convolution model is as follows: when a seismic wavelet passes through the underground reflection interface, it will produce a reflection wave, which can be superimposed to obtain a seismic record. Figure 1 shows the synthetic seismogram convoluted by seismic wavelet (30 Hz) and seismic reflection coefficients. It can be seen that the seismic reflection coefficients are sparse in time domain, and its spectrum is full bandwidth, while the low-frequency information of synthetic seismic record is missing, which is caused by the filtering effect of seismic wavelet. Due to the time-domain sparsity and the frequency-domain full bandwidth of seismic reflection coefficients, Han et al. [22] and Zhang et al. [5] successively proposed the problem of sparse inversion with structural L1 norm constraints. In the frequency domain, the full-bandwidth spectrum is reconstructed through a limited bandwidth, thereby recovering the full-bandwidth energy to a certain extent.

Basic Theory
Suppose w is seismic wavelet, r is underground reflection coefficient, g is noisy seismic data, and n is random noise; g can be described as: Formula (7) is transformed to the frequency domain by Fourier transform: where F is the basis function of Fourier transform. G, W, R and N are sparse representations of g, w, r and n in the frequency domain.
The reflection coefficient r in convolution model is sparse in time domain, which conforms to the precondition of CS theory. According to CS theory, the L 1 norm optimization problem for solving the underdetermined linear equation in Equation (8) is as follows: where • 2 2 is the error function of L2 norm constraint; W is the measurement matrix constructed by seismic wavelet w, W ij = w i−τ j +1 , τ j represents the reflection coefficient location; F is the sparse representation matrix; • 1 1 is the exact solution of L 1 norm constraint; λ is a regularization parameter, which is used to adjust the weights of the two norms to ensure the balance between sparsity and data fidelity. The larger the λ is, the larger the weight of the L 1 norm is, and the sparser the reflection coefficient is. The first term of Equation (9) is constrained by L 2 norm and converges to complete seismic data in the process of solution. The second term of Equation (9) is constrained by the L 1 norm to obtain more sparse reflection coefficient r.
Based on the obtained sparse reflection coefficient r, the full band spectrum R can be expressed as: R is the spectrum of the seismic reflection coefficient in the full-frequency band obtained from limited bandwidth seismic data. Assuming that the original data spectrum is G, the frequency spectrum after compensation is G. The specific method of using R to expand the low frequency information in seismic data is as follows: (1) Select the point where R and G are equal; (2) The spectrum value in R before the point and the spectrum value in G after the point are extracted and combined into the spectrum after low-frequency expansion; the formula is described as follows: where L and H are the operators to extract the low-frequency and high-frequency spectrum values, respectively.
(3) The seismic data after low-frequency expansion can be obtained by inverse Fourier transform of G.
The low-frequency expansion process based on CS is shown in Figure 2. Note that in Equation (8), the sampling matrix W formed by the Fourier transform of w is not a completely random function, so the missing frequency information of the seismic data is also recovered to a certain extent.

Well Constrained Low-Frequency Expansion Method
In the framework of CS theory, low-frequency information can be recovered in the frequency domain by using the sparsity of reflection coefficients in the time domain and frequency domain. This method is greatly affected by noise; that is, the reflection coefficient is more accurate when noise is weak. Otherwise, the accuracy of the reflection coefficient will be greatly reduced. In this section, well data are introduced into the cost function to constrain the inversion process so as to improve the accuracy of the reflection coefficient inversion.
For the noisy seismic data G, the cost function with well constraint is constructed as follows: Note that in Equation (8), the sampling matrix W formed by the Fourier transform of w is not a completely random function, so the missing frequency information of the seismic data is also recovered to a certain extent.

Well Constrained Low-Frequency Expansion Method
In the framework of CS theory, low-frequency information can be recovered in the frequency domain by using the sparsity of reflection coefficients in the time domain and frequency domain. This method is greatly affected by noise; that is, the reflection coefficient is more accurate when noise is weak. Otherwise, the accuracy of the reflection coefficient will be greatly reduced. In this section, well data are introduced into the cost function to constrain the inversion process so as to improve the accuracy of the reflection coefficient inversion.
For the noisy seismic data G, the cost function with well constraint is constructed as follows: min where σ is the noise standard deviation estimation of reflection coefficient, r w is the reflection coefficient at the well location, andr w is the reflection coefficient estimated from the well information.
The cost function can also be solved by a greedy algorithm, base pursuit algorithm or iterative threshold algorithm. The iterative threshold method has a simple structure, less calculation and good stability, so this method is adopted to solve the cost function. The iterative formula is as follows: where WF is the measurement matrix, which is composed of wavelet operator and Fourier transform operator. (WF) T is the conjugate transpose of WF. R w is the sampling operator to extract the location information of the corresponding well. R T w is the conjugate transpose of R w .
The reflection coefficient of anti-noise interference can be obtained by iterating Formula (13), and then the expanded seismic data can be obtained according to Equations (10) and (11).

Algorithm Implementation
In the framework of CS theory, low-frequency information can be recovered in the frequency domain by using the sparsity of reflection coefficients in the time domain and frequency domain. This method is greatly affected by noise; that is, the reflection coefficient is more accurate when noise is weak. Otherwise, the accuracy of reflection coefficient will be greatly reduced. In this section, well data are introduced into the cost function to constrain the inversion process.
The implementation process of the algorithm after introducing the well data is shown in Figure 3.
Formula (13) is solved by gradient descent algorithm, and the specific implementation algorithm is as follows: Initializing: ∀r 0 , L in (maximum number of iterations in the inner loop), L out (maximum number of iterations in the outer loop), ε (residual error), k = 1; Well data processing: the reflection coefficientr w is estimated by well data; Iteration start: r 0 = r 0 ; External circulation: judging G − WFr k 2 > ε and k ≤ L out r 0 = r k−1 Internal circulation: Calculate the spectrum R after frequency expansion: R = F r; Calculate the spectrum G after frequency expansion: G = L( R) + H(G); Obtain the seismic data after low-frequency expansion: g = F −1 ( G).
In the above algorithm, T λ is the threshold function; R w is the sampling operator to extract the location information of the corresponding well; α and λ are the algorithm parameters; α must be greater than the maximum eigenvalue of (WF) T WF; λ adjusts the sparsity of the reflection coefficient; L and H are operators for extracting low-frequency and high-frequency spectrum values, respectively. Appl. Sci. 2021, 11, 5028 8 of 23 Formula (13) is solved by gradient descent algorithm, and the specific implementation algorithm is as follows: Initializing: ∀r0, Lin (maximum number of iterations in the inner loop), Lout (maximum number of iterations in the outer loop), ε (residual error), k = 1; Well data processing: the reflection coefficient ˆw r is estimated by well data; Iteration start: r 0 = r0; External circulation: judging Internal circulation: j = 0 to Lin-1, 1  Figure 4a. The corresponding reflection coefficients and spectrum are shown in Figure 4b,c, respectively. Fourier transform is performed on the composite record to obtain the reflection coefficient of the whole frequency band, and the low-frequency part of the reflection coefficient and the high frequency part of the composite record are taken to form a new spectrum as shown in Figure 5b. Figure 5a is the seismic record after low-frequency expansion. Comparing the results before (Figure 4c) and after frequency expansion (Figure 5b), it can be seen that the Appl. Sci. 2021, 11, 5028 9 of 21 low-frequency energy (0-7 Hz) has been significantly improved from −29 dB to −22 dB, especially the energy at 5 Hz is improved obviously. Clearly, the low-frequency expansion method based on the CS theory can extend the low-frequency of the data well.   Different intensity noise is added to the synthetic seismogram to test the anti-noise ability of this method. The noise data with different SNR (14.72, 9.41, 5.10 and 3.14) are shown in Figures 6a-9a. It can be seen that the effective signal is gradually submerged in the noise with the decrease in SNR. For these noisy seismic records, the sparse inversion method based on compressed sensing and the proposed method are used, respectively. Different intensity noise is added to the synthetic seismogram to test the anti-noise ability of this method. The noise data with different SNR (14.72, 9.41, 5.10 and 3.14) are shown in Figures 6a, 7a, 8a and 9a. It can be seen that the effective signal is gradually submerged in the noise with the decrease in SNR. For these noisy seismic records, the sparse inversion method based on compressed sensing and the proposed method are used, respectively.
When the SNR of seismic record is 14.72, the reflection coefficients are calculated by the CS method and the well constrained CS method, and then the spectrum of seismic record after frequency expansion is compared. Figure 6a is the noisy seismic record, Figure 6b is the spectrum of the noisy data. Figure 6c,e shows the reflection coefficients obtained by the two methods, respectively. It can be seen that both methods can accurately reverse the reflection coefficient. Comparing the spectrum before ( Figure 6b) and after frequency expansion (Figure 6d,f), it can be seen that the low-frequency energy (especially in 0-5 Hz) have been significantly improved from −50 dB to −22 dB. It shows that the lowfrequency information has been extended and the frequency band has also been broadened to a certain extent.  The SNR is reduced to 9.41 (Figure 7a), and Figure 7b is its spectrogram. It can be seen that the energy of low-frequency energy in Figure 7b is higher than that in Figure 5b. This is because the added noise contains low-frequency components. The noisy data are processed by the low-frequency expansion method based on CS and the proposed method. Figure 7c,e shows the reflection coefficients obtained by the two methods, respectively. It can be seen that accuracy of the reflection coefficient reversed by CS theory is reduced, but the accuracy of reflection coefficient reversed by the proposed method is not affected. Comparing the spectrum before ( Figure 7b) and after frequency expansion (Figure 7d,f), it can be seen that both methods still can overcome the noise interference and improve the low-frequency information energy.  The SNR continues to reduce to 5.10, Figure 8a,b is the noisy seismic record and its spectrum. It can be seen that the low-frequency information shows increasing interference with the increase in the noise. Two low-frequency expansion methods are utilized to calculate the reflection coefficient and compensate the frequency. Figure 8c,e shows the reflection coefficients obtained by the two methods, respectively. It can be seen that the accuracy of the reflection coefficients inversed by the compressive sensing theory is further reduced, and the accuracy of the reflection coefficients inversed by the proposed method is basically unaffected. Comparing the spectrum before ( Figure 8b) and after frequency expansion (Figure 8d,f), the low-frequency information based on CS can be better extended before 5 Hz, and the effect is poor at around 10 Hz (pointed by the arrow in Figure 8d). The frequency expansion effect of the proposed method is still not affected by the reduction of SNR; the low-frequency information is extended, and the frequency band is broadened to a certain extent.   Figure 9 shows the contrast of the two methods when the SNR is 3.14. Figure 9a,b is the noisy seismic data and their spectrum, respectively. Figure 9c,e is reflection coefficients obtained by the two methods, respectively. It can be seen that the accuracy of the reflection coefficient based on compressed sensing theory is lower, and the accuracy of the reflection coefficient obtained by the proposed method is still unaffected by noise. Comparing the spectrum before and after frequency expansion (Figure 9d,f), the low-frequency information expansion effect based on CS is poor, especially for the low-frequency information between 4 and 10 Hz (indicated by the arrow in Figure 9d), which is not extended. The proposed method is still not affected by the reduction of SNR. Low-frequency information is extended, and the frequency band is broadened to a certain extent. When the SNR of seismic record is 14.72, the reflection coefficients are calculated by the CS method and the well constrained CS method, and then the spectrum of seismic record after frequency expansion is compared. Figure 6a is the noisy seismic record, Figure 6b is the spectrum of the noisy data. Figure 6c,e shows the reflection coefficients obtained by the two methods, respectively. It can be seen that both methods can accurately reverse the reflection coefficient. Comparing the spectrum before ( Figure 6b) and after frequency expansion (Figure 6d,f), it can be seen that the low-frequency energy (especially in 0-5 Hz) have been significantly improved from −50 dB to −22 dB. It shows that the lowfrequency information has been extended and the frequency band has also been broadened to a certain extent.
The SNR is reduced to 9.41 (Figure 7a), and Figure 7b is its spectrogram. It can be seen that the energy of low-frequency energy in Figure 7b is higher than that in Figure 5b. This is because the added noise contains low-frequency components. The noisy data are processed by the low-frequency expansion method based on CS and the proposed method. Figure 7c,e shows the reflection coefficients obtained by the two methods, respectively. It can be seen that accuracy of the reflection coefficient reversed by CS theory is To further illustrate the effectiveness of the proposed approach, the Power Similarity method is calculated to prove the superiority of the proposed method. The Power Similarity is described as follows: where P c is the extended spectrum of the seismic data without noise, P n is the extended spectrum of the seismic data with noise, generally, P c , P n ∈ [0 10 Hz]. The smaller Q is, the greater the similarity is. Figure 10 quantitatively compares the Power Similarity of the two methods under different SNR. It can be seen that the similarity of low-frequency spectrum energy extended by the two methods are gradually decreases with the decrease of SNR.
The similarity decline speed of the proposed method is slower, which indicates that the effect of low-frequency expansion is better.  In summary, when the SNR of seismic data is high, both methods can accurately estimate the reflection coefficient and effectively extend the low-frequency information. With the decrease in the SNR, the accuracy of the reflection coefficient obtained by the inversion method based on CS gradually reduced. However, accuracy of the reflection coefficient obtained by the proposed method is less affected and can estimate the reflection coefficient more accurately. The expansion effect of low-frequency information is better; that is, its anti-interference ability is stronger.

Field Data Example
The field data consists of 100 traces, and the sampling interval is 2 ms. The low frequency expansion method based on CS and the proposed method are used to extend the low-frequency. Figure 11a,b shows the actual single piece of trace data and its spectrum. It can be seen that the energy of the seismic record is relatively weak at about 5 Hz. The low-frequency expansion method based on CS can effectively extend the low-frequency information, and the energy of the signal about 5 Hz is improved to a certain extent. The proposed method can effectively extend the low-frequency information of 0-10 Hz as a whole.
Due to the limited energy of the extended low-frequency information, it is not easy to show in the full frequency section, so the seismic records are displayed by low frequency filtering. Figure 12a shows the filtering results of the original seismic profile at 2-10 Hz. Figure 12b,c shows the low-frequency expansion results processed by the method based on CS and the proposed method. Through comparative analysis, the following points can be seen: (1) The seismic event at position A in the original section is discontinuous, so it is difficult to judge the basic trend of the formation. After the expansion of the two methods, energy of the seismic event at this position is significantly enhanced, and the basic trend of the formation is clear. (2) The seismic event at B in the original section are also discontinuous, and the continuity has been improved to a certain extend after CS low-frequency expansion.
In the proposed method, the inversion process of the reflection coefficient is constrained by the well data, and the continuity of the seismic event after expansion is substantially improved. (3) The seismic event at the position C is not obvious. After expansion processing of CS method, the low-frequency energy is improved, but the continuity is still poor.
After the expansion processing of the proposed method, the continuity of the event is significantly improved, and the effect of the low-frequency expansion is the best. Due to the limited energy of the extended low-frequency information, it is not easy to show in the full frequency section, so the seismic records are displayed by low frequency filtering. Figure 12a shows the filtering results of the original seismic profile at 2-10 Hz. Figure 12b,c shows the low-frequency expansion results processed by the method based on CS and the proposed method. Through comparative analysis, the following points can be seen: (1) The seismic event at position A in the original section is discontinuous, so it is difficult to judge the basic trend of the formation. After the expansion of the two methods, energy of the seismic event at this position is significantly enhanced, and the basic trend of the formation is clear. (2) The seismic event at B in the original section are also discontinuous, and the continuity has been improved to a certain extend after CS low-frequency expansion. In the proposed method, the inversion process of the reflection coefficient is constrained by the well data, and the continuity of the seismic event after expansion is substantially improved. (3) The seismic event at the position C is not obvious. After expansion processing of CS method, the low-frequency energy is improved, but the continuity is still poor. After Figure 11. Comparisons of low-frequency expansion analysis for actual single trace record. (a) Noisy data, (b) frequency spectrum of original data, (c) data processed by CS, (d) frequency spectrum of (c), (e) data processed by the proposed method, (f) frequency spectrum of (e).
To highlight the effectiveness of the proposed approach, the segment B is taken out and shown in Figure 13. As shown in Figure 13a, the continuity of the seismic event between traces 25 and 34 in the profile is poor. After expansion processing of CS method, the energy of the seismic event between trace 25 and 34 is improved, and the continuity of the seismic event is enhanced. The proposed method can further improve the energy of the seismic event between traces 25 and 34, so as to better reflect the trend of the formation. the expansion processing of the proposed method, the continuity of the event is significantly improved, and the effect of the low-frequency expansion is the best. To highlight the effectiveness of the proposed approach, the segment B is taken out and shown in Figure 13. As shown in Figure 13a, the continuity of the seismic event between traces 25 and 34 in the profile is poor. After expansion processing of CS method, the energy of the seismic event between trace 25 and 34 is improved, and the continuity of the seismic event is enhanced. The proposed method can further improve the energy of the seismic event between traces 25 and 34, so as to better reflect the trend of the formation. Similarly, 1-3 Hz information is taken from the original section (as shown in Figure 14), where a, b and c are the 1-3 Hz information in the original seismic section of the work area, the low-frequency expansion result based on CS and the low-frequency expansion result of the proposed method. The three seismic events framed in the figure have poor continuity. After using the CS method, the energy of low-frequency information is enhanced, and the continuity of the event is improved. After using the proposed method, the energy of low-frequency information is further enhanced, and the resolution has better promotion. Similarly, 1-3 Hz information is taken from the original section (as shown in Figure  14), where a, b and c are the 1-3 Hz information in the original seismic section of the work area, the low-frequency expansion result based on CS and the low-frequency expansion result of the proposed method. The three seismic events framed in the figure have poor continuity. After using the CS method, the energy of low-frequency information is enhanced, and the continuity of the event is improved. After using the proposed method, the energy of low-frequency information is further enhanced, and the resolution has better promotion.

Conclusions
Low-frequency information can reflect the basic trend of the formation and improve the accuracy of velocity analysis and deep imaging. Conventional acquisition technology shows that low-frequency information below 10 Hz in the data is seriously polluted by

Conclusions
Low-frequency information can reflect the basic trend of the formation and improve the accuracy of velocity analysis and deep imaging. Conventional acquisition technology shows that low-frequency information below 10 Hz in the data is seriously polluted by noise, and some low-frequency information will be lost during processing. The lowfrequency expansion method based on CS can reasonably compensate the low-frequency components to improve the quality of the data. However, this method is greatly affected by noise, and can effectively extend the low-frequency information only when the SNR is high. In order to improve the anti-interference ability of the method, well information is introduced in the construction of the objective function to constrain the inversion process of estimating the reflection coefficient. Low-frequency information of the original data is extended by extracting the low-frequency information from the reflection coefficient. Model and actual data processing results show that the proposed method, restricting the reflection coefficients inversion process by the logging data, can obtain the accurate reflection coefficients, improve the anti-noise interference ability of the inversion process, and improve the extension effect of the low-frequency information even if the noise energy is large. In this paper, it is ideal to select the point where R and G are equal. In fact, the slope of the equal point needs to be considered, which is also the problem to be solved in the further.