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.展开更多
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.展开更多
To reduce the dependence of EM inversion on the choice of initial model and to obtain the global minimum, we apply transdimensional Bayesian inversion to time-domain airborne electromagnetic data. The transdimensional...To reduce the dependence of EM inversion on the choice of initial model and to obtain the global minimum, we apply transdimensional Bayesian inversion to time-domain airborne electromagnetic data. The transdimensional Bayesian inversion uses the Monte Carlo method to search the model space and yields models that simultaneously satisfy the acceptance probability and data fitting requirements. Finally, we obtain the probability distribution and uncertainty of the model parameters as well as the maximum probability. Because it is difficult to know the height of the transmitting source during flight, we consider a fixed and a variable flight height. Furthermore, we introduce weights into the prior probability density function of the resistivity and adjust the constraint strength in the inversion model by changing the weighing coefficients. This effectively solves the problem of unsatisfactory inversion results in the middle high-resistivity layer. We validate the proposed method by inverting synthetic data with 3% Gaussian noise and field survey data.展开更多
In this paper, the principle of multi-point forming (MPF) technique is presented. One of the most serious defects, wrinkling, during the multi-point forming process of a shallow rectangle cup is discussed by means of ...In this paper, the principle of multi-point forming (MPF) technique is presented. One of the most serious defects, wrinkling, during the multi-point forming process of a shallow rectangle cup is discussed by means of numerical simulation on the shallow rectangle cup forming process. The effects of thickness, material of sheet metal and the pressure of the blank holder are investigated. Based on the simulation results, the reasons and control methods of wrinkling are pointed out. Moreover, the experiment on the multi-point die forming of the shallow rectangle cup by the MPF machine is done to validate the efficiency of the numerical simulation, and the result proves that the application of an elastic cushion in the forming can restrain wrinkling efficiently.展开更多
Based on the absolute and relative gravity observations in North China from 2009 to 2014,spatial dynamic variations of the regional gravity field are obtained. We employed the Euler deconvolution method and the theore...Based on the absolute and relative gravity observations in North China from 2009 to 2014,spatial dynamic variations of the regional gravity field are obtained. We employed the Euler deconvolution method and the theoretical model to get the best estimates of parameters. Gravity field change caused by the depth and distribution in North China is calculated by back analysis. The results show the structural index that equals 1 is suitable for inversion of the gravity variation data. The inversion results indicate that the depths of anomaly field sources are spread over the Hetao fault. The research method of this paper can be used in the quantitative study on the field source and may shed new light on the interpretations of gravity change, and also provide quantitative basis for earthquake prediction index criterions based on the gravity change.展开更多
基金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.
基金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.
基金This paper was financially supported by the Key National Research Project of China (Nos. 2017YFC0601900 and 2016YFC0303100), and the Key Program of National Natural Science Foundation of China (No. 41530320) and Surface Project (No. 41774125).
文摘To reduce the dependence of EM inversion on the choice of initial model and to obtain the global minimum, we apply transdimensional Bayesian inversion to time-domain airborne electromagnetic data. The transdimensional Bayesian inversion uses the Monte Carlo method to search the model space and yields models that simultaneously satisfy the acceptance probability and data fitting requirements. Finally, we obtain the probability distribution and uncertainty of the model parameters as well as the maximum probability. Because it is difficult to know the height of the transmitting source during flight, we consider a fixed and a variable flight height. Furthermore, we introduce weights into the prior probability density function of the resistivity and adjust the constraint strength in the inversion model by changing the weighing coefficients. This effectively solves the problem of unsatisfactory inversion results in the middle high-resistivity layer. We validate the proposed method by inverting synthetic data with 3% Gaussian noise and field survey data.
文摘In this paper, the principle of multi-point forming (MPF) technique is presented. One of the most serious defects, wrinkling, during the multi-point forming process of a shallow rectangle cup is discussed by means of numerical simulation on the shallow rectangle cup forming process. The effects of thickness, material of sheet metal and the pressure of the blank holder are investigated. Based on the simulation results, the reasons and control methods of wrinkling are pointed out. Moreover, the experiment on the multi-point die forming of the shallow rectangle cup by the MPF machine is done to validate the efficiency of the numerical simulation, and the result proves that the application of an elastic cushion in the forming can restrain wrinkling efficiently.
基金funded by the Natural Science Foundation of China(61627824,41274083)the Youth Foundation of Earthquake Prediction(2017010227)
文摘Based on the absolute and relative gravity observations in North China from 2009 to 2014,spatial dynamic variations of the regional gravity field are obtained. We employed the Euler deconvolution method and the theoretical model to get the best estimates of parameters. Gravity field change caused by the depth and distribution in North China is calculated by back analysis. The results show the structural index that equals 1 is suitable for inversion of the gravity variation data. The inversion results indicate that the depths of anomaly field sources are spread over the Hetao fault. The research method of this paper can be used in the quantitative study on the field source and may shed new light on the interpretations of gravity change, and also provide quantitative basis for earthquake prediction index criterions based on the gravity change.