Strong heterogeneity and complex pore systems of carbonate reservoir rock make its rock physics model building and fluid substitution difficult and complex. However, rock physics models connect reservoir parameters wi...Strong heterogeneity and complex pore systems of carbonate reservoir rock make its rock physics model building and fluid substitution difficult and complex. However, rock physics models connect reservoir parameters with seismic parameters and fluid substitution is the most effective tool for reservoir prediction and quantitative characterization. On the basis of analyzing complex carbonate reservoir pore structures and heterogeneity at seismic scale, we use the gridding method to divide carbonate rock into homogeneous blocks with independent rock parameters and calculate the elastic moduli of dry rock units step by step using different rock physics models based on pore origin and structural feature. Then, the elastic moduli of rocks saturated with different fluids are obtained using fluid substitution based on different pore connectivity. Based on the calculated elastic moduli of rock units, the Hashin-Shtrikman-Walpole elastic boundary theory is adopted to calculate the carbonate elastic parameters at seismic scale. The calculation and analysis of carbonate models with different combinations of pore types demonstrate the effects of pore type on rock elastic parameters. The simulated result is consistent with our knowledge of real data.展开更多
Shaley sandstone is heterogeneous at a seismic scale. Gassmann＇s equation is suited for fluid substitution in a homogeneous medium. To study the difference between shaley sandstone effective elastic moduli calculated...Shaley sandstone is heterogeneous at a seismic scale. Gassmann＇s equation is suited for fluid substitution in a homogeneous medium. To study the difference between shaley sandstone effective elastic moduli calculated by mean porosity as a homogeneous medium, and those calculated directly from the sub-volumes of the volume as a heterogeneous medium, computational experiments are conducted on Han＇s shaley sand model, the soft-sand model, the stiff-sand model, and their combination under the assumption that the shaley sandstone volume is made up of separate homogenous sub-volumes with independent porosity and clay content. Fluid substitutions are conducted by Gassmann＇s equation on rock volume and sub-volumes respectively. The computational data show that at seismic scale, there are minor differences between fluid substitution on rock volume and that on sub-volumes using Gassmann＇s equation. But fluid substitution on sub-volumes can take consideration of the effects of low porosity and low permeability sub-volumes, which can get more reasonable data, especially for low porosity reservoirs.展开更多
The Qilian Orogen Zone(QOZ), located in the north margin of the Tibetan Plateau, is the key area for understanding the deformation and dynamics process of Tibet. Numerous geological and geophysical studies have been c...The Qilian Orogen Zone(QOZ), located in the north margin of the Tibetan Plateau, is the key area for understanding the deformation and dynamics process of Tibet. Numerous geological and geophysical studies have been carried out on the mechanics of the Tibetan Plateau deformation and uplift; however, the detailed structure and deformation style of the Qilian Orogen Zone have remained uncertain due to poor geophysical data coverage and limited resolution power of inversion algorithms. In this study, we analyze the P-wave velocity structure beneath the Qilian Orogen Zone, obtained by applying multi-scale seismic tomography technique to P-wave arrival time data recorded by regional seismic networks. The seismic tomography algorithm used in this study employs sparsity constraints on the wavelet representation of the velocity model via L1-norm regularization. This algorithm can deal efficiently with uneven-sampled volumes, and can obtain multi-scale images of the velocity model. Our results can be summarized as follows:(1) The crustal velocity structure is strongly inhomogeneous and consistent with the surface geological setting. Significant low-velocity anomalies exist in the crust of northeastern Tibet, and slight high-velocity anomalies exist beneath the Qaidam Basin and Alxa terrane.(2)The Qilian Orogen Zone can be divided into two main parts by the Laji Shan Faults: the northwestern part with a low-velocity feature, and the southeastern part with a high-velocity feature at the upper and middle crust.(3) Our tomographic images suggest that northwestern and southeastern Qilian Orogen Zones have undergone different tectonic processes. In the northwest Qilian Orogen Zone, the deformation and growth of the Northern Tibetan Plateau has extended to the Heli Shan and Beida Shan region by northward overthrusting at the upper crust and thickening in the lower crust. We speculate that in the southeast Qilian Orogen Zone the deformation and growth of the Northern Tibet Plateau were of strike-slip style at the upper crust; in the lower crust, the evidence suggests ductile shear extrusion style and active frontage extension to the Alxa terrane.(4) The multi-scale seismic tomography technique provides multiscale analysis and sparse constraints, which has allowed to us obtain stable, high-resolution results.展开更多
Seismic energy decays while propagating subsurface, which may reduce the resolution of seismic data. This paper studies the method of seismic energy dispersion compensation which provides the basic principles for mult...Seismic energy decays while propagating subsurface, which may reduce the resolution of seismic data. This paper studies the method of seismic energy dispersion compensation which provides the basic principles for multi-scale morphology and the spectrum simulation method. These methods are applied in seismic energy compensation. First of all, the seismic data is decomposed into multiple scales and the effective frequency bandwidth is selectively broadened for some scales by using a spectrum simulation method. In this process, according to the amplitude spectrum of each scale, the best simulation range is selected to simulate the middle and low frequency components to ensure the authenticity of the simulation curve which is calculated by the median method, and the high frequency component is broadened. Finally, these scales are reconstructed with reasonable coefficients, and the compensated seismic data can be obtained. Examples are shown to illustrate the feasibility of the energy compensation method.展开更多
Bouguer gravity anomaly in North China is decomposed with multi scale decomposition technique of wavelet transform. Gravity anomalies produced by anomalous density bodies of various scales are revealed from surface to...Bouguer gravity anomaly in North China is decomposed with multi scale decomposition technique of wavelet transform. Gravity anomalies produced by anomalous density bodies of various scales are revealed from surface to Moho. Characteristics of anomalies of different orders and corresponding structural features are discussed. The result shows that details of wavelet transform of different orders reflect the distribution features of rock density at different depths and in various scales. In most cases, the two sides of a fault especially a deep and large fault in North China differ greatly in rock density. This difference records the history of the formation and evolution of the crust. Deep structural setting for the \%M\%s≥7.0 strong earthquakes in this region is also discussed.展开更多
基金sponsored jointly by the National Natural Science Foundation of China(No.41074098)the Key State Science and Technology Project(2011ZX05023-005-005)China University of Petroleum(Beijing) Fund(KYJJ2012-05-08)
文摘Strong heterogeneity and complex pore systems of carbonate reservoir rock make its rock physics model building and fluid substitution difficult and complex. However, rock physics models connect reservoir parameters with seismic parameters and fluid substitution is the most effective tool for reservoir prediction and quantitative characterization. On the basis of analyzing complex carbonate reservoir pore structures and heterogeneity at seismic scale, we use the gridding method to divide carbonate rock into homogeneous blocks with independent rock parameters and calculate the elastic moduli of dry rock units step by step using different rock physics models based on pore origin and structural feature. Then, the elastic moduli of rocks saturated with different fluids are obtained using fluid substitution based on different pore connectivity. Based on the calculated elastic moduli of rock units, the Hashin-Shtrikman-Walpole elastic boundary theory is adopted to calculate the carbonate elastic parameters at seismic scale. The calculation and analysis of carbonate models with different combinations of pore types demonstrate the effects of pore type on rock elastic parameters. The simulated result is consistent with our knowledge of real data.
基金supported by National Natural Science Function of China (No. 41074098)National 973 Basic Research Program (No. 2007CB209606)
文摘Shaley sandstone is heterogeneous at a seismic scale. Gassmann＇s equation is suited for fluid substitution in a homogeneous medium. To study the difference between shaley sandstone effective elastic moduli calculated by mean porosity as a homogeneous medium, and those calculated directly from the sub-volumes of the volume as a heterogeneous medium, computational experiments are conducted on Han＇s shaley sand model, the soft-sand model, the stiff-sand model, and their combination under the assumption that the shaley sandstone volume is made up of separate homogenous sub-volumes with independent porosity and clay content. Fluid substitutions are conducted by Gassmann＇s equation on rock volume and sub-volumes respectively. The computational data show that at seismic scale, there are minor differences between fluid substitution on rock volume and that on sub-volumes using Gassmann＇s equation. But fluid substitution on sub-volumes can take consideration of the effects of low porosity and low permeability sub-volumes, which can get more reasonable data, especially for low porosity reservoirs.
基金supported by the National Natural Science Foundation of China(41574045,41590862)State Key Laboratory of Earthquake Dynamics(LED2013A06)
文摘The Qilian Orogen Zone(QOZ), located in the north margin of the Tibetan Plateau, is the key area for understanding the deformation and dynamics process of Tibet. Numerous geological and geophysical studies have been carried out on the mechanics of the Tibetan Plateau deformation and uplift; however, the detailed structure and deformation style of the Qilian Orogen Zone have remained uncertain due to poor geophysical data coverage and limited resolution power of inversion algorithms. In this study, we analyze the P-wave velocity structure beneath the Qilian Orogen Zone, obtained by applying multi-scale seismic tomography technique to P-wave arrival time data recorded by regional seismic networks. The seismic tomography algorithm used in this study employs sparsity constraints on the wavelet representation of the velocity model via L1-norm regularization. This algorithm can deal efficiently with uneven-sampled volumes, and can obtain multi-scale images of the velocity model. Our results can be summarized as follows:(1) The crustal velocity structure is strongly inhomogeneous and consistent with the surface geological setting. Significant low-velocity anomalies exist in the crust of northeastern Tibet, and slight high-velocity anomalies exist beneath the Qaidam Basin and Alxa terrane.(2)The Qilian Orogen Zone can be divided into two main parts by the Laji Shan Faults: the northwestern part with a low-velocity feature, and the southeastern part with a high-velocity feature at the upper and middle crust.(3) Our tomographic images suggest that northwestern and southeastern Qilian Orogen Zones have undergone different tectonic processes. In the northwest Qilian Orogen Zone, the deformation and growth of the Northern Tibetan Plateau has extended to the Heli Shan and Beida Shan region by northward overthrusting at the upper crust and thickening in the lower crust. We speculate that in the southeast Qilian Orogen Zone the deformation and growth of the Northern Tibet Plateau were of strike-slip style at the upper crust; in the lower crust, the evidence suggests ductile shear extrusion style and active frontage extension to the Alxa terrane.(4) The multi-scale seismic tomography technique provides multiscale analysis and sparse constraints, which has allowed to us obtain stable, high-resolution results.
文摘Seismic energy decays while propagating subsurface, which may reduce the resolution of seismic data. This paper studies the method of seismic energy dispersion compensation which provides the basic principles for multi-scale morphology and the spectrum simulation method. These methods are applied in seismic energy compensation. First of all, the seismic data is decomposed into multiple scales and the effective frequency bandwidth is selectively broadened for some scales by using a spectrum simulation method. In this process, according to the amplitude spectrum of each scale, the best simulation range is selected to simulate the middle and low frequency components to ensure the authenticity of the simulation curve which is calculated by the median method, and the high frequency component is broadened. Finally, these scales are reconstructed with reasonable coefficients, and the compensated seismic data can be obtained. Examples are shown to illustrate the feasibility of the energy compensation method.
文摘Bouguer gravity anomaly in North China is decomposed with multi scale decomposition technique of wavelet transform. Gravity anomalies produced by anomalous density bodies of various scales are revealed from surface to Moho. Characteristics of anomalies of different orders and corresponding structural features are discussed. The result shows that details of wavelet transform of different orders reflect the distribution features of rock density at different depths and in various scales. In most cases, the two sides of a fault especially a deep and large fault in North China differ greatly in rock density. This difference records the history of the formation and evolution of the crust. Deep structural setting for the \%M\%s≥7.0 strong earthquakes in this region is also discussed.
