Variation of reservoir physical properties can cause changes in its elastic parameters. However, this is not a simple linear relation. Furthermore, the lack of observations, data overlap, noise interference, and ideal...Variation of reservoir physical properties can cause changes in its elastic parameters. However, this is not a simple linear relation. Furthermore, the lack of observations, data overlap, noise interference, and idealized models increases the uncertainties of the inversion result. Thus, we propose an inversion method that is different from traditional statistical rock physics modeling. First, we use deterministic and stochastic rock physics models considering the uncertainties of elastic parameters obtained by prestack seismic inversion and introduce weighting coefficients to establish a weighted statistical relation between reservoir and elastic parameters. Second, based on the weighted statistical relation, we use Markov chain Monte Carlo simulations to generate the random joint distribution space of reservoir and elastic parameters that serves as a sample solution space of an objective function. Finally, we propose a fast solution criterion to maximize the posterior probability density and obtain reservoir parameters. The method has high efficiency and application potential.展开更多
Spectral decomposition has been widely used in the detection and identifi cation of underground anomalous features(such as faults,river channels,and karst caves).However,the conventional spectral decomposition method ...Spectral decomposition has been widely used in the detection and identifi cation of underground anomalous features(such as faults,river channels,and karst caves).However,the conventional spectral decomposition method is restrained by the window function,and hence,it mostly has low time–frequency focusing and resolution,thereby hampering the fi ne interpretation of seismic targets.To solve this problem,we investigated the sparse inverse spectral decomposition constrained by the lp norm(0<p≤1).Using a numerical model,we demonstrated the higher time–frequency resolution of this method and its capability for improving the seismic interpretation for thin layers.Moreover,given the actual underground geology that can be often complex,we further propose a p-norm constrained inverse spectral attribute interpretation method based on multiresolution time–frequency feature fusion.By comprehensively analyzing the time–frequency spectrum results constrained by the diff erent p-norms,we can obtain more refined interpretation results than those obtained by the traditional strategy,which incorporates a single norm constraint.Finally,the proposed strategy was applied to the processing and interpretation of actual three-dimensional seismic data for a study area covering about 230 km^(2) in western China.The results reveal that the surface water system in this area is characterized by stepwise convergence from a higher position in the north(a buried hill)toward the south and by the development of faults.We thus demonstrated that the proposed method has huge application potential in seismic interpretation.展开更多
基金supported by the National Science and Technology Major Project(No.2011 ZX05007-006)the 973 Program of China(No.2013CB228604)the Major Project of Petrochina(No.2014B-0610)
文摘Variation of reservoir physical properties can cause changes in its elastic parameters. However, this is not a simple linear relation. Furthermore, the lack of observations, data overlap, noise interference, and idealized models increases the uncertainties of the inversion result. Thus, we propose an inversion method that is different from traditional statistical rock physics modeling. First, we use deterministic and stochastic rock physics models considering the uncertainties of elastic parameters obtained by prestack seismic inversion and introduce weighting coefficients to establish a weighted statistical relation between reservoir and elastic parameters. Second, based on the weighted statistical relation, we use Markov chain Monte Carlo simulations to generate the random joint distribution space of reservoir and elastic parameters that serves as a sample solution space of an objective function. Finally, we propose a fast solution criterion to maximize the posterior probability density and obtain reservoir parameters. The method has high efficiency and application potential.
基金supported by National Natural Science Foundation of China (Grant No. 41974140)the PetroChina Prospective,Basic,and Strategic Technology Research Project (No. 2021DJ0606)
文摘Spectral decomposition has been widely used in the detection and identifi cation of underground anomalous features(such as faults,river channels,and karst caves).However,the conventional spectral decomposition method is restrained by the window function,and hence,it mostly has low time–frequency focusing and resolution,thereby hampering the fi ne interpretation of seismic targets.To solve this problem,we investigated the sparse inverse spectral decomposition constrained by the lp norm(0<p≤1).Using a numerical model,we demonstrated the higher time–frequency resolution of this method and its capability for improving the seismic interpretation for thin layers.Moreover,given the actual underground geology that can be often complex,we further propose a p-norm constrained inverse spectral attribute interpretation method based on multiresolution time–frequency feature fusion.By comprehensively analyzing the time–frequency spectrum results constrained by the diff erent p-norms,we can obtain more refined interpretation results than those obtained by the traditional strategy,which incorporates a single norm constraint.Finally,the proposed strategy was applied to the processing and interpretation of actual three-dimensional seismic data for a study area covering about 230 km^(2) in western China.The results reveal that the surface water system in this area is characterized by stepwise convergence from a higher position in the north(a buried hill)toward the south and by the development of faults.We thus demonstrated that the proposed method has huge application potential in seismic interpretation.