中扬子海相碳酸盐岩中方解石脉成岩环境研究

通过对中扬子京山县城周边露头碳酸盐岩中方解石脉的普通薄片、阴极发光以及碳氧同位素等资料分析,认为该区露头三叠系裂隙中充填放射状和等轴粒状的方解石且阴极发光强度有差异,二叠系和奥陶系裂隙中主要为等轴粒状方解石,阴极发光与围岩相似。方解石脉的δ13CPDB位于-6.76‰～4.01‰之间,δ18 OPDB在-17.95‰～-5.67‰之间,基本为海水潜流带—混合水潜流带沉积环境。三叠系方解石脉基本上为沉积碳酸盐岩,二叠系和奥陶系部分方解石脉受后期成岩流体的影响,其中二叠系有机质成熟产生的流体对该层系方解石脉的形成有一定影响。

Effect of pore structure on seismic rock-physics characteristics of dense carbonates

The Ordovician carbonate rocks of the Yingshan formation in the Tarim Basin have a complex pore structure owing to diagenetic and secondary structures. Seismic elastic parameters（e.g., wave velocity） depend on porosity and pore structure. We estimated the average specific surface, average pore-throat radius, pore roundness, and average aspect ratio of carbonate rocks from the Tazhong area. High P-wave velocity samples have small average specific surface, small average pore-throat radius, and large average aspect ratio. Differences in the pore structure of dense carbonate samples lead to fluid-related velocity variability. However, the relation between velocity dispersion and average specific surface, or the average aspect ratio, is not linear. For large or small average specific surface, the pore structure of the rock samples becomes uniform, which weakens squirt fl ow and minimizes the residuals of ultrasonic data and predictions with the Gassmann equation. When rigid dissolved（casting mold） pores coexist with less rigid microcracks, there are significant P-wave velocity differences between measurements and predictions.

A rock-physical modeling method for carbonate reservoirs at seismic scale

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.

Research on the Seismic Wave Field of Karst Cavern Reservoirs near Deep Carbonate Weathered Crusts

Fracture and cavern hydrocarbon reservoirs in carbonates are an important pool type worldwide. The karst cavern reservoirs are easiest to identify on seismic reflection data. The prediction, exploration, and development of this type of reservoir require theoretical research on seismic wave fields reflected from complex inhomogeneous media. We compute synthetic seismic sections for fluidfilled cavern reservoirs of various heights and widths using random media models and inhomogeneous media elastic wave equations. Results indicate that even caverns significantly smaller than 1/ 4 wavelength are detectible on conventional band-width seismic sections as diffractions migrated into bead-type events. Diffraction amplitude is a function of cavern height and width. We introduce a width-amplitude factor which can be used to calculate the diffraction amplitude of a cavern with a limited width from the diffraction amplitude computed for an infinitely wide cavern.

BISQ model based on a Kelvin-Voigt viscoelastic frame in a partially saturated porous medium

Taking into account three important porous media mechanisms during wave propagation （the Biot-flow, squirt-flow, and solid-skeleton viscoelastic mechanisms）, we introduce water saturation into the dynamic governing equations of wave propagation by analyzing the effective medium theory and then providing a viscoelastic Biot/squirt （BISQ） model which can analyze the wave propagation problems in a partially viscous pore fluid saturated porous media. In this model, the effects of pore fluid distribution patterns on the effective bulk modulus at different frequencies are considered. Then we derive the wave dynamic equations in the time-space domain. The phase velocity and the attenuation coefficient equations of the viscoelatic BISQ model in the frequency-wavenumber domain are deduced through a set of plane harmonic solution assumptions. Finally, by means of numerical simulations, we investigate the effects of water saturation, permeability, and frequency on compressional wave velocity and attenuation. Based on tight sandstone and carbonate experimental observed data, the compressional wave velocities of partially saturated reservoir rocks are calculated. The compressional wave velocity in carbonate reservoirs is more sensitive to gas saturation than in sandstone reservoirs.

Application of Prediction Techniques in Carbonate Karst Reservoir in Tarim Basin

Carbonate karst reservoir is the emphases of Tarim＇s carbonate exploration. However, it is buried at a large depth, which results in Weak seismic reflection signal and low S/N ratio. In addition, the karst reservoir contains great heterogeneity, so reservoir prediction is very difficult. Through many years of research and exploration, we have established a suite of comprehensive evaluation technology for carbonate karst reservoir using geophysical characteristics and a geological concept model, including a technique for reconstructing the paleogeomorphology of buried hills based on a sequence framework, seismic description of the karst reservoir, and strain variant analysis for fracture estimation. The evaluation technology has been successfully applied in the Tabei and Tazhong areas, and commercial production of oil and gas has been achieved. We show the application of this technology in the Lunguxi area in North Tarim in this paper.

Calculation of saturation in carbonate reservoirs based on electrical efficiency

We derived an equation for saturation in carbonate reservoirs based on the electrical efficiency model in the case of lacking core data. Owing to the complex pore structure and strong heterogeneity in carbonate reservoirs, the relation between electrical efficiency and water porosity is either complex or linear. We proposed an electrical efficiency equation that accounts for the relation between electrical efficiency and water porosity. We also proposed a power-law relation between electrical efficiency and deep-formation resistivity and analyzed the factors controlling the error in the water saturation computations. We concluded that the calculation accuracy of the electrical efficiency is critical to the application of the saturation equation. The saturation equation was applied to the carbonate reservoirs of three wells in Iraq and Indonesia. For relative rock electrical efficiency error below 0.1, the water saturation absolute error is also below 0.1. Therefore, we infer that the proposed saturation equation generally satisfies the evaluation criteria for carbonate reservoirs.

Application of Comprehensive Geophysical Techniques to Predict Carbonate Fractured Reservoirs inQingxi Oilfield

Inter-crystalline pores, cavities, and fractures created from diagenetic shrinkage of dolomite are inter-connected each other, forming fine oil- and gas-bearing reservoirs. It is hard to predict these complex fracture-cavity reservoirs because of their random distribution, different growth timing, and so on. Taking the lacustrine dolomite fracture-pore reservoir in the Lower Cretaceous Xiagou Formation in the Qingxi oilfield within the Jiuquan basin as an example, we put forward a comprehensive geophysical method to predict carbonate fractures.

Anisotropic rock physics models for interpreting pore structures in carbonate reservoirs

We developed an anisotropic effective theoretical model for modeling the elastic behavior of anisotropic carbonate reservoirs by combining the anisotropic self-consistent approximation and differential effective medium models.By analyzing the measured data from carbonate samples in the TL area,a carbonate pore-structure model for estimating the elastic parameters of carbonate rocks is proposed,which is a prerequisite in the analysis of carbonate reservoirs.A workflow for determining elastic properties of carbonate reservoirs is established in terms of the anisotropic effective theoretical model and the pore-structure model.We performed numerical experiments and compared the theoretical prediction and measured data.The result of the comparison suggests that the proposed anisotropic effective theoretical model can account for the relation between velocity and porosity in carbonate reservoirs.The model forms the basis for developing new tools for predicting and evaluating the properties of carbonate reservoirs.

The Segmentation of FMI Image Based on 2-D Dyadic Wavelet Transform

A key aspect in extracting quantitative information from FMI logs is to segment the FMI image to get images of pores, vugs and fractures. A segmentation method based on the dyadic wavelet transform in 2-D is introduced in this paper. The first step is to find all the edge pixels of the FMI image using the 2-D wavelet transform. The second step is to calculate a segmentation threshold based on the average value of the edge pixels. Field data processing examples show that sub-images of vugs and fractures can be correctly separated from original FMI data continuously and automatically along the depth axis. The image segmentation lays the foundation for in-situ parameter calculation.

Seismic Prediction of Prolific Oil Zones in Carbonate Reservoirs with Extremely Low Porosity and Permeability under Salt

The Carboniferous reservoir in KJ oilfield is a carbonate reservoir with extremely low porosity and permeability and high-pressure. The reservoir has severe heterogeneity, is deeply buried, has complex master control factors, is covered with thick salt, all of which result in the serious distortion of reflection time and amplitudes under the salt, the poor seismic imaging, and the low S/N ratio and resolution. The key to developing this kind of reservoir is to correctly predict the distribution of highly profitable oil zones. In this paper we start by analyzing the master control factors, perform seismic-log calibration, optimize the seismic attributes indicating the lithofacies, karst, petrophysical properties, and fractures, and combine these results with the seismic, geology, log, oil reservoir engineering, and well data. We decompose the seismic prediction into six key areas： structural interpretation, prediction of lithofacies, karst, petrophysical properties, fractures, and then perform an integrated assessment. First, based on building the models of faults and fractures, sedimentary facies, and karst, we predict the distribution of the most favorable reservoir zones qualitatively. Then, using multi-parameter inversion and integrated multi-attribute analysis, we predict the favorable reservoir distribution quantitatively and semi-quantitatively to clarify the distribution of high-yield zones. We finally have a reliable basis for optimal selection of exploration and development targets.

Porous-grain–upper-boundary model and its application to Tarim Basin carbonates

Most of the carbonates in the Tarim Basin in northwest China are low-porosity and low-permeability rocks. Owing to the complexity of porosity in carbonates, conventional rock- physics models do not describe the relation between velocity and porosity for the Tarim Basin carbonates well. We propose the porous-grain-upper-boundary （PGU） model for estimating the relation between velocity and porosity for low-porosity carbonates. In this model, the carbonate sediments are treated as packed media of porous elastic grains, and the carbonate pores are divided into isolated and connected pores The PGU model is modified from the porous-grain-stiff-sand （PGST） model by replacing the critical porosity with the more practical isolated porosity. In the implementation, the effective elastic constants of the porous grains are calculated by using the differential effective medium （DEM） model. Then, the elastic constants of connected porous grains in dry rocks are calculated by using the modified upper Hashin-Shtrikman bound. The application to the Tarim carbonates shows that relative to other conventional effective medium models the PGU model matches the well log data well.

The constructing of pore structure factor in carbonate rocks and the inversion of reservoir parameters

With a more complex pore structure system compared with clastic rocks, carbonate rocks have not yet been well described by existing conventional rock physical models concerning the pore structure vagary as well as the influence on elastic rock properties. We start with a discussion and an analysis about carbonate rock pore structure utilizing rock slices. Then, given appropriate assumptions, we introduce a new approach to modeling carbonate rocks and construct a pore structure algorithm to identify pore structure mutation with a basis on the Gassmann equation and the Eshelby-Walsh ellipsoid inclusion crack theory. Finally, we compute a single well＇s porosity using this new approach with full wave log data and make a comparison with the predicted result of traditional method and simultaneously invert for reservoir parameters. The study results reveal that the rock pore structure can significantly influence the rocks＇ elastic properties and the predicted porosity error of the new modeling approach is merely 0.74%. Therefore, the approach we introduce can effectively decrease the predicted error of reservoir parameters.

The bound weighted average method(BWAM)for predicting S-wave velocity

The shear-wave velocity is a very important parameter in oil and gas seismic exploration, and vital in prestack elastic-parameters inversion and seismic attribute analysis. However, sheafing-velocity logging is seldom carried out because it is expensive. This paper presents a simple method for predicting S-wave velocity which covers the basic factors that influence seismic wave propagation velocity in rocks. The elastic modulus of a rock is expressed here as a weighted arithmetic average between Voigt and Reuss bounds, where the weighting factor, w, is a measurement of the geometric details of the pore space and mineral grains. The S-wave velocity can be estimated from w, which is derived from the P-wave modulus. The method is applied to process well-logging data for a carbonate reservoir in Sichuan Basin, and shows the predicted S-wave velocities agree well with the measured S-wave velocities.

Prediction of Fracture-Cavity System in Carbonate Reservoir： A Case Study in the Tahe Oilfield

The carbonate rocks in Tahe oilfield, which suffered from multi-period polycycle karstification and structure deformation, are heterogeneous reservoirs that are rich in pores, cavities,and fractures. The reservoirs are diversified in scale, space configuration, and complex in filling. For this kind of reservoir, a suite of seismic prestack or poststack prediction techniques has been developed based on the separation of seismic wave fields. Through cross-verification of the estimated results,a detailed description of palaeogeomorphology and structural features such as pores, cavities, and fractures in unaka has been achieved, the understanding of the spatial distribution of reservoir improved.

Pedogenic Carbonate and Soil Dehydrogenase Activity in Response to Soil Organic Matter in Artemisia ordosica Community

Little attention has been paid to the role of soil organic matter （OM） in the formation of pedogenic carbonate in desert soils. The relationships among soil OM, soil dehydrogenase activity （DHA）, and soil CaCO3 in a plant community dominated by Artemisia ordosica, located on the eastern boundary of Tcngger Desert in the Alxa League, Inner Mongolia, China, were studied to understand whether OM was directly involved in the formation of pedogenic carbonate. The results showed that DHA and CuCO3 positively correlated with OM content, and DHA, OM, and CaCO3 were correlated with each other in their spatial distribution, indicating that abundant OM content contributed to the formation of CaCO3. Therefore, the formation of pedogenic CaCO3 was a biotic process in the plant community dominated by A. ordosica.

Numerical simulation of the dual laterolog for carbonate cave reservoirs and response characteristics

Cave carbonate formations are characterized by heterogeneity, which makes electrical log prediction difficult. It is currently important to know how to use the dual laterolog to accurately identify and quantitatively evaluate caves. Using numerical simulation to calculate electrical log responses can provide a theoretical basis for cave identification and evaluation. In this paper, based on the dual laterolog principles, we first study different size spherical cave models using the finite element method （FEM）, determine a relation between resistivity and cave filling after comprehensively studying the log responses of cave models with different filling material, and finally study the dual laterolog responses on caves filled with shale, limestone, conglomerate, and thin laminated formation of sand and shale. The numerical results provide a theoretical basis for identification and evaluation of carbonate cave reservoirs.

Effects of micrite microtextures on the elastic and petrophysical properties of carbonate reservoirs

Apparent differences in sedimentation and diagenesis exist between carbonate reservoirs in different areas and affect their petrophysical and elastic properties.To elucidate the relevant mechanism,we study and analyze the characteristics of rock microstructure and elastic properties of carbonates and their variation regularity using 89 carbonate samples from the different areas The results show that the overall variation regularities of the physical and elastic properties of the carbonate rocks are controlled by the microtextures of the microcrystalline calcite,whereas the traditional classification of rock-and pore-structures is no longer applicable.The micrite microtextures can be divided,with respect to their morphological features,into porous micrite,compact micrite,and tight micrite.As the micrites evolves from the first to the last type,crystal boundaries are observed with increasingly close coalescence,the micritic intercrystalline porosity and pore-throat radius gradually decrease;meanwhile,the rigidity of the calcite microcrystalline particle boundary and elastic homogeneity are enhanced.As a result,the seismic elastic characteristics,such as permeability and velocity of samples,show a general trend of decreasing with the increase of porosity.For low-porosity rock samples(φ<5%)dominated by tight micrite,the micritic pores have limited contributions to porosity and permeability and the micrite elastic properties are similar to those of the rock matrix.In such cases,the macroscopic physical and elastic properties are more susceptible to the formation of cracks and dissolution pores,but these features are controlled by the pore structure.The pore aspect ratio can be used as a good indication of pore types.The bulk modulus aspect ratio for dissolution pores is greater than 0.2,whereas that of the intergranular pores ranges from 0.1 to 0.2.The porous and compact micrites are observed to have a bulk modulus aspect ratio less than 0.1,whereas the ratio of the tight micrite approaches 0.2。

Visco-elastic Properties of VES Diverting Acid for Carbonate Reservoirs

Storage modulus and loss modulus is the main performance index of visco-elastic properties.In this paper the storage modulus and loss modulus of a new diverting acid and their influencing factors were systematically investigated.Besides,the constitutive equations of the diverting acid at different temperatures were elicited from shearing experiments,which show that the visco-elastic surfactant(VES)acid system is a non-Newtonian power law fluid at low temperature and a Newtonian fluid at high temperature.The storage modulus and loss modulus at different temperatures,pH,and VES content in the acid are critical for the design of acid stimulation for oil well,especially when the VES acid is used in this field only on trial and the basic data are in urgent needed for the design and construction of the acidification stimulation.

Effects of carbon anhydrase on utilization of bicarbonate in microalgae:a case study in Lake Hongfeng

A bidirectional labeling method was established to distinguish the proportions of HCO3- and CO2 utiliza- tion pathways of microalgae in Lake Hongfeng. The method was based on microalgae cultured in a medium by adding equal concentrations of NaH13CO3 with different 613C values simultaneously. The inorganic carbon sources were quantified according to the stable carbon isotope composition in the treated microalgae. The effects of extracellular carbonic anhydrase （CAex） on the HCO3 and CO2 utilization pathways were distinguished using acetazolamide, a potent membrane-impermeable carbonic anhydrase inhibitor. The results show utilization of the added HCO3- was only 8% of the total carbon sources in karst lake. The proportion of the HCO3- utilization path- way was 52% of total inorganic carbon assimilation. Therefore, in the natural water of the karst area, the microalgae used less bicarbonate that preexisted in the aqueous medium than CO2 derived from the atmosphere. CAex increased the utilization of inorganic carbon from the atmosphere. The microalgae with CAex had greater carbon sequestration capacity in this karst area.