A critical porosity model is often used to calculate the dry frame elastic modulus by the rock critical porosity value which is affected by many factors. In practice it is hard for us to obtain an accurate critical po...A critical porosity model is often used to calculate the dry frame elastic modulus by the rock critical porosity value which is affected by many factors. In practice it is hard for us to obtain an accurate critical porosity value and we can generally take only an empirical critical porosity value which often causes errors. In this paper, we propose a method to obtain the rock critical porosity value by inverting P-wave velocity and applying it to predict S-wave velocity. The applications of experiment and log data both show that the critical porosity inversion method can reduce the uncertainty resulting from using an empirical value in the past and provide the accurate critical porosity value for predicting S-wave velocity which significantly improves the prediction accuracy.展开更多
Rock-physics models are constructed for hydrate-bearing sediments in the Qilian Mountains permafrost region using the K–T equation model, and modes I and II of the effective medium model. The K–T equation models the...Rock-physics models are constructed for hydrate-bearing sediments in the Qilian Mountains permafrost region using the K–T equation model, and modes I and II of the effective medium model. The K–T equation models the seismic wave propagation in a two-phase medium to determine the elastic moduli of the composite medium. In the effective medium model, mode I, the hydrate is a component of the pore inclusions in mode I and in mode II it is a component of the matrix. First, the P-wave velocity, S-wave velocity, density, bulk modulus, and shear modulus of the sediment matrix are extracted from logging data.. Second, based on the physical properties of the main components of the sediments, rock-physics model is established using the K–T equation, and two additional rock-physics models are established assuming different hydrate-filling modes for the effective medium. The model and actual velocity data for the hydrate-bearing sediments are compared and it is found that the rock-physics model for the hydrate-filling mode II well reproduces the actual data.展开更多
By substituting rock skeleton modulus expressions into Gassmann approximate fluid equation, we obtain a seismic porosity inversion equation. However, conventional rock skeleton models and their expressions are quite d...By substituting rock skeleton modulus expressions into Gassmann approximate fluid equation, we obtain a seismic porosity inversion equation. However, conventional rock skeleton models and their expressions are quite different from each other, resuling in different seismic porosity inversion equations, potentially leading to difficulties in correctly applying them and evaluating their results. In response to this, a uniform relation with two adjusting parameters suitable for all rock skeleton models is established from an analysis and comparison of various conventional rock skeleton models and their expressions including the Eshelby-Walsh, Pride, Geertsma, Nur, Keys-Xu, and Krief models. By giving the two adjusting parameters specific values, different rock skeleton models with specific physical characteristics can be generated. This allows us to select the most appropriate rock skeleton model based on geological and geophysical conditions, and to develop more wise seismic porosity inversion. As an example of using this method for hydrocarbon prediction and fluid identification, we apply this improved porosity inversion, associated with rock physical data and well log data, to the ZJ basin. Research shows that the existence of an abundant hydrocarbon reservoir is dependent on a moderate porosity range, which means we can use the results of seismic porosity inversion to identify oil reservoirs and dry or water-saturated reservoirs. The seismic inversion results are closely correspond to well log porosity curves in the ZJ area, indicating that the uniform relations and inversion methods proposed in this paper are reliable and effective.展开更多
Matrix porosity calculations of fractured and vuggy reservoirs, such as volcanics and weathered dolomite, are one of the problems urgently needed to solve in well-log evaluation. In this paper, we first compare the an...Matrix porosity calculations of fractured and vuggy reservoirs, such as volcanics and weathered dolomite, are one of the problems urgently needed to solve in well-log evaluation. In this paper, we first compare the an empirical formula for porosity calculation from full diameter rhyolite core experiments with the matrix porosity formulas commonly used. We discuss the applicability of the empirical formula in fractured and vuggy reservoirs, such as intermediate-basic volcanics and weathered dolomite. Based on core analysis data, the error distribution of the calculated porosity of our empirical formula and the other porosity formulas in these reservoirs are given. The statistical error analysis indicates that the our empirical formula provides a higher precision than the other porosity formulas. When the porosity is between 1.5% and 15%, the acoustic experiment formula can be used not only for acidic volcanics but also in other fractured and vuggy reservoirs, such as intermediate-basic volcanics and weathered dolomite. Moreover, the formula can reduce the effects of borehole enlargement and rock alteration on porosity computation.展开更多
Homogeneity and heterogeneity are two totally different concepts in nature.At the particle length scale,rocks exhibit strong heterogeneity in their constituents and porosities.When the heterogeneity of porosity obeys ...Homogeneity and heterogeneity are two totally different concepts in nature.At the particle length scale,rocks exhibit strong heterogeneity in their constituents and porosities.When the heterogeneity of porosity obeys the random uniform distribution,both the mean value and the variance of porosities in the heterogeneous porosity field can be used to reflect the overall heterogeneous characteristics of the porosity field.The main purpose of this work is to investigate the effects of porosity heterogeneity on chemical dissolution front instability in fluid-saturated rocks by the computational simulation method.The related computational simulation results have demonstrated that:1) since the propagation speed of a chemical dissolution front is inversely proportional to the difference between the final porosity and the mean value of porosities in the initial porosity field,an increase in the extent of the porosity heterogeneity can cause an increase in the mean value of porosities in the initial porosity field and an increase in the propagation speed of the chemical dissolution front.2) An increase in the variance of porosities in the initial porosity field can cause an increase in the instability probability of the chemical dissolution front in the fluid-saturated rock.3) The greater the mean value of porosities in the initial porosity field,the quicker the irregular morphology of the chemical dissolution front changes in the supercritical chemical dissolution systems.This means that the irregular morphology of a chemical dissolution front grows quicker in a porosity field of heterogeneity than it does in that of homogeneity when the chemical dissolution system is at a supercritical stage.展开更多
基金sponsored by Important National Science and Technology Specifi c Projects of China (No.2011ZX05001)
文摘A critical porosity model is often used to calculate the dry frame elastic modulus by the rock critical porosity value which is affected by many factors. In practice it is hard for us to obtain an accurate critical porosity value and we can generally take only an empirical critical porosity value which often causes errors. In this paper, we propose a method to obtain the rock critical porosity value by inverting P-wave velocity and applying it to predict S-wave velocity. The applications of experiment and log data both show that the critical porosity inversion method can reduce the uncertainty resulting from using an empirical value in the past and provide the accurate critical porosity value for predicting S-wave velocity which significantly improves the prediction accuracy.
基金supported by the Institute of Geophysical and Geochemical Exploration(IGGE)CAGS of China(No.WH201207)
文摘Rock-physics models are constructed for hydrate-bearing sediments in the Qilian Mountains permafrost region using the K–T equation model, and modes I and II of the effective medium model. The K–T equation models the seismic wave propagation in a two-phase medium to determine the elastic moduli of the composite medium. In the effective medium model, mode I, the hydrate is a component of the pore inclusions in mode I and in mode II it is a component of the matrix. First, the P-wave velocity, S-wave velocity, density, bulk modulus, and shear modulus of the sediment matrix are extracted from logging data.. Second, based on the physical properties of the main components of the sediments, rock-physics model is established using the K–T equation, and two additional rock-physics models are established assuming different hydrate-filling modes for the effective medium. The model and actual velocity data for the hydrate-bearing sediments are compared and it is found that the rock-physics model for the hydrate-filling mode II well reproduces the actual data.
基金supported by the National Nature Science Foundation of China(Grant No.41174114)Important National Science and Technology Specific Projects(Grant No.2011ZX05025-005-010)
文摘By substituting rock skeleton modulus expressions into Gassmann approximate fluid equation, we obtain a seismic porosity inversion equation. However, conventional rock skeleton models and their expressions are quite different from each other, resuling in different seismic porosity inversion equations, potentially leading to difficulties in correctly applying them and evaluating their results. In response to this, a uniform relation with two adjusting parameters suitable for all rock skeleton models is established from an analysis and comparison of various conventional rock skeleton models and their expressions including the Eshelby-Walsh, Pride, Geertsma, Nur, Keys-Xu, and Krief models. By giving the two adjusting parameters specific values, different rock skeleton models with specific physical characteristics can be generated. This allows us to select the most appropriate rock skeleton model based on geological and geophysical conditions, and to develop more wise seismic porosity inversion. As an example of using this method for hydrocarbon prediction and fluid identification, we apply this improved porosity inversion, associated with rock physical data and well log data, to the ZJ basin. Research shows that the existence of an abundant hydrocarbon reservoir is dependent on a moderate porosity range, which means we can use the results of seismic porosity inversion to identify oil reservoirs and dry or water-saturated reservoirs. The seismic inversion results are closely correspond to well log porosity curves in the ZJ area, indicating that the uniform relations and inversion methods proposed in this paper are reliable and effective.
基金sponsored by the Science Research and Technology Development Project of Petrochina Company Limited "Well Logging Interpretation and Integrative Evaluation of the Complex Lithology"(Grant No.2008A-2705)
文摘Matrix porosity calculations of fractured and vuggy reservoirs, such as volcanics and weathered dolomite, are one of the problems urgently needed to solve in well-log evaluation. In this paper, we first compare the an empirical formula for porosity calculation from full diameter rhyolite core experiments with the matrix porosity formulas commonly used. We discuss the applicability of the empirical formula in fractured and vuggy reservoirs, such as intermediate-basic volcanics and weathered dolomite. Based on core analysis data, the error distribution of the calculated porosity of our empirical formula and the other porosity formulas in these reservoirs are given. The statistical error analysis indicates that the our empirical formula provides a higher precision than the other porosity formulas. When the porosity is between 1.5% and 15%, the acoustic experiment formula can be used not only for acidic volcanics but also in other fractured and vuggy reservoirs, such as intermediate-basic volcanics and weathered dolomite. Moreover, the formula can reduce the effects of borehole enlargement and rock alteration on porosity computation.
基金Project(11272359)supported by the National Natural Science Foundation of China
文摘Homogeneity and heterogeneity are two totally different concepts in nature.At the particle length scale,rocks exhibit strong heterogeneity in their constituents and porosities.When the heterogeneity of porosity obeys the random uniform distribution,both the mean value and the variance of porosities in the heterogeneous porosity field can be used to reflect the overall heterogeneous characteristics of the porosity field.The main purpose of this work is to investigate the effects of porosity heterogeneity on chemical dissolution front instability in fluid-saturated rocks by the computational simulation method.The related computational simulation results have demonstrated that:1) since the propagation speed of a chemical dissolution front is inversely proportional to the difference between the final porosity and the mean value of porosities in the initial porosity field,an increase in the extent of the porosity heterogeneity can cause an increase in the mean value of porosities in the initial porosity field and an increase in the propagation speed of the chemical dissolution front.2) An increase in the variance of porosities in the initial porosity field can cause an increase in the instability probability of the chemical dissolution front in the fluid-saturated rock.3) The greater the mean value of porosities in the initial porosity field,the quicker the irregular morphology of the chemical dissolution front changes in the supercritical chemical dissolution systems.This means that the irregular morphology of a chemical dissolution front grows quicker in a porosity field of heterogeneity than it does in that of homogeneity when the chemical dissolution system is at a supercritical stage.
