Total organic carbon (TOC) prediction with elastic parameter inversions has been widely used in the identification and evaluation of source rocks. However, the elastic parameters used to predict TOC are not only deter...Total organic carbon (TOC) prediction with elastic parameter inversions has been widely used in the identification and evaluation of source rocks. However, the elastic parameters used to predict TOC are not only determined by TOC but also depend on the other physical properties of source rocks. Besides, the TOC prediction with the elastic parameters inversion is an indirect method based on the statistical relationship obtained from well logs and experiment data. Therefore, we propose a rock physics model and define a TOC indicator mainly affected by TOC to predict TOC directly. The proposed rock physics model makes the equivalent elastic moduli of source rocks parameterized by the TOC indicator. Combining the equivalent elastic moduli of source rocks and Gray’s approximation leads to a novel linearized approximation of the P-wave reflection coefficient incorporating the TOC indicator. Model examples illustrate that the novel reflectivity approximation well agrees with the exact Zoeppritz equation until incident angles reach 40°. Convoluting the novel P-wave reflection approximation with seismic wavelets as the forward solver, an AVO inversion method based on the Bayesian theory is proposed to invert the TOC indicator with seismic data. The synthetic examples and field tests validate the feasibility and stability of the proposed AVO inversion approach. Using the inversion results of the TOC indicator, TOC is directly and accurately estimated in the target area.展开更多
Rock physics inversion is to use seismic elastic properties of underground strata for predicting reservoir petrophysical parameters.The Markov chain Monte Carlo(MCMC)algorithm is commonly used to solve rock physics in...Rock physics inversion is to use seismic elastic properties of underground strata for predicting reservoir petrophysical parameters.The Markov chain Monte Carlo(MCMC)algorithm is commonly used to solve rock physics inverse problems.However,all the parameters to be inverted are iterated simultaneously in the conventional MCMC algorithm.What is obtained is an optimal solution of combining the petrophysical parameters with being inverted.This study introduces the alternating direction(AD)method into the MCMC algorithm(i.e.the optimized MCMC algorithm)to ensure that each petrophysical parameter can get the optimal solution and improve the convergence of the inversion.Firstly,the Gassmann equations and Xu-White model are used to model shaly sandstone,and the theoretical relationship between seismic elastic properties and reservoir petrophysical parameters is established.Then,in the framework of Bayesian theory,the optimized MCMC algorithm is used to generate a Markov chain to obtain the optimal solution of each physical parameter to be inverted and obtain the maximum posterior density of the physical parameter.The proposed method is applied to actual logging and seismic data and the results show that the method can obtain more accurate porosity,saturation,and clay volume.展开更多
Kerogen content and kerogen porosity play a significant role in elastic properties of organic-rich shales. We construct a rock physics model for organic-rich shales to quantify the effect of kerogen content and keroge...Kerogen content and kerogen porosity play a significant role in elastic properties of organic-rich shales. We construct a rock physics model for organic-rich shales to quantify the effect of kerogen content and kerogen porosity using the Kuster and Toks6z theory and the selfconsistent approximation method. Rock physics modeling results show that with the increase of kerogen content and kerogen-related porosity, the velocity and density of shales decrease, and the effect of kerogen porosity becomes more obvious only for higher kerogen content. We also find that the Poisson's ratio of the shale is not sensitive to kerogen porosity for the case of gas saturation. Finally, for the seismic reflection responses of an organic-rich shale layer, forward modeling results indicate the fifth type AVO re- sponses which correspond to a negative intercept and a positive gradient. The absolute values of intercept and gradient increase with kerogen content and kerogen porosity, and present predictable variations associated with velocities and density.展开更多
Due to the huge differences between the unconventional shale and conventional sand reservoirs in many aspects such as the types and the characteristics of minerals,matrix pores and fluids,the construction of shale roc...Due to the huge differences between the unconventional shale and conventional sand reservoirs in many aspects such as the types and the characteristics of minerals,matrix pores and fluids,the construction of shale rock physics model is significant for the exploration and development of shale reservoirs.To make a better characterization of shale gas-bearing reservoirs,we first propose a new but more suitable rock physics model to characterize the reservoirs.We then use a well A to demonstrate the feasibility and reliability of the proposed rock physics model of shale gas-bearing reservoirs.Moreover,we propose a new brittleness indicator for the high-porosity and organic-rich shale gas-bearing reservoirs.Based on the parameter analysis using the constructed rock physics model,we finally compare the new brittleness indicator with the commonly used Young’s modulus in the content of quartz and organic matter,the matrix porosity,and the types of filled fluids.We also propose a new shale brittleness index by integrating the proposed new brittleness indicator and the Poisson’s ratio.Tests on real data sets demonstrate that the new brittleness indicator and index are more sensitive than the commonly used Young’s modulus and brittleness index for the high-porosity and high-brittleness shale gas-bearing reservoirs.展开更多
Digital rock physics(DRP)is a paramount technology to improve the economic benefits of oil and gas fields,devise more scientific oil and gas field development plans,and create digital oil and gas fields.Currently,a si...Digital rock physics(DRP)is a paramount technology to improve the economic benefits of oil and gas fields,devise more scientific oil and gas field development plans,and create digital oil and gas fields.Currently,a significant gap is present between DRP theory and practical applications.Conventional digital-core construction focuses only on simple cores,and the recognition and segmentation effect of fractures and pores of complex cores is poor.The identification of rock minerals is inaccurate,which leads to the difference between the digital and actual cores.To promote the application of DRP in developing oil and gas fields,based on the high-precision X-ray computed tomography scanning technology,the U-Net deep learning model of the full convolution neural network is used to segment the pores,fractures,and matrix from the complex rock core with natural fractures innovatively.Simultaneously,the distribution of rock minerals is divided,and the distribution of rock conditions is corrected by X-ray diffraction.A pore—fracture network model is established based on the equivalent radius,which lays the foundation for fluid seepage simulation.Finally,the accuracy of the established a digital core is verified by the porosity measured via nuclear magnetic resonance technology,which is of great significance to the development and application of DRP in oil and gas fields.展开更多
Tight sandstone gas reservoirs have characteristics of low porosity and low permeability,complex pore structure.In the paper,we consider the pore and micro-fracture systems in tight sandstone reservoir and perform roc...Tight sandstone gas reservoirs have characteristics of low porosity and low permeability,complex pore structure.In the paper,we consider the pore and micro-fracture systems in tight sandstone reservoir and perform rock physics modeling of sandstone gas reservoirs based on Chapman.Then we analyze the effects of relaxation time,and difference fluid types on dispersion and attention,and seismic responses.We found that the magnitude of the relaxation time is related to the gas saturation and the viscosity coefficient of the fluid.The magnitude of the relaxation time affects the frequency range of the velocity dispersion and attenuation,and the type of fluid affects the magnitude of the velocity dispersion and attenuation.As the relaxation time decreases,the calculated amplitude deceases.And the amplitude and the waveform of the synthetic seismic record are related to different fluid types and formation thickness.展开更多
Heavy oil has high density and viscosity, and exhibits viscoelasticity. Gassmann's theory is not suitable for materials saturated with viscoelastic fluids. Directly applying such model leads to unreliable results ...Heavy oil has high density and viscosity, and exhibits viscoelasticity. Gassmann's theory is not suitable for materials saturated with viscoelastic fluids. Directly applying such model leads to unreliable results for seismic inversion of heavy oil reservoir. To describe the viscoelastic behavior of heavy oil, we modeled the elastic properties of heavy oil with varying viscosity and frequency using the Cole-Cole-Maxwell (CCM) model. Then, we used a CCoherent Potential Approximation (CPA) instead of the Gassmann equations to account for the fluid effect, by extending the single-phase fluid condition to two-phase fluid (heavy oil and water) condition, so that partial saturation of heavy oil can be considered. This rock physics model establishes the relationship between the elastic modulus of reservoir rock and viscosity, frequency and saturation. The viscosity of the heavy oil and the elastic moduli and porosity of typical reservoir rock samples were measured in laboratory, which were used for calibration of the rock physics model. The well-calibrated frequency-variant CPA model was applied to the prediction of the P- and S-wave velocities in the seismic frequency range (1–100 Hz) and the inversion of petrophysical parameters for a heavy oil reservoir. The pre-stack inversion results of elastic parameters are improved compared with those results using the CPA model in the sonic logging frequency (∼10 kHz), or conventional rock physics model such as the Xu-Payne model. In addition, the inversion of the porosity of the reservoir was conducted with the simulated annealing method, and the result fits reasonably well with the logging curve and depicts the location of the heavy oil reservoir on the time slice. The application of the laboratory-calibrated CPA model provides better results with the velocity dispersion correction, suggesting the important role of accurate frequency dependent rock physics models in the seismic prediction of heavy oil reservoirs.展开更多
Seismic characterizing of tight gas sandstone (TGS) reservoirs is essential for identifying promising gas-bearing regions. However, exploring the petrophysical significance of seismic-inverted elastic properties is ch...Seismic characterizing of tight gas sandstone (TGS) reservoirs is essential for identifying promising gas-bearing regions. However, exploring the petrophysical significance of seismic-inverted elastic properties is challenging due to the complex microstructures in TGSs. Meanwhile, interbedded structures of sandstone and mudstone intensify the difficulty in accurately extracting the crucial tight sandstone properties. An integrated rock-physics-based framework is proposed to estimate the reservoir quality of TGSs from seismic data. TGSs with complex pore structures are modeled using the double-porosity model, providing a practical tool to compute rock physics templates for reservoir parameter estimation. The VP/VS ratio is utilized to predict the cumulative thickness of the TGS reservoirs within the target range via the threshold value evaluated from wireline logs for lithology discrimination. This approach also facilitates better capturing the elastic properties of the TGSs for quantitative seismic interpretation. Total porosity is estimated from P-wave impedance using the correlation obtained based on wireline log analysis. After that, the three-dimensional rock-physics templates integrated with the estimated total porosity are constructed to interpret microfracture porosity and gas saturation from velocity ratio and bulk modulus. The integrated framework can optimally estimate the parameters dominating the reservoir quality. The results of the indicator proposed based on the obtained parameters are in good agreement with the gas productions and can be utilized to predict promising TGS reservoirs. Moreover, the results suggest that considering microfracture porosity allows a more accurate prediction of high-quality reservoirs, further validating the applicability of the proposed method in the studied region.展开更多
Knowledge about the seismic elastic modulus dispersion,and associated attenuation,in fluid-saturated rocks is essential for better interpretation of seismic observations taken as part of hydrocarbon identification and...Knowledge about the seismic elastic modulus dispersion,and associated attenuation,in fluid-saturated rocks is essential for better interpretation of seismic observations taken as part of hydrocarbon identification and time-lapse seismic surveillance of both conventional and unconventional reservoir and overburden performances.A Seismic Elastic Moduli Module has been developed,based on the forced-oscillations method,to experimentally investigate the frequency dependence of Young's modulus and Poisson's ratio,as well as the inferred attenuation,of cylindrical samples under different confining pressure conditions.Calibration with three standard samples showed that the measured elastic moduli were consistent with the published data,indicating that the new apparatus can operate reliably over a wide frequency range of f∈[1-2000,10^(6)]Hz.The Young's modulus and Poisson's ratio of the shale and the tight sandstone samples were measured under axial stress oscillations to assess the frequency-and pressure-dependent effects.Under dry condition,both samples appear to be nearly frequency independent,with weak pressure dependence for the shale and significant pressure dependence for the sandstone.In particular,it was found that the tight sandstone with complex pore microstructure exhibited apparent dispersion and attenuation under brine or glycerin saturation conditions,the levels of which were strongly influenced by the increased effective pressure.In addition,the measured Young's moduli results were compared with the theoretical predictions from a scaled poroelastic model with a reasonably good agreement,revealing that the combined fluid flow mechanisms at both mesoscopic and microscopic scales possibly responsible for the measured dispersion.展开更多
A considerable effort has been made in the literature for quality assurance (QA) and quality checking (QC) of the petrophysical log data for computation of reservoir rock property parameters. Well log data plays an in...A considerable effort has been made in the literature for quality assurance (QA) and quality checking (QC) of the petrophysical log data for computation of reservoir rock property parameters. Well log data plays an integral role in building a rock physics model for quantitative interpretation (QI) work. A poor-quality rock physics model may lead to significant financial and HSSE implications by drilling wells in undesired locations. Historically, a variety of techniques have been used including histograms and cross plots for reviewing the feasibility of petrophysical logs for QI work. However, no attempt has ever been made to introduce a simplified workflow. This paper serves two-fold. It provides a simplified step by step approach for building a petrophysics/rock physics model. A case study has been presented to compare the synthetic seismogram generated from the simplified workflow with the actual seismic trace at well locations. Secondly, the paper shows how a few key cross plots and rock property parameters provide adequate information to validate petrophysical data, distinguish overburden and reservoir sections, and to help identify fluids and saturation trends within the reservoir sands. In the mentioned case study, the robustness of the simplified rock physics model has helped seismic data to successfully distinguish hydrocarbon bearing reservoir sands from non-reservoir shales.展开更多
Seismic Rock physics plays a bridge role between the rock moduli and physical properties of the hydrocarbon reservoirs.Prestack seismic inversion is an important method for the quantitative characterization of elastic...Seismic Rock physics plays a bridge role between the rock moduli and physical properties of the hydrocarbon reservoirs.Prestack seismic inversion is an important method for the quantitative characterization of elasticity,physical properties,lithology and fluid properties of subsurface reservoirs.In this paper,a high order approximation of rock physics model for clastic rocks is established and one seismic AVO reflection equation characterized by the high order approximation(Jacobian and Hessian matrix)of rock moduli is derived.Besides,the contribution of porosity,shale content and fluid saturation to AVO reflectivity is analyzed.The feasibility of the proposed AVO equation is discussed in the direct estimation of rock physical properties.On the basis of this,one probabilistic AVO inversion based on differential evolution-Markov chain Monte Carlo stochastic model is proposed on the premise that the model parameters obey Gaussian mixture probability prior model.The stochastic model has both the global optimization characteristics of the differential evolution algorithm and the uncertainty analysis ability of Markov chain Monte Carlo model.Through the cross parallel of multiple Markov chains,multiple stochastic solutions of the model parameters can be obtained simultaneously,and the posterior probability density distribution of the model parameters can be simulated effectively.The posterior mean is treated as the optimal solution of the model to be inverted.Besides,the variance and confidence interval are utilized to evaluate the uncertainties of the estimated results,so as to realize the simultaneous estimation of reservoir elasticity,physical properties,discrete lithofacies and dry rock skeleton.The validity of the proposed approach is verified by theoretical tests and one real application case in eastern China.展开更多
This paper reviews the recent progress in the studies of experimental rock mechanics and tectonophysics concerning seismology and physics of the Earths interior in China. The progress is presented in the following asp...This paper reviews the recent progress in the studies of experimental rock mechanics and tectonophysics concerning seismology and physics of the Earths interior in China. The progress is presented in the following aspects: a) A lot of results of experiment and numerical simulation enrich our knowledge of the brittle fracturing process under the condition with heterogeneity in material and structure; b) Some new results on frictional behavior of non-homogeneous faults reveal the complexity of faulting behavior; c) Some new results on the brittle-plastic transition and plastic flow are obtained; especially the important progress is obtained on rheological properties of rocks in the lower crust and the upper mantle; d) A lot of experimental results are obtained on rock physics at high temperature and pressure and have been used in study of material composition and state. These results provide useful information for understanding the physical properties and deformation mechanisms of material of the Earths interior and earthquake physics.展开更多
Simultaneous waveform inversion was used to predict lithofacies and fluid type across the field. Very often, characterizing reservoirs in terms of lithology and fluid type using conventional methods is replete with un...Simultaneous waveform inversion was used to predict lithofacies and fluid type across the field. Very often, characterizing reservoirs in terms of lithology and fluid type using conventional methods is replete with uncertainties, especially in marginal fields. An approach is employed in this study that integrated rock physics and waveform inverse modelling for lithology and fluid-type characterization to appropriately identify potential hydrocarbon saturated zones and their corresponding lithology. Seismic and well-log data were analyzed using Hampson Russel software. The method adopted includes lithofacies and fluid content analysis using rock physics parameters and seismic simultaneous inverse modelling. Rock physics analysis identified 2 broad reservoirs namely: HDZ1 and HDZ2 reservoirs. Results from the inverse modelling showed that low values of acoustic impedance from 19,743 to 20,487 (ft/s)(g/cc) reflect hydrocarbon-bearing reservoirs while medium to high values shows brine and shale respectively, with brine zone ranging from 20,487 to 22,531 (ft/s)(g/cc) and shale above 22,531 (ft/s)(g/cc). Two lithofacies were identified from inversion analysis of Vp/Vs and Mu-Rho, namely: sand and shale with VpVs 1.95 values respectively. Mu-Rho > 12.29 (GPa)(g/cc) and <12.29 (GPa) (g/cc) represent sand and shale respectively. From 3D volume, it was observed that a high accumulation of hydrocarbon was observed to be saturated at the north to the eastern part of the field forming a meandering channel. Sands were mainly distributed around the northeastern to the southwestern part of the field, that tends to be away from Well 029. This was also validated by the volume of rigidity modulus (Mu-Rho) showing high values indicating sands fall within the northeastern part of the field.展开更多
To establish the relationship among reservoir characteristics and rock physical parameters,we construct the well-bore rock physical models firstly,considering the influence factors,such as mineral composition,shale co...To establish the relationship among reservoir characteristics and rock physical parameters,we construct the well-bore rock physical models firstly,considering the influence factors,such as mineral composition,shale content,porosity,fluid type and saturation.Then with analyzing the change rules of elastic parameters along with the above influence factors and the cross-plots among elastic parameters,the sensitive elastic parameters of tight sandstone reservoir are determined,and the rock physics template of sweet spot is constructed to guide pre-stack seismic inversion.The results show that velocity ratio and Poisson impedance are the most sensitive elastic parameters to indicate the lithologic and gas-bearing properties of sweet spot in tight sandstone reservoir.The high-quality sweet spot is characterized by the lower velocity ratio and Poisson impedance.Finally,the actual seismic data are selected to predict the sweet spots in tight sandstone gas reservoirs,so as to verify the validity of the rock physical simulation results.The significant consistency between the relative logging curves and inversion results in different wells implies that the utilization of well-bore rock physical simulation can guide the prediction of sweet spot in tight sandstone gas reservoirs.展开更多
The relationship between fluid density and saturation and sonic wave velocity of rock samples taken from the WXS Depression in the South China Sea was studied by an oil-water replacement experiment under simulated in-...The relationship between fluid density and saturation and sonic wave velocity of rock samples taken from the WXS Depression in the South China Sea was studied by an oil-water replacement experiment under simulated in-situ temperature and pressure conditions.Two kinds of low-density oils(0.691 and 0.749 g/cm^3) and two kinds of high-density oils(0.834 and 0.873 g/cm^3) were used to saturate the rock samples at different oil-saturation states,and the saturated P- and S-wave velocities were measured.Through Gassmann's equation,the theoretical P- and S-wave velocities were also calculated by the fluid replacement method.With the comparison of the measured values and the theoretical values, this study comes to the following conclusions.(1) With the increase of oil saturation and the decrease of water saturation,the P-wave velocity of rock samples saturated by low-density oil increases and the changing rule is in accord with the effective fluid theory;the P-wave velocity of rock samples saturated by high-density oil decreases and the changing rule goes against the theory.(2) With the increase of oil density(namely 0.691→0.749→0.834→0.873 g/cm^3) when oil saturation is unchanged,P-wave velocity increases gradually.(3) The S-wave velocity is always stable and is not affected by the change of oil density and saturation.The results can be used to constrain pre-stack seismic inversion,and the variation rule of sonic wave velocity is valuable for hydrocarbon identification in the study area.展开更多
During the process of rock compression up to failure, the microwave radiation of rocks was observed and studied by remote sensing. The microwave radiometers for wave lengths of 8 mm, 2 cm and 10 cm were used. The expe...During the process of rock compression up to failure, the microwave radiation of rocks was observed and studied by remote sensing. The microwave radiometers for wave lengths of 8 mm, 2 cm and 10 cm were used. The experiments show that the luminance temperature of microwave radiation of rocks increases with increasing stress during loading, and that the increase speeds up just before rock failure. This phenomenon can be considered as a new precursor of rock failure. The concept of remote sensing rock mechanics (or remote sensing rock physics) advanced on the experimental basis of infrared remote sensing observation during rock failure has been enriched by these new experimental results. Its contents have been extended from infrared remote sensing to microwave remote sensing, and the prospects of applying remote sensing rock mechanics in the study of earthquake and rock burst prediction have also been broadened.展开更多
The electrophysical property of saturated rocks is very important for reservoir identification and evaluation. In this paper, the lattice Boltzmann method (LBM) was used to study the electrophysical property of rock...The electrophysical property of saturated rocks is very important for reservoir identification and evaluation. In this paper, the lattice Boltzmann method (LBM) was used to study the electrophysical property of rock saturated with fluid because of its advantages over conventional numerical approaches in handling complex pore geometry and boundary conditions. The digital core model was constructed through the accumulation of matrix grains based on their radius distribution obtained by the measurements of core samples. The flow of electrical current through the core model saturated with oil and water was simulated on the mesoscopic scale to reveal the non-Archie relationship between resistivity index and water saturation (I-Sw). The results from LBM simulation and laboratory measurements demonstrated that the I-Sw relation in the range of low water saturation was generally not a straight line in the log-log coordinates as described by the Archie equation. We thus developed a new equation based on numerical simulation and physical experiments. This new equation was used to fit the data from laboratory core measurements and previously published data. Determination of fluid saturation and reservoir evaluation could be significantly improved by using the new equation.展开更多
Prediction of seismic attenuation and dispersion that are inherently sensitive to hydraulic and elastic properties of the medium of interest in the presence of mesoscopic fractures and pores,is of great interest in th...Prediction of seismic attenuation and dispersion that are inherently sensitive to hydraulic and elastic properties of the medium of interest in the presence of mesoscopic fractures and pores,is of great interest in the characterization of fractured formations.This has been very difficult,however,considering that stress interactions between fractures and pores,related to their spatial distributions,tend to play a crucial role on affecting overall dynamic elastic properties that are largely unexplored.We thus choose to quantitatively investigate frequency-dependent P-wave characteristics in fractured porous rocks at the scale of a representative sample using a numerical scale-up procedure via performing finite element modelling.Based on 2-D numerical quasi-static experiments,effects of fracture and fluid properties on energy dissipation in response to wave-induced fluid flow at the mesoscopic scale are quantified via solving Biot's equations of consolidation.We show that numerical results are sensitive to some key characteristics of probed synthetic rocks containing unconnected and connected fractures,demonstrating that connectivity,aperture and inclination of fractures as well as fracture infills exhibit strong impacts on the two manifestations of WIFF mechanisms in the connected scenario,and on resulting total wave attenuation and phase velocity.This,in turn,illustrates the importance of these two WIFF mechanisms in fractured rocks and thus,a deeper understanding of them may eventually allow for a better characterization of fracture systems using seismic methods.Moreover,this presented work combines rock physics predictions with seismic numerical simulations in frequency domain to illustrate the sensitivity of seismic signatures on the monitoring of an idealized geologic CO_(2) sequestration in fractured reservoirs.The simulation demonstrates that these two WIFF mechanisms can strongly modify seismic records and hence,indicating that incorporating the two energy dissipation mechanisms in the geophysical interpretation can potentially improving the monitoring and surveying of fluid variations in fractured formations.展开更多
The identification of hydrocarbons using seismic methods is critical in the prediction of shale oil res-ervoirs.However,delineating shales of high oil saturation is challenging owing to the similarity in the elastic p...The identification of hydrocarbons using seismic methods is critical in the prediction of shale oil res-ervoirs.However,delineating shales of high oil saturation is challenging owing to the similarity in the elastic properties of oil-and water-bearing shales.The complexity of the organic matter properties associated with kerogen and hydrocarbon further complicates the characterization of shale oil reservoirs using seismic methods.Nevertheless,the inelastic shale properties associated with oil saturation can enable the utilization of velocity dispersion for hydrocarbon identification in shales.In this study,a seismic inversion scheme based on the fluid dispersion attribute was proposed for the estimation of hydrocarbon enrichment.In the proposed approach,the conventional frequency-dependent inversion scheme was extended by incorporating the PP-wave reflection coefficient presented in terms of the effective fluid bulk modulus.A rock physics model for shale oil reservoirs was constructed to describe the relationship between hydrocarbon saturation and shale inelasticity.According to the modeling results,the hydrocarbon sensitivity of the frequency-dependent effective fluid bulk modulus is superior to the traditional compressional wave velocity dispersion of shales.Quantitative analysis of the inversion re-sults based on synthetics also reveals that the proposed approach identifies the oil saturation and related hydrocarbon enrichment better than the above-mentioned conventional approach.Meanwhile,in real data applications,actual drilling results validate the superiority of the proposed fluid dispersion attribute as a useful hydrocarbon indicator in shale oil reservoirs.展开更多
The major storage space types in the carbonate reservoir in the Ordovician in the TZ45 area are secondary dissolution caves.For the prediction of caved carbonate reservoir,post-stack methods are commonly used in the o...The major storage space types in the carbonate reservoir in the Ordovician in the TZ45 area are secondary dissolution caves.For the prediction of caved carbonate reservoir,post-stack methods are commonly used in the oilfield at present since pre-stack inversion is always limited by poor seismic data quality and insufficient logging data.In this paper,based on amplitude preserved seismic data processing and rock-physics analysis,pre-stack inversion is employed to predict the caved carbonate reservoir in TZ45 area by seriously controlling the quality of inversion procedures.These procedures mainly include angle-gather conversion,partial stack,wavelet estimation,low-frequency model building and inversion residual analysis.The amplitude-preserved data processing method can achieve high quality data based on the principle that they are very consistent with the synthetics.Besides,the foundation of pre-stack inversion and reservoir prediction criterion can be established by the connection between reservoir property and seismic reflection through rock-physics analysis.Finally,the inversion result is consistent with drilling wells in most cases.It is concluded that integrated with amplitude-preserved processing and rock-physics,pre-stack inversion can be effectively applied in the caved carbonate reservoir prediction.展开更多
基金The authors acknowledge the sponsorship of National Natural Science Foundation of China(42174139,41974119,42030103)Laoshan Laboratory Science and Technology Innovation Program(LSKj202203406)Science Foundation from Innovation and Technology Support Program for Young Scientists in Colleges of Shandong Province and Ministry of Science and Technology of China(2019RA2136).
文摘Total organic carbon (TOC) prediction with elastic parameter inversions has been widely used in the identification and evaluation of source rocks. However, the elastic parameters used to predict TOC are not only determined by TOC but also depend on the other physical properties of source rocks. Besides, the TOC prediction with the elastic parameters inversion is an indirect method based on the statistical relationship obtained from well logs and experiment data. Therefore, we propose a rock physics model and define a TOC indicator mainly affected by TOC to predict TOC directly. The proposed rock physics model makes the equivalent elastic moduli of source rocks parameterized by the TOC indicator. Combining the equivalent elastic moduli of source rocks and Gray’s approximation leads to a novel linearized approximation of the P-wave reflection coefficient incorporating the TOC indicator. Model examples illustrate that the novel reflectivity approximation well agrees with the exact Zoeppritz equation until incident angles reach 40°. Convoluting the novel P-wave reflection approximation with seismic wavelets as the forward solver, an AVO inversion method based on the Bayesian theory is proposed to invert the TOC indicator with seismic data. The synthetic examples and field tests validate the feasibility and stability of the proposed AVO inversion approach. Using the inversion results of the TOC indicator, TOC is directly and accurately estimated in the target area.
基金supported by the National Natural Science Foundation of China(No.42174146)CNPC major forwardlooking basic science and technology projects(No.2021DJ0204).
文摘Rock physics inversion is to use seismic elastic properties of underground strata for predicting reservoir petrophysical parameters.The Markov chain Monte Carlo(MCMC)algorithm is commonly used to solve rock physics inverse problems.However,all the parameters to be inverted are iterated simultaneously in the conventional MCMC algorithm.What is obtained is an optimal solution of combining the petrophysical parameters with being inverted.This study introduces the alternating direction(AD)method into the MCMC algorithm(i.e.the optimized MCMC algorithm)to ensure that each petrophysical parameter can get the optimal solution and improve the convergence of the inversion.Firstly,the Gassmann equations and Xu-White model are used to model shaly sandstone,and the theoretical relationship between seismic elastic properties and reservoir petrophysical parameters is established.Then,in the framework of Bayesian theory,the optimized MCMC algorithm is used to generate a Markov chain to obtain the optimal solution of each physical parameter to be inverted and obtain the maximum posterior density of the physical parameter.The proposed method is applied to actual logging and seismic data and the results show that the method can obtain more accurate porosity,saturation,and clay volume.
基金supported by the National Natural Science Foundation of China under Grants U1262208the National Natural Science Foundation of China under Grants 41404090
文摘Kerogen content and kerogen porosity play a significant role in elastic properties of organic-rich shales. We construct a rock physics model for organic-rich shales to quantify the effect of kerogen content and kerogen porosity using the Kuster and Toks6z theory and the selfconsistent approximation method. Rock physics modeling results show that with the increase of kerogen content and kerogen-related porosity, the velocity and density of shales decrease, and the effect of kerogen porosity becomes more obvious only for higher kerogen content. We also find that the Poisson's ratio of the shale is not sensitive to kerogen porosity for the case of gas saturation. Finally, for the seismic reflection responses of an organic-rich shale layer, forward modeling results indicate the fifth type AVO re- sponses which correspond to a negative intercept and a positive gradient. The absolute values of intercept and gradient increase with kerogen content and kerogen porosity, and present predictable variations associated with velocities and density.
文摘Due to the huge differences between the unconventional shale and conventional sand reservoirs in many aspects such as the types and the characteristics of minerals,matrix pores and fluids,the construction of shale rock physics model is significant for the exploration and development of shale reservoirs.To make a better characterization of shale gas-bearing reservoirs,we first propose a new but more suitable rock physics model to characterize the reservoirs.We then use a well A to demonstrate the feasibility and reliability of the proposed rock physics model of shale gas-bearing reservoirs.Moreover,we propose a new brittleness indicator for the high-porosity and organic-rich shale gas-bearing reservoirs.Based on the parameter analysis using the constructed rock physics model,we finally compare the new brittleness indicator with the commonly used Young’s modulus in the content of quartz and organic matter,the matrix porosity,and the types of filled fluids.We also propose a new shale brittleness index by integrating the proposed new brittleness indicator and the Poisson’s ratio.Tests on real data sets demonstrate that the new brittleness indicator and index are more sensitive than the commonly used Young’s modulus and brittleness index for the high-porosity and high-brittleness shale gas-bearing reservoirs.
基金Science and Technology Cooperation Project of the CNPC-SWPU Innovation Alliance(No.2020CX010501)National Science and Technology Major ProjectNational Natural Science Foundation of China Petrochemical Joint Fund Project(U1762107)
文摘Digital rock physics(DRP)is a paramount technology to improve the economic benefits of oil and gas fields,devise more scientific oil and gas field development plans,and create digital oil and gas fields.Currently,a significant gap is present between DRP theory and practical applications.Conventional digital-core construction focuses only on simple cores,and the recognition and segmentation effect of fractures and pores of complex cores is poor.The identification of rock minerals is inaccurate,which leads to the difference between the digital and actual cores.To promote the application of DRP in developing oil and gas fields,based on the high-precision X-ray computed tomography scanning technology,the U-Net deep learning model of the full convolution neural network is used to segment the pores,fractures,and matrix from the complex rock core with natural fractures innovatively.Simultaneously,the distribution of rock minerals is divided,and the distribution of rock conditions is corrected by X-ray diffraction.A pore—fracture network model is established based on the equivalent radius,which lays the foundation for fluid seepage simulation.Finally,the accuracy of the established a digital core is verified by the porosity measured via nuclear magnetic resonance technology,which is of great significance to the development and application of DRP in oil and gas fields.
基金sponsored by the CNOOC research project (No.YXKY-2019-ZY-04)。
文摘Tight sandstone gas reservoirs have characteristics of low porosity and low permeability,complex pore structure.In the paper,we consider the pore and micro-fracture systems in tight sandstone reservoir and perform rock physics modeling of sandstone gas reservoirs based on Chapman.Then we analyze the effects of relaxation time,and difference fluid types on dispersion and attention,and seismic responses.We found that the magnitude of the relaxation time is related to the gas saturation and the viscosity coefficient of the fluid.The magnitude of the relaxation time affects the frequency range of the velocity dispersion and attenuation,and the type of fluid affects the magnitude of the velocity dispersion and attenuation.As the relaxation time decreases,the calculated amplitude deceases.And the amplitude and the waveform of the synthetic seismic record are related to different fluid types and formation thickness.
基金supported by NSFC(41930425)Science Foundation of China University of Petroleum,Beijing(No.2462020YXZZ008)+1 种基金R&D Department of China National Petroleum Corporation(Investigations on fundamental experiments and advanced theoretical methods in geophysical prospecting applications(2022DQ0604-01)the Strategic Cooperation Technology Projects of CNPC and CUPB(ZLZX2020-03)and NSFC(42274142).
文摘Heavy oil has high density and viscosity, and exhibits viscoelasticity. Gassmann's theory is not suitable for materials saturated with viscoelastic fluids. Directly applying such model leads to unreliable results for seismic inversion of heavy oil reservoir. To describe the viscoelastic behavior of heavy oil, we modeled the elastic properties of heavy oil with varying viscosity and frequency using the Cole-Cole-Maxwell (CCM) model. Then, we used a CCoherent Potential Approximation (CPA) instead of the Gassmann equations to account for the fluid effect, by extending the single-phase fluid condition to two-phase fluid (heavy oil and water) condition, so that partial saturation of heavy oil can be considered. This rock physics model establishes the relationship between the elastic modulus of reservoir rock and viscosity, frequency and saturation. The viscosity of the heavy oil and the elastic moduli and porosity of typical reservoir rock samples were measured in laboratory, which were used for calibration of the rock physics model. The well-calibrated frequency-variant CPA model was applied to the prediction of the P- and S-wave velocities in the seismic frequency range (1–100 Hz) and the inversion of petrophysical parameters for a heavy oil reservoir. The pre-stack inversion results of elastic parameters are improved compared with those results using the CPA model in the sonic logging frequency (∼10 kHz), or conventional rock physics model such as the Xu-Payne model. In addition, the inversion of the porosity of the reservoir was conducted with the simulated annealing method, and the result fits reasonably well with the logging curve and depicts the location of the heavy oil reservoir on the time slice. The application of the laboratory-calibrated CPA model provides better results with the velocity dispersion correction, suggesting the important role of accurate frequency dependent rock physics models in the seismic prediction of heavy oil reservoirs.
基金supported by the National Natural Science Foundation of China(Grant numbers 42274160 and 42074153).
文摘Seismic characterizing of tight gas sandstone (TGS) reservoirs is essential for identifying promising gas-bearing regions. However, exploring the petrophysical significance of seismic-inverted elastic properties is challenging due to the complex microstructures in TGSs. Meanwhile, interbedded structures of sandstone and mudstone intensify the difficulty in accurately extracting the crucial tight sandstone properties. An integrated rock-physics-based framework is proposed to estimate the reservoir quality of TGSs from seismic data. TGSs with complex pore structures are modeled using the double-porosity model, providing a practical tool to compute rock physics templates for reservoir parameter estimation. The VP/VS ratio is utilized to predict the cumulative thickness of the TGS reservoirs within the target range via the threshold value evaluated from wireline logs for lithology discrimination. This approach also facilitates better capturing the elastic properties of the TGSs for quantitative seismic interpretation. Total porosity is estimated from P-wave impedance using the correlation obtained based on wireline log analysis. After that, the three-dimensional rock-physics templates integrated with the estimated total porosity are constructed to interpret microfracture porosity and gas saturation from velocity ratio and bulk modulus. The integrated framework can optimally estimate the parameters dominating the reservoir quality. The results of the indicator proposed based on the obtained parameters are in good agreement with the gas productions and can be utilized to predict promising TGS reservoirs. Moreover, the results suggest that considering microfracture porosity allows a more accurate prediction of high-quality reservoirs, further validating the applicability of the proposed method in the studied region.
基金The authors would like to acknowledge financial support from NSFC Basic Research Program on Deep Petroleum Resource Accumulation and Key Engineering Technologies(U19B6003-04-03)National Natural Science Foundation of China(41930425)+2 种基金Beijing Natural Science Foundation(8222073),R&D Department of China National Petroleum Corporation(Investigations on fundamental experiments and advanced theoretical methods in geophysical prospecting applications,2022DQ0604-01)Scientific Research and Technology Development Project of PetroChina(2021DJ1206)National Key Research and Development Program of China(2018YFA0702504).
文摘Knowledge about the seismic elastic modulus dispersion,and associated attenuation,in fluid-saturated rocks is essential for better interpretation of seismic observations taken as part of hydrocarbon identification and time-lapse seismic surveillance of both conventional and unconventional reservoir and overburden performances.A Seismic Elastic Moduli Module has been developed,based on the forced-oscillations method,to experimentally investigate the frequency dependence of Young's modulus and Poisson's ratio,as well as the inferred attenuation,of cylindrical samples under different confining pressure conditions.Calibration with three standard samples showed that the measured elastic moduli were consistent with the published data,indicating that the new apparatus can operate reliably over a wide frequency range of f∈[1-2000,10^(6)]Hz.The Young's modulus and Poisson's ratio of the shale and the tight sandstone samples were measured under axial stress oscillations to assess the frequency-and pressure-dependent effects.Under dry condition,both samples appear to be nearly frequency independent,with weak pressure dependence for the shale and significant pressure dependence for the sandstone.In particular,it was found that the tight sandstone with complex pore microstructure exhibited apparent dispersion and attenuation under brine or glycerin saturation conditions,the levels of which were strongly influenced by the increased effective pressure.In addition,the measured Young's moduli results were compared with the theoretical predictions from a scaled poroelastic model with a reasonably good agreement,revealing that the combined fluid flow mechanisms at both mesoscopic and microscopic scales possibly responsible for the measured dispersion.
文摘A considerable effort has been made in the literature for quality assurance (QA) and quality checking (QC) of the petrophysical log data for computation of reservoir rock property parameters. Well log data plays an integral role in building a rock physics model for quantitative interpretation (QI) work. A poor-quality rock physics model may lead to significant financial and HSSE implications by drilling wells in undesired locations. Historically, a variety of techniques have been used including histograms and cross plots for reviewing the feasibility of petrophysical logs for QI work. However, no attempt has ever been made to introduce a simplified workflow. This paper serves two-fold. It provides a simplified step by step approach for building a petrophysics/rock physics model. A case study has been presented to compare the synthetic seismogram generated from the simplified workflow with the actual seismic trace at well locations. Secondly, the paper shows how a few key cross plots and rock property parameters provide adequate information to validate petrophysical data, distinguish overburden and reservoir sections, and to help identify fluids and saturation trends within the reservoir sands. In the mentioned case study, the robustness of the simplified rock physics model has helped seismic data to successfully distinguish hydrocarbon bearing reservoir sands from non-reservoir shales.
基金supported by the National Grand Project for Science and Technology(Grant Nos.2016ZX05024-004,2017ZX05009-001,2017ZX05036-005)the Science Foundation from SINOPEC Key Laboratory of Geophysics(Grant No.WTYJY-WX2019-0104)。
文摘Seismic Rock physics plays a bridge role between the rock moduli and physical properties of the hydrocarbon reservoirs.Prestack seismic inversion is an important method for the quantitative characterization of elasticity,physical properties,lithology and fluid properties of subsurface reservoirs.In this paper,a high order approximation of rock physics model for clastic rocks is established and one seismic AVO reflection equation characterized by the high order approximation(Jacobian and Hessian matrix)of rock moduli is derived.Besides,the contribution of porosity,shale content and fluid saturation to AVO reflectivity is analyzed.The feasibility of the proposed AVO equation is discussed in the direct estimation of rock physical properties.On the basis of this,one probabilistic AVO inversion based on differential evolution-Markov chain Monte Carlo stochastic model is proposed on the premise that the model parameters obey Gaussian mixture probability prior model.The stochastic model has both the global optimization characteristics of the differential evolution algorithm and the uncertainty analysis ability of Markov chain Monte Carlo model.Through the cross parallel of multiple Markov chains,multiple stochastic solutions of the model parameters can be obtained simultaneously,and the posterior probability density distribution of the model parameters can be simulated effectively.The posterior mean is treated as the optimal solution of the model to be inverted.Besides,the variance and confidence interval are utilized to evaluate the uncertainties of the estimated results,so as to realize the simultaneous estimation of reservoir elasticity,physical properties,discrete lithofacies and dry rock skeleton.The validity of the proposed approach is verified by theoretical tests and one real application case in eastern China.
文摘This paper reviews the recent progress in the studies of experimental rock mechanics and tectonophysics concerning seismology and physics of the Earths interior in China. The progress is presented in the following aspects: a) A lot of results of experiment and numerical simulation enrich our knowledge of the brittle fracturing process under the condition with heterogeneity in material and structure; b) Some new results on frictional behavior of non-homogeneous faults reveal the complexity of faulting behavior; c) Some new results on the brittle-plastic transition and plastic flow are obtained; especially the important progress is obtained on rheological properties of rocks in the lower crust and the upper mantle; d) A lot of experimental results are obtained on rock physics at high temperature and pressure and have been used in study of material composition and state. These results provide useful information for understanding the physical properties and deformation mechanisms of material of the Earths interior and earthquake physics.
文摘Simultaneous waveform inversion was used to predict lithofacies and fluid type across the field. Very often, characterizing reservoirs in terms of lithology and fluid type using conventional methods is replete with uncertainties, especially in marginal fields. An approach is employed in this study that integrated rock physics and waveform inverse modelling for lithology and fluid-type characterization to appropriately identify potential hydrocarbon saturated zones and their corresponding lithology. Seismic and well-log data were analyzed using Hampson Russel software. The method adopted includes lithofacies and fluid content analysis using rock physics parameters and seismic simultaneous inverse modelling. Rock physics analysis identified 2 broad reservoirs namely: HDZ1 and HDZ2 reservoirs. Results from the inverse modelling showed that low values of acoustic impedance from 19,743 to 20,487 (ft/s)(g/cc) reflect hydrocarbon-bearing reservoirs while medium to high values shows brine and shale respectively, with brine zone ranging from 20,487 to 22,531 (ft/s)(g/cc) and shale above 22,531 (ft/s)(g/cc). Two lithofacies were identified from inversion analysis of Vp/Vs and Mu-Rho, namely: sand and shale with VpVs 1.95 values respectively. Mu-Rho > 12.29 (GPa)(g/cc) and <12.29 (GPa) (g/cc) represent sand and shale respectively. From 3D volume, it was observed that a high accumulation of hydrocarbon was observed to be saturated at the north to the eastern part of the field forming a meandering channel. Sands were mainly distributed around the northeastern to the southwestern part of the field, that tends to be away from Well 029. This was also validated by the volume of rigidity modulus (Mu-Rho) showing high values indicating sands fall within the northeastern part of the field.
基金supported by the National Key R&D Program of China(Grant No.2018YFC1405900)the Major Projects of National Science and Technology(Grant Nos.2016ZX05011-002,2016ZX05027-002-005)+3 种基金the National Natural Science Foundation of China(Grant No.41806073)the Natural Science Foundation of Shandong Province(Grant No.ZR2017BD014)Shandong Provincial Key Laboratory of Depositional Mineralization and Sedimentary Minerals,Shandong University of Science and Technology(Grant No.DMSM2017042)the Fundamental Research Funds for the Central Universities(Grant No.201964016)
文摘To establish the relationship among reservoir characteristics and rock physical parameters,we construct the well-bore rock physical models firstly,considering the influence factors,such as mineral composition,shale content,porosity,fluid type and saturation.Then with analyzing the change rules of elastic parameters along with the above influence factors and the cross-plots among elastic parameters,the sensitive elastic parameters of tight sandstone reservoir are determined,and the rock physics template of sweet spot is constructed to guide pre-stack seismic inversion.The results show that velocity ratio and Poisson impedance are the most sensitive elastic parameters to indicate the lithologic and gas-bearing properties of sweet spot in tight sandstone reservoir.The high-quality sweet spot is characterized by the lower velocity ratio and Poisson impedance.Finally,the actual seismic data are selected to predict the sweet spots in tight sandstone gas reservoirs,so as to verify the validity of the rock physical simulation results.The significant consistency between the relative logging curves and inversion results in different wells implies that the utilization of well-bore rock physical simulation can guide the prediction of sweet spot in tight sandstone gas reservoirs.
基金supported by the Science and Technology Research Key Project of Ministry of Education (Grant No.109035)the National Natural Science Foundation Key Project(Grant No.40830423)the Key Project of Students Extra-curricular Science and Technology Research Program of Schlumberger(Grant No.SLBX0908)
文摘The relationship between fluid density and saturation and sonic wave velocity of rock samples taken from the WXS Depression in the South China Sea was studied by an oil-water replacement experiment under simulated in-situ temperature and pressure conditions.Two kinds of low-density oils(0.691 and 0.749 g/cm^3) and two kinds of high-density oils(0.834 and 0.873 g/cm^3) were used to saturate the rock samples at different oil-saturation states,and the saturated P- and S-wave velocities were measured.Through Gassmann's equation,the theoretical P- and S-wave velocities were also calculated by the fluid replacement method.With the comparison of the measured values and the theoretical values, this study comes to the following conclusions.(1) With the increase of oil saturation and the decrease of water saturation,the P-wave velocity of rock samples saturated by low-density oil increases and the changing rule is in accord with the effective fluid theory;the P-wave velocity of rock samples saturated by high-density oil decreases and the changing rule goes against the theory.(2) With the increase of oil density(namely 0.691→0.749→0.834→0.873 g/cm^3) when oil saturation is unchanged,P-wave velocity increases gradually.(3) The S-wave velocity is always stable and is not affected by the change of oil density and saturation.The results can be used to constrain pre-stack seismic inversion,and the variation rule of sonic wave velocity is valuable for hydrocarbon identification in the study area.
文摘During the process of rock compression up to failure, the microwave radiation of rocks was observed and studied by remote sensing. The microwave radiometers for wave lengths of 8 mm, 2 cm and 10 cm were used. The experiments show that the luminance temperature of microwave radiation of rocks increases with increasing stress during loading, and that the increase speeds up just before rock failure. This phenomenon can be considered as a new precursor of rock failure. The concept of remote sensing rock mechanics (or remote sensing rock physics) advanced on the experimental basis of infrared remote sensing observation during rock failure has been enriched by these new experimental results. Its contents have been extended from infrared remote sensing to microwave remote sensing, and the prospects of applying remote sensing rock mechanics in the study of earthquake and rock burst prediction have also been broadened.
基金sponsored by the project No.50404001 from the National Natural Science Foundation of Chinathe National Key Fundamental Research & Development Project(Grant No.2007CB209601)+1 种基金the China National PetroleumCorporation Fundamental Research Program (Grant No.06A30102)the China Postdoctoral Science Foundation(Project No.2004035350)
文摘The electrophysical property of saturated rocks is very important for reservoir identification and evaluation. In this paper, the lattice Boltzmann method (LBM) was used to study the electrophysical property of rock saturated with fluid because of its advantages over conventional numerical approaches in handling complex pore geometry and boundary conditions. The digital core model was constructed through the accumulation of matrix grains based on their radius distribution obtained by the measurements of core samples. The flow of electrical current through the core model saturated with oil and water was simulated on the mesoscopic scale to reveal the non-Archie relationship between resistivity index and water saturation (I-Sw). The results from LBM simulation and laboratory measurements demonstrated that the I-Sw relation in the range of low water saturation was generally not a straight line in the log-log coordinates as described by the Archie equation. We thus developed a new equation based on numerical simulation and physical experiments. This new equation was used to fit the data from laboratory core measurements and previously published data. Determination of fluid saturation and reservoir evaluation could be significantly improved by using the new equation.
文摘Prediction of seismic attenuation and dispersion that are inherently sensitive to hydraulic and elastic properties of the medium of interest in the presence of mesoscopic fractures and pores,is of great interest in the characterization of fractured formations.This has been very difficult,however,considering that stress interactions between fractures and pores,related to their spatial distributions,tend to play a crucial role on affecting overall dynamic elastic properties that are largely unexplored.We thus choose to quantitatively investigate frequency-dependent P-wave characteristics in fractured porous rocks at the scale of a representative sample using a numerical scale-up procedure via performing finite element modelling.Based on 2-D numerical quasi-static experiments,effects of fracture and fluid properties on energy dissipation in response to wave-induced fluid flow at the mesoscopic scale are quantified via solving Biot's equations of consolidation.We show that numerical results are sensitive to some key characteristics of probed synthetic rocks containing unconnected and connected fractures,demonstrating that connectivity,aperture and inclination of fractures as well as fracture infills exhibit strong impacts on the two manifestations of WIFF mechanisms in the connected scenario,and on resulting total wave attenuation and phase velocity.This,in turn,illustrates the importance of these two WIFF mechanisms in fractured rocks and thus,a deeper understanding of them may eventually allow for a better characterization of fracture systems using seismic methods.Moreover,this presented work combines rock physics predictions with seismic numerical simulations in frequency domain to illustrate the sensitivity of seismic signatures on the monitoring of an idealized geologic CO_(2) sequestration in fractured reservoirs.The simulation demonstrates that these two WIFF mechanisms can strongly modify seismic records and hence,indicating that incorporating the two energy dissipation mechanisms in the geophysical interpretation can potentially improving the monitoring and surveying of fluid variations in fractured formations.
基金supported by the National Natural Science Foundation of China(Grant numbers 42074153 and 42274160)the Open Research Fund of SINOPEC Key Laboratory of Geophysics(Grant number 33550006-20-ZC0699-0006).
文摘The identification of hydrocarbons using seismic methods is critical in the prediction of shale oil res-ervoirs.However,delineating shales of high oil saturation is challenging owing to the similarity in the elastic properties of oil-and water-bearing shales.The complexity of the organic matter properties associated with kerogen and hydrocarbon further complicates the characterization of shale oil reservoirs using seismic methods.Nevertheless,the inelastic shale properties associated with oil saturation can enable the utilization of velocity dispersion for hydrocarbon identification in shales.In this study,a seismic inversion scheme based on the fluid dispersion attribute was proposed for the estimation of hydrocarbon enrichment.In the proposed approach,the conventional frequency-dependent inversion scheme was extended by incorporating the PP-wave reflection coefficient presented in terms of the effective fluid bulk modulus.A rock physics model for shale oil reservoirs was constructed to describe the relationship between hydrocarbon saturation and shale inelasticity.According to the modeling results,the hydrocarbon sensitivity of the frequency-dependent effective fluid bulk modulus is superior to the traditional compressional wave velocity dispersion of shales.Quantitative analysis of the inversion re-sults based on synthetics also reveals that the proposed approach identifies the oil saturation and related hydrocarbon enrichment better than the above-mentioned conventional approach.Meanwhile,in real data applications,actual drilling results validate the superiority of the proposed fluid dispersion attribute as a useful hydrocarbon indicator in shale oil reservoirs.
基金supported by National Basic Research Program(2006CB202304)of Chinaco-supported by the National Basic Research Program of China(Grant No.2011CB201103)the National Science and Technology Major Project of China(Grant No.2011ZX05004003)
文摘The major storage space types in the carbonate reservoir in the Ordovician in the TZ45 area are secondary dissolution caves.For the prediction of caved carbonate reservoir,post-stack methods are commonly used in the oilfield at present since pre-stack inversion is always limited by poor seismic data quality and insufficient logging data.In this paper,based on amplitude preserved seismic data processing and rock-physics analysis,pre-stack inversion is employed to predict the caved carbonate reservoir in TZ45 area by seriously controlling the quality of inversion procedures.These procedures mainly include angle-gather conversion,partial stack,wavelet estimation,low-frequency model building and inversion residual analysis.The amplitude-preserved data processing method can achieve high quality data based on the principle that they are very consistent with the synthetics.Besides,the foundation of pre-stack inversion and reservoir prediction criterion can be established by the connection between reservoir property and seismic reflection through rock-physics analysis.Finally,the inversion result is consistent with drilling wells in most cases.It is concluded that integrated with amplitude-preserved processing and rock-physics,pre-stack inversion can be effectively applied in the caved carbonate reservoir prediction.