Multi-component ocean-bottom seismic data offer several advantages,including a high signal-to-noise ratio and comprehensive PP and PS wavefield information.However,the strong reflection characteristics of the sea surf...Multi-component ocean-bottom seismic data offer several advantages,including a high signal-to-noise ratio and comprehensive PP and PS wavefield information.However,the strong reflection characteristics of the sea surface generate various types of ghost waves that severely degrade data quality.The source-side ghost waves not only change the effective wave waveforms but also induce the notch effect in frequency spectra and even produce false structures in stacking profiles.Here,the optimization of the updown deconvolution method is addressed to suppress the source-side ghost waves.Based on this method,a complete source-side ghost suppression process is constructed.In traditional methods,the upgoing and downgoing wavefields are derived through vertical wavenumber calibration.However,the wavenumber-domain method is computationally demanding and relies on precise source-receiver geometry for accurate determination of vertical wavenumbers.To address these problems,the upgoing and downgoing wavefields are extracted using the matching method in this study.Then,the up-down deconvolution process is applied to improve the effectiveness of source-side ghost wave suppression.Finally,synthetic and field data case studies,with detailed discussion,are presented to illustrate that the algorithm of our optimized up-down deconvolution is effective.展开更多
Alzheimer’s disease(AD)is the most common origin of sporadic dementia.Rare familial forms have identified a central role for toxicity based on aggregation of peptide fragments generated from amyloid precursor protein...Alzheimer’s disease(AD)is the most common origin of sporadic dementia.Rare familial forms have identified a central role for toxicity based on aggregation of peptide fragments generated from amyloid precursor protein(APP),named amyloid-beta(Aβ),which exists in two common forms,Aβ1-40(Aβ40)and Aβ1-42(Aβ42).The latter is more neurotoxic.A common clinical biomarker measured in blood is the ratio Aβ42/Aβ40.展开更多
Dynamic deconvolution in the time-frequency domain is an effective method to improve seismic resolution.However,the traditional methods usually suffer from the limitation of window width and shape,resulting in insuffi...Dynamic deconvolution in the time-frequency domain is an effective method to improve seismic resolution.However,the traditional methods usually suffer from the limitation of window width and shape,resulting in insufficient time-frequency resolution and phase distortion of seismic signatures,which significantly limits their applicability in complex geological conditions.To overcome this problem,we propose an adaptive local frequency modulation transform(ALFMT),which is directly embedded in the dynamic deconvolution method.ALFMT can dynamically adjust the shape of the analysis window based on the adaptive local frequency information so that the time-frequency energy is effectively focused near the local frequency,thereby improving the focusing and resolution of the time-frequency representation.The synthetic examples and actual seismic application demonstrate that the ALFMT-based dynamic deconvolution can effectively compensate for the seismic amplitude energy,and the compensated data has a higher resolution and clearer reflection characteristics.展开更多
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.展开更多
Optical-resolution photoacoustic microscopy is a novel imaging technique that combines the advantages of optical and ultrasound imaging,enabling high-resolution visualization of biological tissues at the micrometer sc...Optical-resolution photoacoustic microscopy is a novel imaging technique that combines the advantages of optical and ultrasound imaging,enabling high-resolution visualization of biological tissues at the micrometer scale.However,the divergence of the excited Gaussian beam limits the depth-of-field of the system to less than 100μm,which hinders accurate three-dimensional imaging of living tissues and restrictsits applicability in biological research.Therefore,there is an urgent need for an effective method to enhance the depth-of-field without altering the hardware configuration.This paper presents a photoacoustic microscopy depth-of-field extension method and system based on three-dimensional continuity and sparsity deconvolution.This method utilizes a depth-varying point spread function and incorporates continuity and sparsity con-straints into the deconvolution process to mitigate the effect of background noise,enhancing the stability and accuracy of the depth-of-field extension.Experimental results using tungsten wire phantoms suggest that the depth-of-field of system can be extended to 650 pm,which is 7.2 times greater than conventional system,while improving the resolution of the defocused region by an average factor of 3.5.Furthermore,experiments on zebrafish and nude mouse ears with irregular topologies demonstrate that the proposed method successfully overcomes image blurring and the loss of structural information due to limited depth-of-field.All the results suggest that the system with higher lateral resolution and enhanced depth-of-field has significant potential for a wide range of practical biomedical applications.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
Three-dimensional(3 D)gravitational and magnetic exploration is performed using aerial measurement tools,however,this has difficulties with measuring-height design and the construction of a joint-interpretation scheme...Three-dimensional(3 D)gravitational and magnetic exploration is performed using aerial measurement tools,however,this has difficulties with measuring-height design and the construction of a joint-interpretation scheme.At present,the height in such experiments is set according to the measurement scale,and the distribution characteristics of anomalies are not fully considered.Here,we present the idea of using the attenuation characteristics of a singular-value spectrum to evaluate the contributions of various measurement heights and multi-height combinations for inversion to correctly and reasonably design appropriate measuring heights and the number of various measurement heights to be set.The jointgradient Euler-deconvolution method can accurately obtain the distribution of geological bodies from 3 D gravitational and magnetic data at an improved resolution,and experimental tests confirm these findings.Therefore,an actual 3 D aeromagnetic-data-acquisition and inversion test were carried out in the vicinity of the Zhurihe Iron Mine in Inner Mongolia.The fl ight-height diff erence was set to 60 m,and the specifi c distribution of lodes was obtained by the joint-gradient Euler-deconvolution method.This provides a reliable basis for future detailed exploration and proves that the methods presented in this paper have good practicalapplication eff ects and prospects.展开更多
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.展开更多
Airgun arrays are widely used in marine seismic exploration because signatures excited by airgun arrays have high energy and high-peak bubble ratio, whereas the considerable length and width of the array and ghost ref...Airgun arrays are widely used in marine seismic exploration because signatures excited by airgun arrays have high energy and high-peak bubble ratio, whereas the considerable length and width of the array and ghost reflections make the airgun array signature directional. As a result, the relation of the reflection amplitude with the incident and azimuth angles is variable. This means that the directivity of the airgun array results in a nonstationary wavelet and distorts the relation of the amplitude variation with the incident and azimuth angles. To remove the directivity effect, we propose a nonstationary inversion-based directional deconvolution. At fi rst, the signature as a function of take-off angle and azimuth angle is calculated using the spatial configuration of the airgun array and the near-field signatures. Then, based on the velocity model, the time-variant take-off angles are estimated and directional fi lters are designed using the take-off angles. Finally, the directivity-dependent signatures are shaped to the signature right below the airgun array using nonstationary inversion in the directional deconvolution.展开更多
In inertial confinement fusion (ICF), X-ray coded imaging is considered as the most potential means to diagnose the compressed core. The traditional Richardson-Lucy (RL) method has a strong ability to deblur the i...In inertial confinement fusion (ICF), X-ray coded imaging is considered as the most potential means to diagnose the compressed core. The traditional Richardson-Lucy (RL) method has a strong ability to deblur the image where the noise follows the Poisson distribution. However, it always suffers from over-fitting and noise amplification, especially when the signal-to-noise ratio of image is relatively low. In this paper, we propose an improved deconvolution method for X-ray coded imaging. We model the image data as a set of independent Gaussian distributions and derive the iterative solution with a maximum-likelihood scheme. The experimental results on X-ray coded imaging data demonstrate that this method is superior to the RL method in terms of anti-overfitting and noise suppression.展开更多
In order to overcome the defects that the analysis of multi-well typical curves of shale gas reservoirs is rarely applied to engineering,this study proposes a robust production data analysis method based on deconvolut...In order to overcome the defects that the analysis of multi-well typical curves of shale gas reservoirs is rarely applied to engineering,this study proposes a robust production data analysis method based on deconvolution,which is used for multi-well inter-well interference research.In this study,a multi-well conceptual trilinear seepage model for multi-stage fractured horizontal wells was established,and its Laplace solutions under two different outer boundary conditions were obtained.Then,an improved pressure deconvolution algorithm was used to normalize the scattered production data.Furthermore,the typical curve fitting was carried out using the production data and the seepage model solution.Finally,some reservoir parameters and fracturing parameters were interpreted,and the intensity of inter-well interference was compared.The effectiveness of the method was verified by analyzing the production dynamic data of six shale gas wells in Duvernay area.The results showed that the fitting effect of typical curves was greatly improved due to the mutual restriction between deconvolution calculation parameter debugging and seepage model parameter debugging.Besides,by using the morphological characteristics of the log-log typical curves and the time corresponding to the intersection point of the log-log typical curves of two models under different outer boundary conditions,the strength of the interference between wells on the same well platform was well judged.This work can provide a reference for the optimization of well spacing and hydraulic fracturing measures for shale gas wells.展开更多
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.展开更多
In this paper,we consider the use of blind deconvolution for optoacoustic(photoacoustic)imaging and investigate the performance of the method as means for increasing the resolution of the reconstructed image beyond th...In this paper,we consider the use of blind deconvolution for optoacoustic(photoacoustic)imaging and investigate the performance of the method as means for increasing the resolution of the reconstructed image beyond the physical restrictions of the system.The method is demonstrated with optoacoustic measurement obtained from six-day-old mice,imaged in the near-infrared using a broadband hydrophone in a circular scanning configuration.Wefind that estimates of the unknown point spread function,achieved by blind deconvolution,improve the resolution and contrast in the images and show promise for enhancing optoacoustic images.展开更多
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...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.展开更多
This note is concerned with the H-infinity deconvolution filtering problem for linear time-varying discretetime systems described by state space models, The H-infinity deconvolution filter is derived by proposing a ne...This note is concerned with the H-infinity deconvolution filtering problem for linear time-varying discretetime systems described by state space models, The H-infinity deconvolution filter is derived by proposing a new approach in Krein space. With the new approach, it is clearly shown that the central deconvolution filter in an H-infinity setting is the same as the one in an H2 setting associated with one constructed stochastic state-space model. This insight allows us to calculate the complicated H-infinity deconvolution filter in an intuitive and simple way. The deconvolution filter is calculated by performing Riccati equation with the same order as that of the original system.展开更多
基金supported by the National Key Research and Development Program of China(No.2024YFC2813304)supported by the High-Performance Computing Platform of China University of Geosciences,Beijing。
摘要Multi-component ocean-bottom seismic data offer several advantages,including a high signal-to-noise ratio and comprehensive PP and PS wavefield information.However,the strong reflection characteristics of the sea surface generate various types of ghost waves that severely degrade data quality.The source-side ghost waves not only change the effective wave waveforms but also induce the notch effect in frequency spectra and even produce false structures in stacking profiles.Here,the optimization of the updown deconvolution method is addressed to suppress the source-side ghost waves.Based on this method,a complete source-side ghost suppression process is constructed.In traditional methods,the upgoing and downgoing wavefields are derived through vertical wavenumber calibration.However,the wavenumber-domain method is computationally demanding and relies on precise source-receiver geometry for accurate determination of vertical wavenumbers.To address these problems,the upgoing and downgoing wavefields are extracted using the matching method in this study.Then,the up-down deconvolution process is applied to improve the effectiveness of source-side ghost wave suppression.Finally,synthetic and field data case studies,with detailed discussion,are presented to illustrate that the algorithm of our optimized up-down deconvolution is effective.
摘要Alzheimer’s disease(AD)is the most common origin of sporadic dementia.Rare familial forms have identified a central role for toxicity based on aggregation of peptide fragments generated from amyloid precursor protein(APP),named amyloid-beta(Aβ),which exists in two common forms,Aβ1-40(Aβ40)and Aβ1-42(Aβ42).The latter is more neurotoxic.A common clinical biomarker measured in blood is the ratio Aβ42/Aβ40.
基金the sponsorship of the National Science and Technology Major Project of China for New Oil and Gas Exploration and Development(2024ZD1400102)the National Natural Science Foundation of China(U24B2020)+1 种基金the Key Technology for Geophysical Prediction of Ultra-Deep Carbonate Reservoirs(P24240)the Fundamental Research Funds for the Central Universities(23CX10001A).
摘要Dynamic deconvolution in the time-frequency domain is an effective method to improve seismic resolution.However,the traditional methods usually suffer from the limitation of window width and shape,resulting in insufficient time-frequency resolution and phase distortion of seismic signatures,which significantly limits their applicability in complex geological conditions.To overcome this problem,we propose an adaptive local frequency modulation transform(ALFMT),which is directly embedded in the dynamic deconvolution method.ALFMT can dynamically adjust the shape of the analysis window based on the adaptive local frequency information so that the time-frequency energy is effectively focused near the local frequency,thereby improving the focusing and resolution of the time-frequency representation.The synthetic examples and actual seismic application demonstrate that the ALFMT-based dynamic deconvolution can effectively compensate for the seismic amplitude energy,and the compensated data has a higher resolution and clearer reflection characteristics.
基金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 National Key R&D Program of China[Grant No.2022YFC2402400]the National Natural Science Foundation of China[Grant No.62275062]+2 种基金Project of Shandong Innovation and Startup Community of High-end Medical Apparatus and Instruments[Grant Nos.2023-SGTTXM-002 and 2024-SGTTXM-005]the Shandong Province Technology Innovation Guidance Plan(Central Leading Local Science and Technology Development Fund)[Grant No.YDZX2023115]the Taishan Scholar Special Funding Project of Shandong Province,and the Shandong Laboratory of Advanced Biomaterials and Medical Devices in Weihai[Grant No.ZL202402].
摘要Optical-resolution photoacoustic microscopy is a novel imaging technique that combines the advantages of optical and ultrasound imaging,enabling high-resolution visualization of biological tissues at the micrometer scale.However,the divergence of the excited Gaussian beam limits the depth-of-field of the system to less than 100μm,which hinders accurate three-dimensional imaging of living tissues and restrictsits applicability in biological research.Therefore,there is an urgent need for an effective method to enhance the depth-of-field without altering the hardware configuration.This paper presents a photoacoustic microscopy depth-of-field extension method and system based on three-dimensional continuity and sparsity deconvolution.This method utilizes a depth-varying point spread function and incorporates continuity and sparsity con-straints into the deconvolution process to mitigate the effect of background noise,enhancing the stability and accuracy of the depth-of-field extension.Experimental results using tungsten wire phantoms suggest that the depth-of-field of system can be extended to 650 pm,which is 7.2 times greater than conventional system,while improving the resolution of the defocused region by an average factor of 3.5.Furthermore,experiments on zebrafish and nude mouse ears with irregular topologies demonstrate that the proposed method successfully overcomes image blurring and the loss of structural information due to limited depth-of-field.All the results suggest that the system with higher lateral resolution and enhanced depth-of-field has significant potential for a wide range of practical biomedical applications.
基金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 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 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.
基金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.
基金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.
基金supported by the National Key Research and Development Program of China (Nos. 2017YFC0602203,2017YFC0601606,2017YFC0601305 and 2017YFC0602000)National Science and Technology Major Project task (No.2016ZX05027-002-003)+1 种基金National Natural Science Foundation of China (No. 41604098)State Key Program of National Natural Science of China (No. 41430322)。
摘要Three-dimensional(3 D)gravitational and magnetic exploration is performed using aerial measurement tools,however,this has difficulties with measuring-height design and the construction of a joint-interpretation scheme.At present,the height in such experiments is set according to the measurement scale,and the distribution characteristics of anomalies are not fully considered.Here,we present the idea of using the attenuation characteristics of a singular-value spectrum to evaluate the contributions of various measurement heights and multi-height combinations for inversion to correctly and reasonably design appropriate measuring heights and the number of various measurement heights to be set.The jointgradient Euler-deconvolution method can accurately obtain the distribution of geological bodies from 3 D gravitational and magnetic data at an improved resolution,and experimental tests confirm these findings.Therefore,an actual 3 D aeromagnetic-data-acquisition and inversion test were carried out in the vicinity of the Zhurihe Iron Mine in Inner Mongolia.The fl ight-height diff erence was set to 60 m,and the specifi c distribution of lodes was obtained by the joint-gradient Euler-deconvolution method.This provides a reliable basis for future detailed exploration and proves that the methods presented in this paper have good practicalapplication eff ects and prospects.
基金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.
基金the National Natural Science Foundation of China (No.41474109)the China National Petroleum Corporation under grant number 2016A-33.
摘要Airgun arrays are widely used in marine seismic exploration because signatures excited by airgun arrays have high energy and high-peak bubble ratio, whereas the considerable length and width of the array and ghost reflections make the airgun array signature directional. As a result, the relation of the reflection amplitude with the incident and azimuth angles is variable. This means that the directivity of the airgun array results in a nonstationary wavelet and distorts the relation of the amplitude variation with the incident and azimuth angles. To remove the directivity effect, we propose a nonstationary inversion-based directional deconvolution. At fi rst, the signature as a function of take-off angle and azimuth angle is calculated using the spatial configuration of the airgun array and the near-field signatures. Then, based on the velocity model, the time-variant take-off angles are estimated and directional fi lters are designed using the take-off angles. Finally, the directivity-dependent signatures are shaped to the signature right below the airgun array using nonstationary inversion in the directional deconvolution.
基金Project supported by the National High-Tech ICF Committee of China,Foundation of China Academy of Engineering Physics(Grant Nos.2009A0102003 and 2011B0102021)the National Natural Science Foundation of China(Grant No.10905051)
摘要In inertial confinement fusion (ICF), X-ray coded imaging is considered as the most potential means to diagnose the compressed core. The traditional Richardson-Lucy (RL) method has a strong ability to deblur the image where the noise follows the Poisson distribution. However, it always suffers from over-fitting and noise amplification, especially when the signal-to-noise ratio of image is relatively low. In this paper, we propose an improved deconvolution method for X-ray coded imaging. We model the image data as a set of independent Gaussian distributions and derive the iterative solution with a maximum-likelihood scheme. The experimental results on X-ray coded imaging data demonstrate that this method is superior to the RL method in terms of anti-overfitting and noise suppression.
基金financial support from PetroChina Innovation Foundation。
摘要In order to overcome the defects that the analysis of multi-well typical curves of shale gas reservoirs is rarely applied to engineering,this study proposes a robust production data analysis method based on deconvolution,which is used for multi-well inter-well interference research.In this study,a multi-well conceptual trilinear seepage model for multi-stage fractured horizontal wells was established,and its Laplace solutions under two different outer boundary conditions were obtained.Then,an improved pressure deconvolution algorithm was used to normalize the scattered production data.Furthermore,the typical curve fitting was carried out using the production data and the seepage model solution.Finally,some reservoir parameters and fracturing parameters were interpreted,and the intensity of inter-well interference was compared.The effectiveness of the method was verified by analyzing the production dynamic data of six shale gas wells in Duvernay area.The results showed that the fitting effect of typical curves was greatly improved due to the mutual restriction between deconvolution calculation parameter debugging and seepage model parameter debugging.Besides,by using the morphological characteristics of the log-log typical curves and the time corresponding to the intersection point of the log-log typical curves of two models under different outer boundary conditions,the strength of the interference between wells on the same well platform was well judged.This work can provide a reference for the optimization of well spacing and hydraulic fracturing measures for shale gas wells.
基金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.
摘要In this paper,we consider the use of blind deconvolution for optoacoustic(photoacoustic)imaging and investigate the performance of the method as means for increasing the resolution of the reconstructed image beyond the physical restrictions of the system.The method is demonstrated with optoacoustic measurement obtained from six-day-old mice,imaged in the near-infrared using a broadband hydrophone in a circular scanning configuration.Wefind that estimates of the unknown point spread function,achieved by blind deconvolution,improve the resolution and contrast in the images and show promise for enhancing optoacoustic images.
基金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 National Natural Science Foundation of China (No.60574016,60804034)the Natural Science Foundation of Shandong Province (No.Y2007G34)+2 种基金the National Natural Science Foundation for Distinguished Youth Scholars of China (No.60825304)973 Program (No.2009cb320600)the first two authors are also supported by "Taishan Scholarship" Construction Engineering
摘要This note is concerned with the H-infinity deconvolution filtering problem for linear time-varying discretetime systems described by state space models, The H-infinity deconvolution filter is derived by proposing a new approach in Krein space. With the new approach, it is clearly shown that the central deconvolution filter in an H-infinity setting is the same as the one in an H2 setting associated with one constructed stochastic state-space model. This insight allows us to calculate the complicated H-infinity deconvolution filter in an intuitive and simple way. The deconvolution filter is calculated by performing Riccati equation with the same order as that of the original system.