Traditional deconvolution methods based on single-channel inversion do not consider the spatial structural relation between channels,and hence,they yield high-resolution results with the existing transverse inconsiste...Traditional deconvolution methods based on single-channel inversion do not consider the spatial structural relation between channels,and hence,they yield high-resolution results with the existing transverse inconsistency or discontinuity.Therefore,in this study,the local dip angle was used to obtain the structural information and construct the spatial structurally constraint operator.This operator is then introduced into multichannel deconvolution as a regularization operator to improve the resolution and maintain the transverse continuity of seismic data.Model tests and actual seismic data processing have demonstrated the effectiveness and practicability of this method.展开更多
Deconvolution is a commonly employed technique for enhancing image quality in optical imaging methods.Unfortu-nately,its application in optical coherence tomography(OCT)is often hindered by sensitivity to noise,which ...Deconvolution is a commonly employed technique for enhancing image quality in optical imaging methods.Unfortu-nately,its application in optical coherence tomography(OCT)is often hindered by sensitivity to noise,which leads to ad-ditive ringing artifacts.These artifacts considerably degrade the quality of deconvolved images,thereby limiting its effect-iveness in OCT imaging.In this study,we propose a framework that integrates numerical random phase masks into the deconvolution process,effectively eliminating these artifacts and enhancing image clarity.The optimized joint operation of an iterative Richardson-Lucy deconvolution and numerical synthesis of random phase masks(RPM),termed as De-conv-RPM,enables a 2.5-fold reduction in full width at half-maximum(FWHM).We demonstrate that the Deconv-RPM method significantly enhances image clarity,allowing for the discernment of previously unresolved cellular-level details in nonkeratinized epithelial cells ex vivo and moving blood cells in vivo.展开更多
Metal–organic frameworks(MOFs) are crystalline porous materials with tunable properties, exhibiting great potential in gas adsorption, separation and catalysis.[1,2]It is challenging to visualize MOFs with transmissi...Metal–organic frameworks(MOFs) are crystalline porous materials with tunable properties, exhibiting great potential in gas adsorption, separation and catalysis.[1,2]It is challenging to visualize MOFs with transmission electron microscopy(TEM) due to their inherent instability under electron beam irradiation. Here, we employ cryo-electron microscopy(cryoEM) to capture images of MOF ZIF-8, revealing inverted-space structural information at a resolution of up to about 1.7A and enhancing its critical electron dose to around 20 e^(-)/A^(2). In addition, it is confirmed by electron-beam irradiation experiments that the high voltage could effectively mitigate the radiolysis, and the structure of ZIF-8 is more stable along the [100] direction under electron beam irradiation. Meanwhile, since the high-resolution electron microscope images are modulated by contrast transfer function(CTF) and it is difficult to determine the positions corresponding to the atomic columns directly from the images. We employ image deconvolution to eliminate the impact of CTF and obtain the structural images of ZIF-8. As a result, the heavy atom Zn and the organic imidazole ring within the organic framework can be distinguished from structural images.展开更多
Laser-induced fluorescence(LIF)spectroscopy is employed for plasma diagnosis,necessitating the utilization of deconvolution algorithms to isolate the Doppler effect from the raw spectral signal.However,direct deconvol...Laser-induced fluorescence(LIF)spectroscopy is employed for plasma diagnosis,necessitating the utilization of deconvolution algorithms to isolate the Doppler effect from the raw spectral signal.However,direct deconvolution becomes invalid in the presence of noise as it leads to infinite amplification of high-frequency noise components.To address this issue,we propose a deconvolution algorithm based on the maximum entropy principle.We validate the effectiveness of the proposed algorithm by utilizing simulated LIF spectra at various noise levels(signal-to-noise ratio,SNR=20–80 d B)and measured LIF spectra with Xe as the working fluid.In the typical measured spectrum(SNR=26.23 d B)experiment,compared with the Gaussian filter and the Richardson–Lucy(R-L)algorithm,the proposed algorithm demonstrates an increase in SNR of 1.39 d B and 4.66 d B,respectively,along with a reduction in the root-meansquare error(RMSE)of 35%and 64%,respectively.Additionally,there is a decrease in the spectral angle(SA)of 0.05 and 0.11,respectively.In the high-quality spectrum(SNR=43.96 d B)experiment,the results show that the running time of the proposed algorithm is reduced by about98%compared with the R-L iterative algorithm.Moreover,the maximum entropy algorithm avoids parameter optimization settings and is more suitable for automatic implementation.In conclusion,the proposed algorithm can accurately resolve Doppler spectrum details while effectively suppressing noise,thus highlighting its advantage in LIF spectral deconvolution applications.展开更多
The Yutu-2 rover onboard the Chang’E-4 mission performed the first lunar penetrating radar detection on the farside of the Moon.The high-frequency channel presented us with many unprecedented details of the subsurfac...The Yutu-2 rover onboard the Chang’E-4 mission performed the first lunar penetrating radar detection on the farside of the Moon.The high-frequency channel presented us with many unprecedented details of the subsurface structures within a depth of approximately 50 m.However,it was still difficult to identify finer layers from the cluttered reflections and scattering waves.We applied deconvolution to improve the vertical resolution of the radar profile by extending the limited bandwidth associated with the emissive radar pulse.To overcome the challenges arising from the mixed-phase wavelets and the problematic amplification of noise,we performed predictive deconvolution to remove the minimum-phase components from the Chang’E-4 dataset,followed by a comprehensive phase rotation to rectify phase anomalies in the radar image.Subsequently,we implemented irreversible migration filtering to mitigate the noise and diminutive clutter echoes amplified by deconvolution.The processed data showed evident enhancement of the vertical resolution with a widened bandwidth in the frequency domain and better signal clarity in the time domain,providing us with more undisputed details of subsurface structures near the Chang’E-4 landing site.展开更多
The conventional nonstationary convolutional model assumes that the seismic signal is recorded at normal incidence. Raw shot gathers are far from this assumption because of the effects of offsets. Because of such prob...The conventional nonstationary convolutional model assumes that the seismic signal is recorded at normal incidence. Raw shot gathers are far from this assumption because of the effects of offsets. Because of such problems, we propose a novel prestack nonstationary deconvolution approach. We introduce the radial trace （RT） transform to the nonstationary deconvolution, we estimate the nonstationary deconvolution factor with hyperbolic smoothing based on variable-step sampling （VSS） in the RT domain, and we obtain the high-resolution prestack nonstationary deconvolution data. The RT transform maps the shot record from the offset and traveltime coordinates to those of apparent velocity and traveltime. The ray paths of the traces in the RT better satisfy the assumptions of the convolutional model. The proposed method combines the advantages of stationary deconvolution and inverse Q filtering, without prior information for Q. The nonstationary deconvolution in the RT domain is more suitable than that in the space-time （XT） domain for prestack data because it is the generalized extension of normal incidence. Tests with synthetic and real data demonstrate that the proposed method is more effective in compensating for large-offset and deep data.展开更多
The Robinson convolution model is mainly restricted by three inappropriate assumptions, i.e., statistically white reflectivity, minimum-phase wavelet, and stationarity. Modern reflectivity inversion methods（e.g., spa...The Robinson convolution model is mainly restricted by three inappropriate assumptions, i.e., statistically white reflectivity, minimum-phase wavelet, and stationarity. Modern reflectivity inversion methods（e.g., sparsity-constrained deconvolution） generally attempt to suppress the problems associated with the first two assumptions but often ignore that seismic traces are nonstationary signals, which undermines the basic assumption of unchanging wavelet in reflectivity inversion. Through tests on reflectivity series, we confirm the effects of nonstationarity on reflectivity estimation and the loss of significant information, especially in deep layers. To overcome the problems caused by nonstationarity, we propose a nonstationary convolutional model, and then use the attenuation curve in log spectra to detect and correct the influences of nonstationarity. We use Gabor deconvolution to handle nonstationarity and sparsity-constrained deconvolution to separating reflectivity and wavelet. The combination of the two deconvolution methods effectively handles nonstationarity and greatly reduces the problems associated with the unreasonable assumptions regarding reflectivity and wavelet. Using marine seismic data, we show that correcting nonstationarity helps recover subtle reflectivity information and enhances the characterization of details with respect to the geological record.展开更多
In seismic data processing, blind deconvolution is a key technology. Introduced in this paper is a flow of one kind of blind deconvolution. The optimal precondition conjugate gradients (PCG) in Kyrlov subspace is als...In seismic data processing, blind deconvolution is a key technology. Introduced in this paper is a flow of one kind of blind deconvolution. The optimal precondition conjugate gradients (PCG) in Kyrlov subspace is also used to improve the stability of the algorithm. The computation amount is greatly decreased.展开更多
Sparsity constrained deconvolution can improve the resolution of band-limited seismic data compared to conventional deconvolution. However, such deconvolution methods result in nonunique solutions and suppress weak re...Sparsity constrained deconvolution can improve the resolution of band-limited seismic data compared to conventional deconvolution. However, such deconvolution methods result in nonunique solutions and suppress weak reflections. The Cauchy function, modified Cauchy function, and Huber function are commonly used constraint criteria in sparse deconvolution. We used numerical experiments to analyze the ability of sparsity constrained deconvolution to restore reflectivity sequences and protect weak reflections under different constraint criteria. The experimental results demonstrate that the performance of sparsity constrained deconvolution depends on the agreement between the constraint criteria and the probability distribution of the reflectivity sequences; furthermore, the modified Cauchy- constrained criterion protects the weak reflections better than the other criteria. Based on the model experiments, the probability distribution of the reflectivity sequences of carbonate and clastic formations is statistically analyzed by using well-logging data and then the modified Cauchy-constrained deconvolution is applied to real seismic data much improving the resolution.展开更多
The predictive deconvolution algorithm （PD）, which is based on second-order statistics, assumes that the primaries and the multiples are implicitly orthogonal. However, the seismic data usually do not satisfy this a...The predictive deconvolution algorithm （PD）, which is based on second-order statistics, assumes that the primaries and the multiples are implicitly orthogonal. However, the seismic data usually do not satisfy this assumption in practice. Since the seismic data （primaries and multiples） have a non-Gaussian distribution, in this paper we present an improved predictive deconvolution algorithm （IPD） by maximizing the non-Gaussianity of the recovered primaries. Applications of the IPD method on synthetic and real seismic datasets show that the proposed method obtains promising results.展开更多
Deconvolution is widely used to increase the resolution of seismic data.To compare the resolution ability of conventional spectrum whitening deconvolution to thin layers with that of welldriven deconvolution,a complex...Deconvolution is widely used to increase the resolution of seismic data.To compare the resolution ability of conventional spectrum whitening deconvolution to thin layers with that of welldriven deconvolution,a complex sedimentary geological model was designed,and then the simulated seismic data were processed respectively by each of the two methods.The amplitude spectrum of seismic data was almost white after spectrum whitening,but the wavelet resolution was low.The amplitude spectrum after well-driven deconvolution deviated from white spectrum,but the wavelet resolution was high.Further analysis showed that if an actual reflectivity series could not well satisfy the hypothesis of white spectrum,spectrum whitening deconvolution had a potential risk of wavelet distortion,which might lead to a pitfall in high resolution seismic data interpretation.On the other hand,the wavelet after well-driven deconvolution had higher resolution both in the time and frequency domains.It is favorable for high resolution seismic interpretation and reservoir prediction.展开更多
Seismic deconvolution plays an important role in the seismic characterization of thin-layer structures and seismic resolution enhancement.However,the trace-by-trace processing strategy is applied and ignores the spati...Seismic deconvolution plays an important role in the seismic characterization of thin-layer structures and seismic resolution enhancement.However,the trace-by-trace processing strategy is applied and ignores the spatial connection along seismic traces,which gives the deconvolved result strong ambiguity and poor spatial continuity.To alleviate this issue,we developed a structurally constrained deconvolution algorithm.The proposed method extracts the refl ection structure characterization from the raw seismic data and introduces it to the multichannel deconvolution algorithm as a spatial refl ection regularization.Benefi ting from the introduction of the reflection regularization,the proposed method enhances the stability and spatial continuity of conventional deconvolution methods.Synthetic and field data examples confi rm the correctness and feasibility of the proposed method.展开更多
Conventional predictive deconvolution assumes that the reflection coefficients of the earth conform to an uncorrelated white noise sequence. The Wiener-Hopf (WH) equation is constructed to solve the filter and elimina...Conventional predictive deconvolution assumes that the reflection coefficients of the earth conform to an uncorrelated white noise sequence. The Wiener-Hopf (WH) equation is constructed to solve the filter and eliminate the correlated components of the seismic records, attenuate multiples, and improve seismic resolution. However, in practice, the primary refl ectivity series of fi eld data rarely satisfy the white noise sequence assumption, with the result that the correlated components of the primary reflectivity series are also eliminated by traditional deconvolution. This results in signal distortion. To solve this problem, we have proposed an improved method for deconvolution. First, we estimated the wavelet correlation from seismic records using the spectrum-modeling method. Second, this wavelet autocorrelation was used to construct a new autocorrelation function which contains the correlated components caused by the existence of multiples and avoids the correlated components of the primary reflectivity series. Finally, the new autocorrelation function was brought into the WH equation, and the predictive fi lter operator was calculated for deconvolution. In this paper, we have applied this new method to simulated and field data processing, and we have compared its performance with that of traditional predictive deconvolution. Our results show that the new method can adapt to non-white refl ectivity series without changing the statistical characteristics of the primary reflection coefficient series. Compared with traditional predictive deconvolution, the new method reduces processing noise and improves fidelity, all while maintaining the ability to attenuate multiples and enhance seismic resolution.展开更多
Visual perception of humans penetrating turbid medium is hampered by scattering.Various techniques have been prompted recently to recover optical imaging through turbid materials.Among them,speckle correlation based o...Visual perception of humans penetrating turbid medium is hampered by scattering.Various techniques have been prompted recently to recover optical imaging through turbid materials.Among them,speckle correlation based on deconvolution is one of the most attractive methods taking advantage of high imaging quality,robustness,eas-of-use,and ease-of-integration.By exploiting the point spread function(PSF)of the scattering system,large Field-of-View,extended Depth-of-Field,noninvasiveness and spectral resoluation are now available as successful solutions for high quality and multifunctional image reconstruction.In this paper,we review the progress of imaging through a scattering medium based on deconvolution method,including the principle,the breakthrough of the limitation of the optical memory ffect,the improvement of the deconvolution algorithm and innovative applications.展开更多
To improve the anti-noise performance of the time-domain Bregman iterative algorithm,an adaptive frequency-domain Bregman sparse-spike deconvolution algorithm is proposed.By solving the Bregman algorithm in the freque...To improve the anti-noise performance of the time-domain Bregman iterative algorithm,an adaptive frequency-domain Bregman sparse-spike deconvolution algorithm is proposed.By solving the Bregman algorithm in the frequency domain,the influence of Gaussian as well as outlier noise on the convergence of the algorithm is effectively avoided.In other words,the proposed algorithm avoids data noise effects by implementing the calculations in the frequency domain.Moreover,the computational efficiency is greatly improved compared with the conventional method.Generalized cross validation is introduced in the solving process to optimize the regularization parameter and thus the algorithm is equipped with strong self-adaptation.Different theoretical models are built and solved using the algorithms in both time and frequency domains.Finally,the proposed and the conventional methods are both used to process actual seismic data.The comparison of the results confirms the superiority of the proposed algorithm due to its noise resistance and self-adaptation capability.展开更多
To subtract the slit function from the measured spectrum, a wavelet-based deconvolution method is proposed to obtain a regularized solution of the problem. The method includes reconstructing the signal from the wavele...To subtract the slit function from the measured spectrum, a wavelet-based deconvolution method is proposed to obtain a regularized solution of the problem. The method includes reconstructing the signal from the wavelet modulus maxima. For the purpose of maxima selection, the spatially selective noise filtration technique was used to distinguish modulus maxima produced by signal from the one created by noise. To test the method, sodium spectrum measured at a wide slit was deconvolved. He-Ne spectrum measured at the corresponding slit width was used as slit function. Sodium measured at a narrow slit was used as the reference spectrum. The deconvolutton result shows that this method can enhance the resolution of the degraded spectrum greatly.展开更多
The shallow subsurface defects are difficult to be identified and quantified by ultrasonic time-of-flight diffraction(TOFD)due to the low resolution induced by pulse width and beam spreading.In this paper,Sparse-SAFT ...The shallow subsurface defects are difficult to be identified and quantified by ultrasonic time-of-flight diffraction(TOFD)due to the low resolution induced by pulse width and beam spreading.In this paper,Sparse-SAFT is proposed to improve the time resolution and lateral resolution in TOFD imaging by combining sparse deconvolution and synthetic aperture focusing technique(SAFT).The mathematical model in the frequency domain is established based on the l1 and l2 norm constraints,and the optimization problem is solved for enhancing time resolution.On this basis,SAFT is employed to improve lateral resolution by delay-and-sum beamforming.The simulated and experimental results indicate that the lateral wave and tip-diffracted waves can be decoupled with Sparse-SAFT.The shallow subsurface defects with a height of 3.0 mm at the depth of 3.0 mm were detected quantitatively,and the relative measurement errors of flaw heights and depths were no more than 10.3%.Compared to conventional SAFT,the time resolution and lateral resolution are enhanced by 72.5 and 56%with Sparse-SAFT,respectively.Finally,the proposed method is also suitable for improving resolution to detect the defects beyond dead zone.展开更多
In this paper, the inverse problem of reconstructing reflectivity function of a medium is examined within a blind deconvolution framework. The ultrasound pulse is estimated using higher-order statistics, and Wiener fi...In this paper, the inverse problem of reconstructing reflectivity function of a medium is examined within a blind deconvolution framework. The ultrasound pulse is estimated using higher-order statistics, and Wiener filter is used to obtain the ultrasonic reflectivity function through wavelet-based models. A new approach to the parameter estimation of the inverse filtering step is proposed in the nondestructive evaluation field, which is based on the theory of Fourier-Wavelet regularized deconvolution (ForWaRD). This new approach can be viewed as a solution to the open problem of adaptation of the ForWaRD framework to perform the convolution kernel estimation and deconvolution interdependently. The results indicate stable solutions of the esti- mated pulse and an improvement in the radio-frequency (RF) signal taking into account its signal-to-noise ratio (SNR) and axial resolution. Simulations and experiments showed that the proposed approach can provide robust and optimal estimates of the reflectivity function.展开更多
For the multi-sensor linear discrete time-invariant stochastic systems with correlated measurement noises and unknown noise statistics,an on-line noise statistics estimator is obtained using the correlation method.Sub...For the multi-sensor linear discrete time-invariant stochastic systems with correlated measurement noises and unknown noise statistics,an on-line noise statistics estimator is obtained using the correlation method.Substituting it into the optimal weighted fusion steady-state white noise deconvolution estimator based on the Kalman filtering,a self-tuning weighted measurement fusion white noise deconvolution estimator is presented.By the Dynamic Error System Analysis(DESA) method,it proved that the self-tuning fusion white noise deconvolution estimator converges to the steady-state optimal fusion white noise deconvolution estimator in a realization.Therefore,it has the asymptotically global optimality.A simulation example for the tracking system with 3 sensors and the Bernoulli-Gaussian input white noise shows its effectiveness.展开更多
基金supported by the basic and forward-looking project(No.2023YQX302)。
文摘Traditional deconvolution methods based on single-channel inversion do not consider the spatial structural relation between channels,and hence,they yield high-resolution results with the existing transverse inconsistency or discontinuity.Therefore,in this study,the local dip angle was used to obtain the structural information and construct the spatial structurally constraint operator.This operator is then introduced into multichannel deconvolution as a regularization operator to improve the resolution and maintain the transverse continuity of seismic data.Model tests and actual seismic data processing have demonstrated the effectiveness and practicability of this method.
基金supported by the Guangdong Natural Science Fund General Program (2023A1515011289)Singapore Ministry of Health's National Medical Research Council under its Open Fund Individual Research Grant (MOH-OFIRG19may-0009)+2 种基金Ministry of Education Singapore under its Academic Research Fund Tier 1 (RG35/22)Academic Research Funding Tier 2 (MOE-T2EP30120-0001)China-Singapore International Joint Research Institute (203-A022001).
文摘Deconvolution is a commonly employed technique for enhancing image quality in optical imaging methods.Unfortu-nately,its application in optical coherence tomography(OCT)is often hindered by sensitivity to noise,which leads to ad-ditive ringing artifacts.These artifacts considerably degrade the quality of deconvolved images,thereby limiting its effect-iveness in OCT imaging.In this study,we propose a framework that integrates numerical random phase masks into the deconvolution process,effectively eliminating these artifacts and enhancing image clarity.The optimized joint operation of an iterative Richardson-Lucy deconvolution and numerical synthesis of random phase masks(RPM),termed as De-conv-RPM,enables a 2.5-fold reduction in full width at half-maximum(FWHM).We demonstrate that the Deconv-RPM method significantly enhances image clarity,allowing for the discernment of previously unresolved cellular-level details in nonkeratinized epithelial cells ex vivo and moving blood cells in vivo.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.12074409 and 12374021)。
文摘Metal–organic frameworks(MOFs) are crystalline porous materials with tunable properties, exhibiting great potential in gas adsorption, separation and catalysis.[1,2]It is challenging to visualize MOFs with transmission electron microscopy(TEM) due to their inherent instability under electron beam irradiation. Here, we employ cryo-electron microscopy(cryoEM) to capture images of MOF ZIF-8, revealing inverted-space structural information at a resolution of up to about 1.7A and enhancing its critical electron dose to around 20 e^(-)/A^(2). In addition, it is confirmed by electron-beam irradiation experiments that the high voltage could effectively mitigate the radiolysis, and the structure of ZIF-8 is more stable along the [100] direction under electron beam irradiation. Meanwhile, since the high-resolution electron microscope images are modulated by contrast transfer function(CTF) and it is difficult to determine the positions corresponding to the atomic columns directly from the images. We employ image deconvolution to eliminate the impact of CTF and obtain the structural images of ZIF-8. As a result, the heavy atom Zn and the organic imidazole ring within the organic framework can be distinguished from structural images.
文摘Laser-induced fluorescence(LIF)spectroscopy is employed for plasma diagnosis,necessitating the utilization of deconvolution algorithms to isolate the Doppler effect from the raw spectral signal.However,direct deconvolution becomes invalid in the presence of noise as it leads to infinite amplification of high-frequency noise components.To address this issue,we propose a deconvolution algorithm based on the maximum entropy principle.We validate the effectiveness of the proposed algorithm by utilizing simulated LIF spectra at various noise levels(signal-to-noise ratio,SNR=20–80 d B)and measured LIF spectra with Xe as the working fluid.In the typical measured spectrum(SNR=26.23 d B)experiment,compared with the Gaussian filter and the Richardson–Lucy(R-L)algorithm,the proposed algorithm demonstrates an increase in SNR of 1.39 d B and 4.66 d B,respectively,along with a reduction in the root-meansquare error(RMSE)of 35%and 64%,respectively.Additionally,there is a decrease in the spectral angle(SA)of 0.05 and 0.11,respectively.In the high-quality spectrum(SNR=43.96 d B)experiment,the results show that the running time of the proposed algorithm is reduced by about98%compared with the R-L iterative algorithm.Moreover,the maximum entropy algorithm avoids parameter optimization settings and is more suitable for automatic implementation.In conclusion,the proposed algorithm can accurately resolve Doppler spectrum details while effectively suppressing noise,thus highlighting its advantage in LIF spectral deconvolution applications.
基金supported by the National Natural Science Foundation of China(Grant Nos.42325406 and 42304187)the China Postdoctoral Science Foundation(Grant No.2023M733476)+3 种基金the CAS Project for Young Scientists in Basic Research(Grant No.YSBR082)the National Key R&D Program of China(Grant No.2022YFF0503203)the Key Research Program of the Institute of Geology and GeophysicsChinese Academy of Sciences(Grant Nos.IGGCAS-202101 and IGGCAS-202401).
文摘The Yutu-2 rover onboard the Chang’E-4 mission performed the first lunar penetrating radar detection on the farside of the Moon.The high-frequency channel presented us with many unprecedented details of the subsurface structures within a depth of approximately 50 m.However,it was still difficult to identify finer layers from the cluttered reflections and scattering waves.We applied deconvolution to improve the vertical resolution of the radar profile by extending the limited bandwidth associated with the emissive radar pulse.To overcome the challenges arising from the mixed-phase wavelets and the problematic amplification of noise,we performed predictive deconvolution to remove the minimum-phase components from the Chang’E-4 dataset,followed by a comprehensive phase rotation to rectify phase anomalies in the radar image.Subsequently,we implemented irreversible migration filtering to mitigate the noise and diminutive clutter echoes amplified by deconvolution.The processed data showed evident enhancement of the vertical resolution with a widened bandwidth in the frequency domain and better signal clarity in the time domain,providing us with more undisputed details of subsurface structures near the Chang’E-4 landing site.
基金financially supported by the National Science and Technology Major Project of China(No.2011ZX05023-005-005)the National Natural Science Foundation of China(No.41274137)
文摘The conventional nonstationary convolutional model assumes that the seismic signal is recorded at normal incidence. Raw shot gathers are far from this assumption because of the effects of offsets. Because of such problems, we propose a novel prestack nonstationary deconvolution approach. We introduce the radial trace （RT） transform to the nonstationary deconvolution, we estimate the nonstationary deconvolution factor with hyperbolic smoothing based on variable-step sampling （VSS） in the RT domain, and we obtain the high-resolution prestack nonstationary deconvolution data. The RT transform maps the shot record from the offset and traveltime coordinates to those of apparent velocity and traveltime. The ray paths of the traces in the RT better satisfy the assumptions of the convolutional model. The proposed method combines the advantages of stationary deconvolution and inverse Q filtering, without prior information for Q. The nonstationary deconvolution in the RT domain is more suitable than that in the space-time （XT） domain for prestack data because it is the generalized extension of normal incidence. Tests with synthetic and real data demonstrate that the proposed method is more effective in compensating for large-offset and deep data.
基金funded by the National Basic Research Program of China(973 Program)(Grant No.2011CB201100)Major Program of the National Natural Science Foundation of China(Grant No.2011ZX05004003)
文摘The Robinson convolution model is mainly restricted by three inappropriate assumptions, i.e., statistically white reflectivity, minimum-phase wavelet, and stationarity. Modern reflectivity inversion methods（e.g., sparsity-constrained deconvolution） generally attempt to suppress the problems associated with the first two assumptions but often ignore that seismic traces are nonstationary signals, which undermines the basic assumption of unchanging wavelet in reflectivity inversion. Through tests on reflectivity series, we confirm the effects of nonstationarity on reflectivity estimation and the loss of significant information, especially in deep layers. To overcome the problems caused by nonstationarity, we propose a nonstationary convolutional model, and then use the attenuation curve in log spectra to detect and correct the influences of nonstationarity. We use Gabor deconvolution to handle nonstationarity and sparsity-constrained deconvolution to separating reflectivity and wavelet. The combination of the two deconvolution methods effectively handles nonstationarity and greatly reduces the problems associated with the unreasonable assumptions regarding reflectivity and wavelet. Using marine seismic data, we show that correcting nonstationarity helps recover subtle reflectivity information and enhances the characterization of details with respect to the geological record.
基金With the support of the key project of Knowledge Innovation, CAS(KZCX1-y01, KZCX-SW-18), Fund of the China National Natural Sciences and the Daqing Oilfield with Grant No. 49894190
文摘In seismic data processing, blind deconvolution is a key technology. Introduced in this paper is a flow of one kind of blind deconvolution. The optimal precondition conjugate gradients (PCG) in Kyrlov subspace is also used to improve the stability of the algorithm. The computation amount is greatly decreased.
基金supported by the Major Basic Research Development Program of China (973 Program)(No.2013CB228606)the National Science foundation of China (No.41174117)+1 种基金the National Major Science-Technology Project (No.2011ZX05031-001)Innovation Fund of PetroChina (No.2010D-5006-0301)
文摘Sparsity constrained deconvolution can improve the resolution of band-limited seismic data compared to conventional deconvolution. However, such deconvolution methods result in nonunique solutions and suppress weak reflections. The Cauchy function, modified Cauchy function, and Huber function are commonly used constraint criteria in sparse deconvolution. We used numerical experiments to analyze the ability of sparsity constrained deconvolution to restore reflectivity sequences and protect weak reflections under different constraint criteria. The experimental results demonstrate that the performance of sparsity constrained deconvolution depends on the agreement between the constraint criteria and the probability distribution of the reflectivity sequences; furthermore, the modified Cauchy- constrained criterion protects the weak reflections better than the other criteria. Based on the model experiments, the probability distribution of the reflectivity sequences of carbonate and clastic formations is statistically analyzed by using well-logging data and then the modified Cauchy-constrained deconvolution is applied to real seismic data much improving the resolution.
基金National 863 Foundation of China(No.2006AA09A102-10)National Natural Science Foundation of China(No.40874056)NCET Fund
文摘The predictive deconvolution algorithm （PD）, which is based on second-order statistics, assumes that the primaries and the multiples are implicitly orthogonal. However, the seismic data usually do not satisfy this assumption in practice. Since the seismic data （primaries and multiples） have a non-Gaussian distribution, in this paper we present an improved predictive deconvolution algorithm （IPD） by maximizing the non-Gaussianity of the recovered primaries. Applications of the IPD method on synthetic and real seismic datasets show that the proposed method obtains promising results.
基金supported by National 973 Key Basic Research Development Program(No.2007CB209608)National 863 High Technology Research Development Program(No.2007AA06Z218)
文摘Deconvolution is widely used to increase the resolution of seismic data.To compare the resolution ability of conventional spectrum whitening deconvolution to thin layers with that of welldriven deconvolution,a complex sedimentary geological model was designed,and then the simulated seismic data were processed respectively by each of the two methods.The amplitude spectrum of seismic data was almost white after spectrum whitening,but the wavelet resolution was low.The amplitude spectrum after well-driven deconvolution deviated from white spectrum,but the wavelet resolution was high.Further analysis showed that if an actual reflectivity series could not well satisfy the hypothesis of white spectrum,spectrum whitening deconvolution had a potential risk of wavelet distortion,which might lead to a pitfall in high resolution seismic data interpretation.On the other hand,the wavelet after well-driven deconvolution had higher resolution both in the time and frequency domains.It is favorable for high resolution seismic interpretation and reservoir prediction.
基金National Key R&D Program of China(No.2018YFA0702504)the National Natural Science Foundation of China(Nos.42074141,41874141)the Strategic Cooperation Technology Projects of CNPC and CUP(ZLZX2020-03).
文摘Seismic deconvolution plays an important role in the seismic characterization of thin-layer structures and seismic resolution enhancement.However,the trace-by-trace processing strategy is applied and ignores the spatial connection along seismic traces,which gives the deconvolved result strong ambiguity and poor spatial continuity.To alleviate this issue,we developed a structurally constrained deconvolution algorithm.The proposed method extracts the refl ection structure characterization from the raw seismic data and introduces it to the multichannel deconvolution algorithm as a spatial refl ection regularization.Benefi ting from the introduction of the reflection regularization,the proposed method enhances the stability and spatial continuity of conventional deconvolution methods.Synthetic and field data examples confi rm the correctness and feasibility of the proposed method.
基金supported by Scientific Research Foundation of Shandong University of Science and Technology for Recruited Talents(No.2017RCJJ034)
文摘Conventional predictive deconvolution assumes that the reflection coefficients of the earth conform to an uncorrelated white noise sequence. The Wiener-Hopf (WH) equation is constructed to solve the filter and eliminate the correlated components of the seismic records, attenuate multiples, and improve seismic resolution. However, in practice, the primary refl ectivity series of fi eld data rarely satisfy the white noise sequence assumption, with the result that the correlated components of the primary reflectivity series are also eliminated by traditional deconvolution. This results in signal distortion. To solve this problem, we have proposed an improved method for deconvolution. First, we estimated the wavelet correlation from seismic records using the spectrum-modeling method. Second, this wavelet autocorrelation was used to construct a new autocorrelation function which contains the correlated components caused by the existence of multiples and avoids the correlated components of the primary reflectivity series. Finally, the new autocorrelation function was brought into the WH equation, and the predictive fi lter operator was calculated for deconvolution. In this paper, we have applied this new method to simulated and field data processing, and we have compared its performance with that of traditional predictive deconvolution. Our results show that the new method can adapt to non-white refl ectivity series without changing the statistical characteristics of the primary reflection coefficient series. Compared with traditional predictive deconvolution, the new method reduces processing noise and improves fidelity, all while maintaining the ability to attenuate multiples and enhance seismic resolution.
基金supported by National Natural Science Foundation of China(Nos.61705035,61575223,11534017 and 61475038)the Project of Department of Education of Guangdong Province(No.2018KTSCX241)+1 种基金State Key Laboratory of Optoelectronic Materials and Technologies(Sun Yat-sen University)STU Scienti¯c Research Foundation for Talents.
文摘Visual perception of humans penetrating turbid medium is hampered by scattering.Various techniques have been prompted recently to recover optical imaging through turbid materials.Among them,speckle correlation based on deconvolution is one of the most attractive methods taking advantage of high imaging quality,robustness,eas-of-use,and ease-of-integration.By exploiting the point spread function(PSF)of the scattering system,large Field-of-View,extended Depth-of-Field,noninvasiveness and spectral resoluation are now available as successful solutions for high quality and multifunctional image reconstruction.In this paper,we review the progress of imaging through a scattering medium based on deconvolution method,including the principle,the breakthrough of the limitation of the optical memory ffect,the improvement of the deconvolution algorithm and innovative applications.
基金supported by the National Natural Science Foundation of China(No.NSFC 41204101)Open Projects Fund of the State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation(No.PLN201733)+1 种基金Youth Innovation Promotion Association of the Chinese Academy of Sciences(No.2015051)Open Projects Fund of the Natural Gas and Geology Key Laboratory of Sichuan Province(No.2015trqdz03)
文摘To improve the anti-noise performance of the time-domain Bregman iterative algorithm,an adaptive frequency-domain Bregman sparse-spike deconvolution algorithm is proposed.By solving the Bregman algorithm in the frequency domain,the influence of Gaussian as well as outlier noise on the convergence of the algorithm is effectively avoided.In other words,the proposed algorithm avoids data noise effects by implementing the calculations in the frequency domain.Moreover,the computational efficiency is greatly improved compared with the conventional method.Generalized cross validation is introduced in the solving process to optimize the regularization parameter and thus the algorithm is equipped with strong self-adaptation.Different theoretical models are built and solved using the algorithms in both time and frequency domains.Finally,the proposed and the conventional methods are both used to process actual seismic data.The comparison of the results confirms the superiority of the proposed algorithm due to its noise resistance and self-adaptation capability.
文摘To subtract the slit function from the measured spectrum, a wavelet-based deconvolution method is proposed to obtain a regularized solution of the problem. The method includes reconstructing the signal from the wavelet modulus maxima. For the purpose of maxima selection, the spatially selective noise filtration technique was used to distinguish modulus maxima produced by signal from the one created by noise. To test the method, sodium spectrum measured at a wide slit was deconvolved. He-Ne spectrum measured at the corresponding slit width was used as slit function. Sodium measured at a narrow slit was used as the reference spectrum. The deconvolutton result shows that this method can enhance the resolution of the degraded spectrum greatly.
基金National Key Research and Development Program of China(Grant No.2019YFA0709003)National Natural Science Foundation of China(Grant No.51905079)Liaoning Revitalization Talents Program(Grant No.XLYC1902082).
文摘The shallow subsurface defects are difficult to be identified and quantified by ultrasonic time-of-flight diffraction(TOFD)due to the low resolution induced by pulse width and beam spreading.In this paper,Sparse-SAFT is proposed to improve the time resolution and lateral resolution in TOFD imaging by combining sparse deconvolution and synthetic aperture focusing technique(SAFT).The mathematical model in the frequency domain is established based on the l1 and l2 norm constraints,and the optimization problem is solved for enhancing time resolution.On this basis,SAFT is employed to improve lateral resolution by delay-and-sum beamforming.The simulated and experimental results indicate that the lateral wave and tip-diffracted waves can be decoupled with Sparse-SAFT.The shallow subsurface defects with a height of 3.0 mm at the depth of 3.0 mm were detected quantitatively,and the relative measurement errors of flaw heights and depths were no more than 10.3%.Compared to conventional SAFT,the time resolution and lateral resolution are enhanced by 72.5 and 56%with Sparse-SAFT,respectively.Finally,the proposed method is also suitable for improving resolution to detect the defects beyond dead zone.
基金Project (No. PRC 03-41/2003) supported by the Ministry of Con-struction of Cuba
文摘In this paper, the inverse problem of reconstructing reflectivity function of a medium is examined within a blind deconvolution framework. The ultrasound pulse is estimated using higher-order statistics, and Wiener filter is used to obtain the ultrasonic reflectivity function through wavelet-based models. A new approach to the parameter estimation of the inverse filtering step is proposed in the nondestructive evaluation field, which is based on the theory of Fourier-Wavelet regularized deconvolution (ForWaRD). This new approach can be viewed as a solution to the open problem of adaptation of the ForWaRD framework to perform the convolution kernel estimation and deconvolution interdependently. The results indicate stable solutions of the esti- mated pulse and an improvement in the radio-frequency (RF) signal taking into account its signal-to-noise ratio (SNR) and axial resolution. Simulations and experiments showed that the proposed approach can provide robust and optimal estimates of the reflectivity function.
基金Supported by National Natural Science Foundation of China (No.60874063)Key Laboratory of Electronics Engineering,College of Heilongjiang Province (No.DZZD2010-5),and Science and Automatic Control Key Laboratory of Heilongjiang University
文摘For the multi-sensor linear discrete time-invariant stochastic systems with correlated measurement noises and unknown noise statistics,an on-line noise statistics estimator is obtained using the correlation method.Substituting it into the optimal weighted fusion steady-state white noise deconvolution estimator based on the Kalman filtering,a self-tuning weighted measurement fusion white noise deconvolution estimator is presented.By the Dynamic Error System Analysis(DESA) method,it proved that the self-tuning fusion white noise deconvolution estimator converges to the steady-state optimal fusion white noise deconvolution estimator in a realization.Therefore,it has the asymptotically global optimality.A simulation example for the tracking system with 3 sensors and the Bernoulli-Gaussian input white noise shows its effectiveness.
