How to fit a properly nonlinear classification model from conventional well logs to lithofacies is a key problem for machine learning methods.Kernel methods(e.g.,KFD,SVM,MSVM)are effective attempts to solve this issue...How to fit a properly nonlinear classification model from conventional well logs to lithofacies is a key problem for machine learning methods.Kernel methods(e.g.,KFD,SVM,MSVM)are effective attempts to solve this issue due to abilities of handling nonlinear features by kernel functions.Deep mining of log features indicating lithofacies still needs to be improved for kernel methods.Hence,this work employs deep neural networks to enhance the kernel principal component analysis(KPCA)method and proposes a deep kernel method(DKM)for lithofacies identification using well logs.DKM includes a feature extractor and a classifier.The feature extractor consists of a series of KPCA models arranged according to residual network structure.A gradient-free optimization method is introduced to automatically optimize parameters and structure in DKM,which can avoid complex tuning of parameters in models.To test the validation of the proposed DKM for lithofacies identification,an open-sourced dataset with seven con-ventional logs(GR,CAL,AC,DEN,CNL,LLD,and LLS)and lithofacies labels from the Daniudi Gas Field in China is used.There are eight lithofacies,namely clastic rocks(pebbly,coarse,medium,and fine sand-stone,siltstone,mudstone),coal,and carbonate rocks.The comparisons between DKM and three commonly used kernel methods(KFD,SVM,MSVM)show that(1)DKM(85.7%)outperforms SVM(77%),KFD(79.5%),and MSVM(82.8%)in accuracy of lithofacies identification;(2)DKM is about twice faster than the multi-kernel method(MSVM)with good accuracy.The blind well test in Well D13 indicates that compared with the other three methods DKM improves about 24%in accuracy,35%in precision,41%in recall,and 40%in F1 score,respectively.In general,DKM is an effective method for complex lithofacies identification.This work also discussed the optimal structure and classifier for DKM.Experimental re-sults show that(m_(1),m_(2),O)is the optimal model structure and linear svM is the optimal classifier.(m_(1),m_(2),O)means there are m KPCAs,and then m2 residual units.A workflow to determine an optimal classifier in DKM for lithofacies identification is proposed,too.展开更多
In this paper, we propose a hybrid PML (H-PML) combining the normal absorption factor of convolutional PML (C-PML) with tangential absorption factor of Mutiaxial PML (M-PML). The H-PML boundary conditions can be...In this paper, we propose a hybrid PML (H-PML) combining the normal absorption factor of convolutional PML (C-PML) with tangential absorption factor of Mutiaxial PML (M-PML). The H-PML boundary conditions can better suppress the numerical instability in some extreme models, and the computational speed of finite-element method and the dynamic range are greatly increased using this HPML. We use the finite-element method with a hybrid PML to model the acoustic reflection of the interface when wireline and well logging while drilling (LWD), in a formation with a reflector outside the borehole. The simulation results suggests that the PS- and SP- reflected waves arrive at the same time when the inclination between the well and the outer interface is zero, and the difference in arrival times increases with increasing dip angle. When there are fractures outside the well, the reflection signal is clearer in the subsequent reflection waves and may be used to identify the fractured zone. The difference between the dominant wavelength and the model scale shows that LWD reflection logging data are of higher resolution and quality than wireline acoustic reflection logging.展开更多
In order to optimize the wood internal quality detection and evaluation system and improve the comprehensive utilization rate of wood,this paper invented a set of log internal defect detection and visualization system...In order to optimize the wood internal quality detection and evaluation system and improve the comprehensive utilization rate of wood,this paper invented a set of log internal defect detection and visualization system by using the ultrasonic dry coupling agent method.The detection and visualization analysis of internal log defects were realized through log specimen test.The main conclusions show that the accuracy,reliability and practicability of the system for detecting the internal defects of log specimens have been effectively verified.The system can make the edge of the detected image smooth by interpolation algorithm,and the edge detection algorithm can be used to detect and reflect the location of internal defects of logs accurately.The content mentioned above has good application value for meeting the requirement of increasing demand for wood resources and improving the automation level of wood nondestructive testing instruments.展开更多
D-T_(2)two-dimensional nuclear magnetic resonance(2D NMR)logging technology can distinguish pore fluid types intuitively,and it is widely used in oil and gas exploration.Many 2D NMR inversion methods(e.g.,truncated si...D-T_(2)two-dimensional nuclear magnetic resonance(2D NMR)logging technology can distinguish pore fluid types intuitively,and it is widely used in oil and gas exploration.Many 2D NMR inversion methods(e.g.,truncated singular value decomposition(TSVD),Butler-Reds-Dawson(BRD),LM-norm smoothing,and TIST-L1 regularization methods)have been proposed successively,but most are limited to numerical simulations.This study focused on the applicability of different inversion methods for NMR logging data of various acquisition sequences,from which the optimal inversion method was selected based on the comparative analysis.First,the two-dimensional NMR logging principle was studied.Then,these inversion methods were studied in detail,and the precision and computational efficiency of CPMG and diffusion editing(DE)sequences obtained from oil-water and gas-water models were compared,respectively.The inversion results and calculation time of truncated singular value decomposition(TSVD),Butler-Reds-Dawson(BRD),LM-norm smoothing,and TIST-L1 regularization were compared and analyzed through numerical simulations.The inversion method was optimized to process SP mode logging data from the MR Scanner instrument.The results showed that the TIST-regularization and LM-norm smoothing methods were more accurate for the CPMG and DE sequence echo trains of the oil-water and gas-water models.However,the LM-norm smoothing method was less time-consuming,making it more suitable for logging data processing.A case study in well A25 showed that the processing results by the LM-norm smoothing method were consistent with GEOLOG software.This demonstrates that the LM-norm smoothing method is applicable in practical NMR logging processing.展开更多
Tight and unconventional reservoirs have become the focus with the progress of petroleum exploration and development.Micro-fractures in these reservoirs can effectively improve reservoir permeability,and well-develope...Tight and unconventional reservoirs have become the focus with the progress of petroleum exploration and development.Micro-fractures in these reservoirs can effectively improve reservoir permeability,and well-developed micro-fractures can serve to directly improve productivity.Compared with the centered electrical well logging method,the Micro Spherical Focused Logging(MSFL)is more suitable for microfracture identification due to its high resolution and near borehole wall measuring method.In this study,an anisotropic model is used to depict micro-fractured formation.First,a forward model with microfractured formation,borehole,logging instrument and surrounding rock is established.Subsequently,MSFL responses under different micro-fracture porosity,resistivity,dip angle and borehole radius,are calculated based on the finite element method(FEM).Finally,the MSFL responses under different microfracture parameters are analyzed with the response laws clarified.展开更多
In this paper, a new numerical method, the coupling method of spherical harmonic function spectral and streamline diffusion finite element for unsteady Boltzmann equation in the neutron logging field, is discussed. Th...In this paper, a new numerical method, the coupling method of spherical harmonic function spectral and streamline diffusion finite element for unsteady Boltzmann equation in the neutron logging field, is discussed. The convergence and error estimations of this scheme are proved. Its applications in the field of neutron logging show its effectiveness.展开更多
Because of the constraint mode of the inversion objective function in the traditional resistivity-inversion method of electromagnetic-propagation resistivity logging while drilling(EPR-LWD),obvious differences appear ...Because of the constraint mode of the inversion objective function in the traditional resistivity-inversion method of electromagnetic-propagation resistivity logging while drilling(EPR-LWD),obvious differences appear in the radial and vertical investigation characteristics between the amplitude-ratio and phase difference,which affect the practical application of EPR-LWD data.In this paper,according to the EPR-LWD data,a self-adaptive constraint resistivity-inversion method,which adopts a self-adaptive constraint weighted expression in the objective function to balance the contributions of the phase difference and amplitude attenuation,is proposed.A particle swarm optimization algorithm is also introduced to eliminate the dependence of the accuracy and convergence on the initial value of the inversion.According to the inversion results of multiple classical formation models for EPR-LWD,the differences between the adaptive constraint inversion-resistivity logs with the traditional amplitude-ratio and the phase difference of the resistivity logs are discussed in detail.The results demonstrate that the adaptive resistivity logs take into account the advantages of the amplitude-ratio logs in the radial investigation and phase difference logs in the vertical resolution.Further,it is superior in thin-layer identification and invasion-effect appraisal compared with the single-amplitude-ratio and phase difference logs.The inversion results can provide a theoretical reference for research on the resistivity-inversion method of electromagnetic wave LWD.展开更多
日志主要记录软硬件的运行信息,通过查看系统日志,可以找到系统出现的问题及原因,确保系统的稳定性和正常运行。日志解析的目的是将半结构化的原始日志解析为可阅读的日志模板,现有解析方法往往只注重于对原始日志的解析,而忽略了后期...日志主要记录软硬件的运行信息,通过查看系统日志,可以找到系统出现的问题及原因,确保系统的稳定性和正常运行。日志解析的目的是将半结构化的原始日志解析为可阅读的日志模板,现有解析方法往往只注重于对原始日志的解析,而忽略了后期模板处理,导致结果的精度不能进一步提高。自此,提出了一种日志解析方法FMLogs(logs parsing based on frequency and MinHash algorithm)。该方法通过设计正则表达式和调节阈值参数以获得最佳性能,同时采用了字符级频率统计和MinHash方法对长度相同和不同的日志模板进行合并。FMLogs在七个真实数据集上进行了广泛的实验,取得了0.924的平均解析准确率和0.983的F 1-Score。实验结果表明,FMLogs是一种有效的日志解析方法,在解析日志的同时具有较高的准确性和效率,并能保证性能的稳定。展开更多
This paper introduces briefly the tasks and characteristics of China Continent Science Drilling (CCSD) Well Logging Engineering, the logging methods measured with CCSD, the quality control of original logging informat...This paper introduces briefly the tasks and characteristics of China Continent Science Drilling (CCSD) Well Logging Engineering, the logging methods measured with CCSD, the quality control of original logging information, the logging plan of CCSD, the logging engineering management of CCSD, the logging interpretation and the results and reports made with CCSD.展开更多
Nuclear magnetic resonance logging (NMR) is an open well logging method. Drilling mud resistivity, formation resistivity and sodium ions influence its radio frequency (RF) field strength and NMR logging signals. R...Nuclear magnetic resonance logging (NMR) is an open well logging method. Drilling mud resistivity, formation resistivity and sodium ions influence its radio frequency (RF) field strength and NMR logging signals. Research on these effects can provide an important basis for NMR logging data acquisition and interpretation. Three models, water-based drilling mud--water bearing formation, water- based drilling mud--oil bearing formation, oil-based drilling mud--water bearing formation, were studied by finite element method numerical simulation. The influences of drilling mud resistivity and formation resistivity on the NMR logging tool RF field and the influences of sodium ions on the NMR logging signals were simulated numerically. On the basis of analysis, RF field correction and sodium ion correction formulae were proposed and their application range was also discussed. The results indicate that when drilling mud resistivity and formation resistivity are 0.02 Ω·m and 0.2 Ω·m respectively, the attenuation index of centric NMR logging tool is 8.9% and 9.47% respectively. The RF field of an eccentric NMR logging tool is affected mainly by formation resistivity. When formation resistivity is 0.1 Ω·m, the attenuation index is 17.5%. For centric NMR logging tools, the signals coming from sodium ions can be up to 31.8% of total signal. Suggestions are proposed for further research into NMR logging tool correction method and response characteristics.展开更多
In the past, most of the studies for compressional velocities are based on experimental measurements, which lack the support of field data. The purpose of this study is to estimate the compressional velocities based o...In the past, most of the studies for compressional velocities are based on experimental measurements, which lack the support of field data. The purpose of this study is to estimate the compressional velocities based on well log data of delta front subfacies of Lower Tertiary ages of Ji-Dong oil field, China. At initial stage, we have chosen the well log parameters (effect factors) which strongly influence on compressional velocities and established a new modified equation for compressional velocities, which is based on these effect factors. Then Gardner, De-hua Han and this newly established equation were utilized to calculate the compressional velocities in each well. Finally, Least-square regression was carried out to check the fitting of each equation. Regression results clearly indicate that our purposed equation shows better fitting as compared to Gardner and De-hua Han equations.展开更多
Dispersion and attenuation analysis can be used to determine formation anisotropy induced by fractures,or stresses.In this paper,we propose a nonparametric spectrum estimation method to get phase dispersion characteri...Dispersion and attenuation analysis can be used to determine formation anisotropy induced by fractures,or stresses.In this paper,we propose a nonparametric spectrum estimation method to get phase dispersion characteristics and attenuation coefficient.By designing an appropriate vector filter,phase velocity,attenuation coefficient and amplitude can be inverted from the waveform recorded by the receiver array.Performance analysis of this algorithm is compared with Extended Prony Method(EPM)and Forward and Backward Matrix Pencil(FBMP)method.Based on the analysis results,the proposed method is capable of achieving high resolution and precision as the parametric spectrum estimation methods.At the meantime,it also keeps high stability as the other nonparametric spectrum estimation methods.At last,applications to synthetic waveforms modeled using finite difference method and real data show its efficiency.The real data processing results show that the P-wave attenuation log is more sensitive to oil formation compared to S-wave;and the S-wave attenuation log is more sensitive to shale formation compared to P-wave.展开更多
基金supported by the National Natural Science Foundation of China(Grant No.42002134)China Postdoctoral Science Foundation(Grant No.2021T140735)Science Foundation of China University of Petroleum,Beijing(Grant Nos.2462020XKJS02 and 2462020YXZZ004).
文摘How to fit a properly nonlinear classification model from conventional well logs to lithofacies is a key problem for machine learning methods.Kernel methods(e.g.,KFD,SVM,MSVM)are effective attempts to solve this issue due to abilities of handling nonlinear features by kernel functions.Deep mining of log features indicating lithofacies still needs to be improved for kernel methods.Hence,this work employs deep neural networks to enhance the kernel principal component analysis(KPCA)method and proposes a deep kernel method(DKM)for lithofacies identification using well logs.DKM includes a feature extractor and a classifier.The feature extractor consists of a series of KPCA models arranged according to residual network structure.A gradient-free optimization method is introduced to automatically optimize parameters and structure in DKM,which can avoid complex tuning of parameters in models.To test the validation of the proposed DKM for lithofacies identification,an open-sourced dataset with seven con-ventional logs(GR,CAL,AC,DEN,CNL,LLD,and LLS)and lithofacies labels from the Daniudi Gas Field in China is used.There are eight lithofacies,namely clastic rocks(pebbly,coarse,medium,and fine sand-stone,siltstone,mudstone),coal,and carbonate rocks.The comparisons between DKM and three commonly used kernel methods(KFD,SVM,MSVM)show that(1)DKM(85.7%)outperforms SVM(77%),KFD(79.5%),and MSVM(82.8%)in accuracy of lithofacies identification;(2)DKM is about twice faster than the multi-kernel method(MSVM)with good accuracy.The blind well test in Well D13 indicates that compared with the other three methods DKM improves about 24%in accuracy,35%in precision,41%in recall,and 40%in F1 score,respectively.In general,DKM is an effective method for complex lithofacies identification.This work also discussed the optimal structure and classifier for DKM.Experimental re-sults show that(m_(1),m_(2),O)is the optimal model structure and linear svM is the optimal classifier.(m_(1),m_(2),O)means there are m KPCAs,and then m2 residual units.A workflow to determine an optimal classifier in DKM for lithofacies identification is proposed,too.
基金supported by the National Natural Science Foundation of China(No.41204094)Science Foundation of China University of Petroleum,Beijing(No.2462015YQ0506)
文摘In this paper, we propose a hybrid PML (H-PML) combining the normal absorption factor of convolutional PML (C-PML) with tangential absorption factor of Mutiaxial PML (M-PML). The H-PML boundary conditions can better suppress the numerical instability in some extreme models, and the computational speed of finite-element method and the dynamic range are greatly increased using this HPML. We use the finite-element method with a hybrid PML to model the acoustic reflection of the interface when wireline and well logging while drilling (LWD), in a formation with a reflector outside the borehole. The simulation results suggests that the PS- and SP- reflected waves arrive at the same time when the inclination between the well and the outer interface is zero, and the difference in arrival times increases with increasing dip angle. When there are fractures outside the well, the reflection signal is clearer in the subsequent reflection waves and may be used to identify the fractured zone. The difference between the dominant wavelength and the model scale shows that LWD reflection logging data are of higher resolution and quality than wireline acoustic reflection logging.
文摘In order to optimize the wood internal quality detection and evaluation system and improve the comprehensive utilization rate of wood,this paper invented a set of log internal defect detection and visualization system by using the ultrasonic dry coupling agent method.The detection and visualization analysis of internal log defects were realized through log specimen test.The main conclusions show that the accuracy,reliability and practicability of the system for detecting the internal defects of log specimens have been effectively verified.The system can make the edge of the detected image smooth by interpolation algorithm,and the edge detection algorithm can be used to detect and reflect the location of internal defects of logs accurately.The content mentioned above has good application value for meeting the requirement of increasing demand for wood resources and improving the automation level of wood nondestructive testing instruments.
基金sponsored by the National Natural Science Foundation of China(Nos.42174149,41774144)the National Major Projects(No.2016ZX05014-001).
文摘D-T_(2)two-dimensional nuclear magnetic resonance(2D NMR)logging technology can distinguish pore fluid types intuitively,and it is widely used in oil and gas exploration.Many 2D NMR inversion methods(e.g.,truncated singular value decomposition(TSVD),Butler-Reds-Dawson(BRD),LM-norm smoothing,and TIST-L1 regularization methods)have been proposed successively,but most are limited to numerical simulations.This study focused on the applicability of different inversion methods for NMR logging data of various acquisition sequences,from which the optimal inversion method was selected based on the comparative analysis.First,the two-dimensional NMR logging principle was studied.Then,these inversion methods were studied in detail,and the precision and computational efficiency of CPMG and diffusion editing(DE)sequences obtained from oil-water and gas-water models were compared,respectively.The inversion results and calculation time of truncated singular value decomposition(TSVD),Butler-Reds-Dawson(BRD),LM-norm smoothing,and TIST-L1 regularization were compared and analyzed through numerical simulations.The inversion method was optimized to process SP mode logging data from the MR Scanner instrument.The results showed that the TIST-regularization and LM-norm smoothing methods were more accurate for the CPMG and DE sequence echo trains of the oil-water and gas-water models.However,the LM-norm smoothing method was less time-consuming,making it more suitable for logging data processing.A case study in well A25 showed that the processing results by the LM-norm smoothing method were consistent with GEOLOG software.This demonstrates that the LM-norm smoothing method is applicable in practical NMR logging processing.
基金This research is funded by SINOPEC Group Project P20039-2.
文摘Tight and unconventional reservoirs have become the focus with the progress of petroleum exploration and development.Micro-fractures in these reservoirs can effectively improve reservoir permeability,and well-developed micro-fractures can serve to directly improve productivity.Compared with the centered electrical well logging method,the Micro Spherical Focused Logging(MSFL)is more suitable for microfracture identification due to its high resolution and near borehole wall measuring method.In this study,an anisotropic model is used to depict micro-fractured formation.First,a forward model with microfractured formation,borehole,logging instrument and surrounding rock is established.Subsequently,MSFL responses under different micro-fracture porosity,resistivity,dip angle and borehole radius,are calculated based on the finite element method(FEM).Finally,the MSFL responses under different microfracture parameters are analyzed with the response laws clarified.
文摘In this paper, a new numerical method, the coupling method of spherical harmonic function spectral and streamline diffusion finite element for unsteady Boltzmann equation in the neutron logging field, is discussed. The convergence and error estimations of this scheme are proved. Its applications in the field of neutron logging show its effectiveness.
基金supported by the Foundation of Key Laboratory of Exploration Technology for Oil and Gas Resources of the Ministry of Education, Yangtze University, Wuhan (No. K201812)the Foundation of State Key Laboratory of Petroleum Resources and Prospecting, China University of Petroleum, Beijing (No. PRP/open-1704)the Foundation of Education of Hubei Province, China (No. Q20171304)。
文摘Because of the constraint mode of the inversion objective function in the traditional resistivity-inversion method of electromagnetic-propagation resistivity logging while drilling(EPR-LWD),obvious differences appear in the radial and vertical investigation characteristics between the amplitude-ratio and phase difference,which affect the practical application of EPR-LWD data.In this paper,according to the EPR-LWD data,a self-adaptive constraint resistivity-inversion method,which adopts a self-adaptive constraint weighted expression in the objective function to balance the contributions of the phase difference and amplitude attenuation,is proposed.A particle swarm optimization algorithm is also introduced to eliminate the dependence of the accuracy and convergence on the initial value of the inversion.According to the inversion results of multiple classical formation models for EPR-LWD,the differences between the adaptive constraint inversion-resistivity logs with the traditional amplitude-ratio and the phase difference of the resistivity logs are discussed in detail.The results demonstrate that the adaptive resistivity logs take into account the advantages of the amplitude-ratio logs in the radial investigation and phase difference logs in the vertical resolution.Further,it is superior in thin-layer identification and invasion-effect appraisal compared with the single-amplitude-ratio and phase difference logs.The inversion results can provide a theoretical reference for research on the resistivity-inversion method of electromagnetic wave LWD.
文摘日志主要记录软硬件的运行信息,通过查看系统日志,可以找到系统出现的问题及原因,确保系统的稳定性和正常运行。日志解析的目的是将半结构化的原始日志解析为可阅读的日志模板,现有解析方法往往只注重于对原始日志的解析,而忽略了后期模板处理,导致结果的精度不能进一步提高。自此,提出了一种日志解析方法FMLogs(logs parsing based on frequency and MinHash algorithm)。该方法通过设计正则表达式和调节阈值参数以获得最佳性能,同时采用了字符级频率统计和MinHash方法对长度相同和不同的日志模板进行合并。FMLogs在七个真实数据集上进行了广泛的实验,取得了0.924的平均解析准确率和0.983的F 1-Score。实验结果表明,FMLogs是一种有效的日志解析方法,在解析日志的同时具有较高的准确性和效率,并能保证性能的稳定。
文摘This paper introduces briefly the tasks and characteristics of China Continent Science Drilling (CCSD) Well Logging Engineering, the logging methods measured with CCSD, the quality control of original logging information, the logging plan of CCSD, the logging engineering management of CCSD, the logging interpretation and the results and reports made with CCSD.
基金supported by the National Natural Science Foundation of China (No.41074102)China International Science and Technology Cooperation (No.2009DFA61030)
文摘Nuclear magnetic resonance logging (NMR) is an open well logging method. Drilling mud resistivity, formation resistivity and sodium ions influence its radio frequency (RF) field strength and NMR logging signals. Research on these effects can provide an important basis for NMR logging data acquisition and interpretation. Three models, water-based drilling mud--water bearing formation, water- based drilling mud--oil bearing formation, oil-based drilling mud--water bearing formation, were studied by finite element method numerical simulation. The influences of drilling mud resistivity and formation resistivity on the NMR logging tool RF field and the influences of sodium ions on the NMR logging signals were simulated numerically. On the basis of analysis, RF field correction and sodium ion correction formulae were proposed and their application range was also discussed. The results indicate that when drilling mud resistivity and formation resistivity are 0.02 Ω·m and 0.2 Ω·m respectively, the attenuation index of centric NMR logging tool is 8.9% and 9.47% respectively. The RF field of an eccentric NMR logging tool is affected mainly by formation resistivity. When formation resistivity is 0.1 Ω·m, the attenuation index is 17.5%. For centric NMR logging tools, the signals coming from sodium ions can be up to 31.8% of total signal. Suggestions are proposed for further research into NMR logging tool correction method and response characteristics.
文摘In the past, most of the studies for compressional velocities are based on experimental measurements, which lack the support of field data. The purpose of this study is to estimate the compressional velocities based on well log data of delta front subfacies of Lower Tertiary ages of Ji-Dong oil field, China. At initial stage, we have chosen the well log parameters (effect factors) which strongly influence on compressional velocities and established a new modified equation for compressional velocities, which is based on these effect factors. Then Gardner, De-hua Han and this newly established equation were utilized to calculate the compressional velocities in each well. Finally, Least-square regression was carried out to check the fitting of each equation. Regression results clearly indicate that our purposed equation shows better fitting as compared to Gardner and De-hua Han equations.
基金This research was supported by the National Natural Science Foundation of China(No.42274141)Science Foundation of China University of Petroleum,Beijing(No.2462020YXZZ007).
文摘Dispersion and attenuation analysis can be used to determine formation anisotropy induced by fractures,or stresses.In this paper,we propose a nonparametric spectrum estimation method to get phase dispersion characteristics and attenuation coefficient.By designing an appropriate vector filter,phase velocity,attenuation coefficient and amplitude can be inverted from the waveform recorded by the receiver array.Performance analysis of this algorithm is compared with Extended Prony Method(EPM)and Forward and Backward Matrix Pencil(FBMP)method.Based on the analysis results,the proposed method is capable of achieving high resolution and precision as the parametric spectrum estimation methods.At the meantime,it also keeps high stability as the other nonparametric spectrum estimation methods.At last,applications to synthetic waveforms modeled using finite difference method and real data show its efficiency.The real data processing results show that the P-wave attenuation log is more sensitive to oil formation compared to S-wave;and the S-wave attenuation log is more sensitive to shale formation compared to P-wave.