In this paper, we built upon the estimating primaries by sparse inversion （EPSI） method. We use the 3D curvelet transform and modify the EPSI method to the sparse inversion of the biconvex optimization and Ll-norm r...In this paper, we built upon the estimating primaries by sparse inversion （EPSI） method. We use the 3D curvelet transform and modify the EPSI method to the sparse inversion of the biconvex optimization and Ll-norm regularization, and use alternating optimization to directly estimate the primary reflection coefficients and source wavelet. The 3D curvelet transform is used as a sparseness constraint when inverting the primary reflection coefficients, which results in avoiding the prediction subtraction process in the surface-related multiples elimination （SRME） method. The proposed method not only reduces the damage to the effective waves but also improves the elimination of multiples. It is also a wave equation- based method for elimination of surface multiple reflections, which effectively removes surface multiples under complex submarine conditions.展开更多
On the assumption that the seismic wavelet amplitude spectrum is estimated accurately, a group of wavelets with different phase spectra, regarded as estimated wavelets, are used to implement linear least-squares inver...On the assumption that the seismic wavelet amplitude spectrum is estimated accurately, a group of wavelets with different phase spectra, regarded as estimated wavelets, are used to implement linear least-squares inversion. During inversion, except for the wavelet phase, all other factors affecting inversion results are not taken into account. The inversion results of a sparse reflectivity model （or blocky impedance model） show that： （1） although the synthetic data using inversion results matches well with the original seismic data, the inverted reflectivity and acoustic impedance are different from that of the real model. （2） the inversion result reliability is dependent on the estimated wavelet Z transform root distribution. When the estimated wavelet Z transform roots only differ from that of the real wavelet near the unit circle, the inverted reflectivity and impedance are usually consistent with the real model; （3） although the synthetic data matches well with the original data and the Cauchy norm （or modified Cauchy norm） with a constant damping parameter has been optimized, the inverted results are still greatly different from the real model. Finally, we suggest using the L1 norm, Kurtosis, variation, Cauchy norm with adaptive damping parameter or/and modified Cauchy norm with adaptive damping parameter as evaluation criteria to reduce the bad influence of inaccurate wavelet phase estimation and obtain good results in theory.展开更多
Monitoring and evaluating the nutritional status of vegetation under stress from exhausted coal mining sites by hyper-spectral remote sensing is important in future ecological restoration engineering. The Wangpingcun ...Monitoring and evaluating the nutritional status of vegetation under stress from exhausted coal mining sites by hyper-spectral remote sensing is important in future ecological restoration engineering. The Wangpingcun coal mine, located in the Mentougou district of Beijing, was chosen as a case study. The ecological damage was analyzed by 3S technology, field investigation and from chemical data. The derivative spectra of the diagnostic absorption bands are derived from the spectra measured in the field and used as characteristic spectral variables. A correlation analysis was conducted for the nitrogen content of the vegetation samples and the fast derivative spectrum and the estimation model of nitrogen content established by a multiple stepwise linear regression method. The spatial distribution of nitrogen content was extracted by a parameter mapping method from the Hyperion data which revealed the distribution of the nitrogen content. In addition, the estimation model was evaluated for two evaluation indicators which are important for the precision of the model. Experimental results indicate that by linear regression and parameter mapping, the estimation model precision was Very high. The coefficient of determination, R2, was 0.795 and the standard deviation of residual （SDR） 0.19. The nitrogen content of most samples was about 1.03% and the nitrogen content in the study site seems inversely proportional to the distance from the piles of coal waste. Therefore, we can conclude that inversely modeling nitrogen content by hyper-spectral remote sensing in exhausted coal mining sites is feasible and our study can be taken as reference in species selection and in subseauent management and maintenance in ecological restoration.展开更多
Cauchy priori distribution-based Bayesian AVO reflectivity inversion may lead to sparse estimates that are sensitive to large reflectivities. For the inversion, the computation of the covariance matrix and regularized...Cauchy priori distribution-based Bayesian AVO reflectivity inversion may lead to sparse estimates that are sensitive to large reflectivities. For the inversion, the computation of the covariance matrix and regularized terms requires prior estimation of model parameters, which makes the iterative inversion weakly nonlinear. At the same time, the relations among the model parameters are assumed linear. Furthermore, the reflectivities, the results of the inversion, or the elastic parameters with cumulative error recovered by integrating reflectivities are not well suited for detecting hydrocarbons and fuids. In contrast, in Bayesian linear AVO inversion, the elastic parameters can be directly extracted from prestack seismic data without linear assumptions for the model parameters. Considering the advantages of the abovementioned methods, the Bayesian AVO reflectivity inversion process is modified and Cauchy distribution is explored as a prior probability distribution and the time-variant covariance is also considered. Finally, we propose a new method for the weakly nonlinear AVO waveform inversion. Furthermore, the linear assumptions are abandoned and elastic parameters, such as P-wave velocity, S-wave velocity, and density, can be directly recovered from seismic data especially for interfaces with large reflectivities. Numerical analysis demonstrates that all the elastic parameters can be estimated from prestack seismic data even when the signal-to-noise ratio of the seismic data is low.展开更多
Generalized Inversion Method has been used to estimate the spatial variation of site effects,using the digital data of SH-waves recorded by 63 stations in the Capital Circle Region of China from 2001 to 2006.We gained...Generalized Inversion Method has been used to estimate the spatial variation of site effects,using the digital data of SH-waves recorded by 63 stations in the Capital Circle Region of China from 2001 to 2006.We gained the site effects of all stations participating in the calculation.We found that the site effect of rock was stabile and about 1.0 from 1.0Hz to 10.0Hz,while the site effect of deposit was high in low frequencies,about 3 ～ 7 from 1.0Hz to 8.0Hz,and the site effect was protuberant at about 5.0Hz,then fell as the frequency increased.The result shows the shape and intensity of station site effects are mainly influenced by the lithology below the station,and possibly also by the local geological structure.展开更多
针对PM2.5遥感模型对气溶胶细粒子比FMF(Fine mode fraction)参数的需求,结合多光谱偏振传感器对大气探测的优势,基于最优估计OE(Optimal Estimation)反演框架,提出了一种基于线偏振度(Degree of Linear Polarization)测量的FMF最优化...针对PM2.5遥感模型对气溶胶细粒子比FMF(Fine mode fraction)参数的需求,结合多光谱偏振传感器对大气探测的优势,基于最优估计OE(Optimal Estimation)反演框架,提出了一种基于线偏振度(Degree of Linear Polarization)测量的FMF最优化反演方法。采用矢量化的辐射传输模式UNL-VRTM进行地基天空光的线偏振度观测模拟,分析了线偏振度对FMF参数的波段敏感性,并基于仿真数据开展了算法的反演测试。研究结果表明:偏振测量在长波近红外波段对FMF的敏感性高于可见光波段;基于OE框架的FMF反演算法具有良好的闭合性;在地基天顶观测模式下,引入线偏振度测量参与反演能够有效提高FMF的反演精度,FMF反演误差从1.4%下降到了0.18%。最优化反演方法对于气溶胶遥感具有一定的潜力和可行性,有望成为提高PM2.5遥感监测能力的新途径。展开更多
The traditional method of Synthetic Aperture Radar(SAR)wind field retrieval is based on an empirical relation between the near surface winds and the normalized radar backscatter cross section to estimate wind speeds,w...The traditional method of Synthetic Aperture Radar(SAR)wind field retrieval is based on an empirical relation between the near surface winds and the normalized radar backscatter cross section to estimate wind speeds,where this relation is called the geophysical model function(GMF).However,the accuracy rapidly decreases due to the impact of rainfall on the measurement of SAR and the saturation of backscattered intensity under the condition of tropical cyclone.Because of no available instrument synchronously monitoring rain rate on the satellite platform of SAR,we have to derive the precipitation of the SAR observation time from non-simultaneous passive microwave observations of rain in combination with geostationary IR images,and then use the model of rain correction to remove the impact of rain on SAR wind field measurements.For the saturation of radar backscatter cross section in high wind speed conditions,we develop an approach to estimate tropical cyclone parameters and wind fields based on the improved Holland model and the SAR image features of tropical cyclone.To retrieve the low-to-moderate wind speed,the wind direction of tropical cyclone is estimated from the SAR image using wavelet analysis.And then the maximum wind speed and the central pressure of tropical cyclone are calculated by a least square minimization of the difference between the improved Holland model and the low-to-moderate wind speed retrieved from SAR.In addition,wind fields are estimated from the improved Holland model using the above-mentioned parameters of tropical cyclone as input.To evaluate the accuracy of our approach,the SAR images of typhoon Aere,typhoon Khanun,and hurricane Ophelia are used to estimate tropical cyclone parameters and wind fields,which are compared with the best track data and reanalyzed wind fields of the Joint Typhoon Warning Center(JTWC)and the Hurricane Research Division(HRD).The results indicate that the tropical cyclone center,maximum wind speed,and central pressure are generally consistent with the best track data,and wind fields agree well with reanalyzed data from HRD.展开更多
Elastic wave inverse scattering theory plays an important role in parameters estimation of heterogeneous media.Combining inverse scattering theory,perturbation theory and stationary phase approximation,we derive the P...Elastic wave inverse scattering theory plays an important role in parameters estimation of heterogeneous media.Combining inverse scattering theory,perturbation theory and stationary phase approximation,we derive the P-wave seismic scattering coefficient equation in terms of fluid factor,shear modulus and density of background homogeneous media and perturbation media.With this equation as forward solver,a pre-stack seismic Bayesian inversion method is proposed to estimate the fluid factor of heterogeneous media.In this method,Cauchy distribution is utilized to the ratios of fluid factors,shear moduli and densities of perturbation media and background homogeneous media,respectively.Gaussian distribution is utilized to the likelihood function.The introduction of constraints from initial smooth models enhances the stability of the estimation of model parameters.Model test and real data example demonstrate that the proposed method is able to estimate the fluid factor of heterogeneous media from pre-stack seismic data directly and reasonably.展开更多
基金supported by the National Science and Technology Major Project (No.2011ZX05023-005-008)
文摘In this paper, we built upon the estimating primaries by sparse inversion （EPSI） method. We use the 3D curvelet transform and modify the EPSI method to the sparse inversion of the biconvex optimization and Ll-norm regularization, and use alternating optimization to directly estimate the primary reflection coefficients and source wavelet. The 3D curvelet transform is used as a sparseness constraint when inverting the primary reflection coefficients, which results in avoiding the prediction subtraction process in the surface-related multiples elimination （SRME） method. The proposed method not only reduces the damage to the effective waves but also improves the elimination of multiples. It is also a wave equation- based method for elimination of surface multiple reflections, which effectively removes surface multiples under complex submarine conditions.
基金supported by National Key Basic Research Development Program (Grant No. 2007CB209600)National Major Science and Technology Program (Grant No. 2008ZX05010-002)
文摘On the assumption that the seismic wavelet amplitude spectrum is estimated accurately, a group of wavelets with different phase spectra, regarded as estimated wavelets, are used to implement linear least-squares inversion. During inversion, except for the wavelet phase, all other factors affecting inversion results are not taken into account. The inversion results of a sparse reflectivity model （or blocky impedance model） show that： （1） although the synthetic data using inversion results matches well with the original seismic data, the inverted reflectivity and acoustic impedance are different from that of the real model. （2） the inversion result reliability is dependent on the estimated wavelet Z transform root distribution. When the estimated wavelet Z transform roots only differ from that of the real wavelet near the unit circle, the inverted reflectivity and impedance are usually consistent with the real model; （3） although the synthetic data matches well with the original data and the Cauchy norm （or modified Cauchy norm） with a constant damping parameter has been optimized, the inverted results are still greatly different from the real model. Finally, we suggest using the L1 norm, Kurtosis, variation, Cauchy norm with adaptive damping parameter or/and modified Cauchy norm with adaptive damping parameter as evaluation criteria to reduce the bad influence of inaccurate wavelet phase estimation and obtain good results in theory.
文摘Monitoring and evaluating the nutritional status of vegetation under stress from exhausted coal mining sites by hyper-spectral remote sensing is important in future ecological restoration engineering. The Wangpingcun coal mine, located in the Mentougou district of Beijing, was chosen as a case study. The ecological damage was analyzed by 3S technology, field investigation and from chemical data. The derivative spectra of the diagnostic absorption bands are derived from the spectra measured in the field and used as characteristic spectral variables. A correlation analysis was conducted for the nitrogen content of the vegetation samples and the fast derivative spectrum and the estimation model of nitrogen content established by a multiple stepwise linear regression method. The spatial distribution of nitrogen content was extracted by a parameter mapping method from the Hyperion data which revealed the distribution of the nitrogen content. In addition, the estimation model was evaluated for two evaluation indicators which are important for the precision of the model. Experimental results indicate that by linear regression and parameter mapping, the estimation model precision was Very high. The coefficient of determination, R2, was 0.795 and the standard deviation of residual （SDR） 0.19. The nitrogen content of most samples was about 1.03% and the nitrogen content in the study site seems inversely proportional to the distance from the piles of coal waste. Therefore, we can conclude that inversely modeling nitrogen content by hyper-spectral remote sensing in exhausted coal mining sites is feasible and our study can be taken as reference in species selection and in subseauent management and maintenance in ecological restoration.
基金supported by the National High-Tech Research and Development Program of China(863 Program)(No.2008AA093001)
文摘Cauchy priori distribution-based Bayesian AVO reflectivity inversion may lead to sparse estimates that are sensitive to large reflectivities. For the inversion, the computation of the covariance matrix and regularized terms requires prior estimation of model parameters, which makes the iterative inversion weakly nonlinear. At the same time, the relations among the model parameters are assumed linear. Furthermore, the reflectivities, the results of the inversion, or the elastic parameters with cumulative error recovered by integrating reflectivities are not well suited for detecting hydrocarbons and fuids. In contrast, in Bayesian linear AVO inversion, the elastic parameters can be directly extracted from prestack seismic data without linear assumptions for the model parameters. Considering the advantages of the abovementioned methods, the Bayesian AVO reflectivity inversion process is modified and Cauchy distribution is explored as a prior probability distribution and the time-variant covariance is also considered. Finally, we propose a new method for the weakly nonlinear AVO waveform inversion. Furthermore, the linear assumptions are abandoned and elastic parameters, such as P-wave velocity, S-wave velocity, and density, can be directly recovered from seismic data especially for interfaces with large reflectivities. Numerical analysis demonstrates that all the elastic parameters can be estimated from prestack seismic data even when the signal-to-noise ratio of the seismic data is low.
基金sponsored by the Special Foundation of China Earthquake Administration (2007-8-26)
文摘Generalized Inversion Method has been used to estimate the spatial variation of site effects,using the digital data of SH-waves recorded by 63 stations in the Capital Circle Region of China from 2001 to 2006.We gained the site effects of all stations participating in the calculation.We found that the site effect of rock was stabile and about 1.0 from 1.0Hz to 10.0Hz,while the site effect of deposit was high in low frequencies,about 3 ～ 7 from 1.0Hz to 8.0Hz,and the site effect was protuberant at about 5.0Hz,then fell as the frequency increased.The result shows the shape and intensity of station site effects are mainly influenced by the lithology below the station,and possibly also by the local geological structure.
文摘针对PM2.5遥感模型对气溶胶细粒子比FMF(Fine mode fraction)参数的需求,结合多光谱偏振传感器对大气探测的优势,基于最优估计OE(Optimal Estimation)反演框架,提出了一种基于线偏振度(Degree of Linear Polarization)测量的FMF最优化反演方法。采用矢量化的辐射传输模式UNL-VRTM进行地基天空光的线偏振度观测模拟,分析了线偏振度对FMF参数的波段敏感性,并基于仿真数据开展了算法的反演测试。研究结果表明:偏振测量在长波近红外波段对FMF的敏感性高于可见光波段;基于OE框架的FMF反演算法具有良好的闭合性;在地基天顶观测模式下,引入线偏振度测量参与反演能够有效提高FMF的反演精度,FMF反演误差从1.4%下降到了0.18%。最优化反演方法对于气溶胶遥感具有一定的潜力和可行性,有望成为提高PM2.5遥感监测能力的新途径。
基金supported by the National Natural Science Foundation of China(Grant Nos.41201350&41228007)the Knowledge Innovation Program of the Chinese Academy of Sciences(Grant No.Y0S04300KB)
文摘The traditional method of Synthetic Aperture Radar(SAR)wind field retrieval is based on an empirical relation between the near surface winds and the normalized radar backscatter cross section to estimate wind speeds,where this relation is called the geophysical model function(GMF).However,the accuracy rapidly decreases due to the impact of rainfall on the measurement of SAR and the saturation of backscattered intensity under the condition of tropical cyclone.Because of no available instrument synchronously monitoring rain rate on the satellite platform of SAR,we have to derive the precipitation of the SAR observation time from non-simultaneous passive microwave observations of rain in combination with geostationary IR images,and then use the model of rain correction to remove the impact of rain on SAR wind field measurements.For the saturation of radar backscatter cross section in high wind speed conditions,we develop an approach to estimate tropical cyclone parameters and wind fields based on the improved Holland model and the SAR image features of tropical cyclone.To retrieve the low-to-moderate wind speed,the wind direction of tropical cyclone is estimated from the SAR image using wavelet analysis.And then the maximum wind speed and the central pressure of tropical cyclone are calculated by a least square minimization of the difference between the improved Holland model and the low-to-moderate wind speed retrieved from SAR.In addition,wind fields are estimated from the improved Holland model using the above-mentioned parameters of tropical cyclone as input.To evaluate the accuracy of our approach,the SAR images of typhoon Aere,typhoon Khanun,and hurricane Ophelia are used to estimate tropical cyclone parameters and wind fields,which are compared with the best track data and reanalyzed wind fields of the Joint Typhoon Warning Center(JTWC)and the Hurricane Research Division(HRD).The results indicate that the tropical cyclone center,maximum wind speed,and central pressure are generally consistent with the best track data,and wind fields agree well with reanalyzed data from HRD.
基金supported by the National Basic Research Program of China(Grant No.2013CB228604)the National Grand Project for Science and Technology(Grant Nos.2011ZX05030-004-002,2011ZX05019-003&2011ZX05006-002)
文摘Elastic wave inverse scattering theory plays an important role in parameters estimation of heterogeneous media.Combining inverse scattering theory,perturbation theory and stationary phase approximation,we derive the P-wave seismic scattering coefficient equation in terms of fluid factor,shear modulus and density of background homogeneous media and perturbation media.With this equation as forward solver,a pre-stack seismic Bayesian inversion method is proposed to estimate the fluid factor of heterogeneous media.In this method,Cauchy distribution is utilized to the ratios of fluid factors,shear moduli and densities of perturbation media and background homogeneous media,respectively.Gaussian distribution is utilized to the likelihood function.The introduction of constraints from initial smooth models enhances the stability of the estimation of model parameters.Model test and real data example demonstrate that the proposed method is able to estimate the fluid factor of heterogeneous media from pre-stack seismic data directly and reasonably.
