(Multichannel)Singular spectrum analysis is considered as one of the most effective methods for seismic incoherent noise suppression.It utilizes the low-rank feature of seismic signal and regards the noise suppression...(Multichannel)Singular spectrum analysis is considered as one of the most effective methods for seismic incoherent noise suppression.It utilizes the low-rank feature of seismic signal and regards the noise suppression as a low-rank reconstruction problem.However,in some cases the seismic geophones receive some erratic disturbances and the amplitudes are dramatically larger than other receivers.The presence of this kind of noise,called erratic noise,makes singular spectrum analysis(SSA)reconstruction unstable and has undesirable effects on the final results.We robustify the low-rank reconstruction of seismic data by a reweighted damped SSA(RD-SSA)method.It incorporates the damped SSA,an improved version of SSA,into a reweighted framework.The damping operator is used to weaken the artificial disturbance introduced by the low-rank projection of both erratic and random noise.The central idea of the RD-SSA method is to iteratively approximate the observed data with the quadratic norm for the first iteration and the Tukeys bisquare norm for the rest iterations.The RD-SSA method can suppress seismic incoherent noise and keep the reconstruction process robust to the erratic disturbance.The feasibility of RD-SSA is validated via both synthetic and field data examples.展开更多
3D eikonal equation is a partial differential equation for the calculation of first-arrival traveltimes and has been widely applied in many scopes such as ray tracing,source localization,reflection migration,seismic m...3D eikonal equation is a partial differential equation for the calculation of first-arrival traveltimes and has been widely applied in many scopes such as ray tracing,source localization,reflection migration,seismic monitoring and tomographic imaging.In recent years,many advanced methods have been developed to solve the 3D eikonal equation in heterogeneous media.However,there are still challenges for the stable and accurate calculation of first-arrival traveltimes in 3D strongly inhomogeneous media.In this paper,we propose an adaptive finite-difference(AFD)method to numerically solve the 3D eikonal equation.The novel method makes full use of the advantages of different local operators characterizing different seismic wave types to calculate factors and traveltimes,and then the most accurate factor and traveltime are adaptively selected for the convergent updating based on the Fermat principle.Combined with global fast sweeping describing seismic waves propagating along eight directions in 3D media,our novel method can achieve the robust calculation of first-arrival traveltimes with high precision at grid points either near source point or far away from source point even in a velocity model with large and sharp contrasts.Several numerical examples show the good performance of the AFD method,which will be beneficial to many scientific applications.展开更多
Compared with the transverse isotropic(TI)medium,the orthorhombic anisotropic medium has both horizontal and vertical symmetry axes and it can be approximated as a set of vertical fissures developed in a group of hori...Compared with the transverse isotropic(TI)medium,the orthorhombic anisotropic medium has both horizontal and vertical symmetry axes and it can be approximated as a set of vertical fissures developed in a group of horizontal strata.Although the full-elastic orthorhombic anisotropic wave equation can accurately simulate seismic wave propagation in the underground media,a huge computational cost is required in seismic modeling,migration,and inversion.The conventional coupled pseudo-acoustic wave equations based on acoustic approximation can be used to significantly reduce the cost of calculation.However,these equations usually suffer from unwanted shear wave artifacts during wave propagation,and the presence of these artifacts can significantly degrade the imaging quality.To solve these problems,we derived a new pure P-wave equation for orthorhombic media that eliminates shear wave artifacts while compromising computational efficiency and accuracy.In addition,the derived equation involves pseudo-differential operators and it must be solved by 3D FFT algorithms.In order to reduce the number of 3D FFT,we utilized the finite difference and pseudo-spectral methods to conduct 3D forward modeling.Furthermore,we simplified the equation by using elliptic approximation and implemented 3D reverse-time migration(RTM).Forward modeling tests on several homogeneous and heterogeneous models confirm that the accuracy of the new equation is better than that of conventional methods.3D RTM imaging tests on three-layer and SEG/EAGE 3D salt models confirm that the ORT media have better imaging quality.展开更多
Geo-microbial prospecting for hydrocarbons is an exploration method based on the seepage of light gaseous hydrocarbons from oil/gas reservoirs to the surface and their utilization by hydrocarbon oxidizing bacteria. Th...Geo-microbial prospecting for hydrocarbons is an exploration method based on the seepage of light gaseous hydrocarbons from oil/gas reservoirs to the surface and their utilization by hydrocarbon oxidizing bacteria. These bacteria utilize hydrocarbon gases as their sole source of food and are found to be enriched in the near surface soils/sediments above the oil and gas reservoirs. The detection of anomalous populations of n-pentane and n-hexane oxidizing bacteria in the surface soils can help to evaluate the prospects for hydrocarbon exploration. A geo-microbial survey has been carried out in the Bikaner Nagaur basin to investigate the prospects for hydrocarbon exploration. In the present study, bacterial counts for n-pentane utilizing bacteria range between 2.0× 10^2 and 1.26× 10^6 cfu/gm and n-hexane utilizing bacteria range between 2.0× 10^2 and 1.21 × 10^6 cfu/gm. The bacterial concentration distribution maps show four distinct anomalies in the study area. The possibility of discovering oil or gas reservoirs using the microbiological method is emphasized by the tact that the hydrocarbon-oxidizing bacteria range between 103 and 106 cfu/gm in soil/sediment receiving hydrocarbon micro-seepages. In the present study area of the Bikaner Nagaur basin, n-pentane and n-hexane utilizing bacteria are found between 105 and 106 cfu/gm of soil sample, which is significant and thereby substantiates the seepage of lighter hydrocarbon accumulations from oil and gas reservoirs. Geo-microbial prospecting studies suggest that hydrocarbon micro-seepage of subsurface origin is present in the study area and indicate that the area has positive prospects for petroleum exploration.展开更多
Full-waveform inversion(FWI)is a powerful tool to reconstruct subsurface geophysical parameters with high resolution.As3 D surveys become widely implemented,corresponding 3 D processing techniques are required to solv...Full-waveform inversion(FWI)is a powerful tool to reconstruct subsurface geophysical parameters with high resolution.As3 D surveys become widely implemented,corresponding 3 D processing techniques are required to solve complex geological cases,while a large amount of computation is the most challenging problem.We propose an adaptive variable-grid 3 D FWI on graphics processing unit devices to improve computational efficiency without losing accuracy.The irregular-grid discretization strategy is based on a dispersion relation,and the grid size adapts to depth,velocity,and frequency automatically.According to the transformed grid coordinates,we derive a modified acoustic wave equation and apply it to full wavefield simulation.The 3 D variable-grid modeling is conducted on several 3 D models to validate its feasibility,accuracy and efficiency.Then we apply the proposed modeling method to full-waveform inversion for source and residual wavefield propagation.It is demonstrated that the adaptive variable-grid FWI is capable of decreasing computing time and memory requirements.From the inversion results of the 3 D SEG/EAGE overthrust model,our method retains inversion accuracy when recovering both thrust and channels.展开更多
Envelope inversion(El)is an efficient tool to mitigate the nonlinearity of conventional full waveform inversion(FWI)by utilizing the ultralow-frequency component in the seismic data.However,the performance of envelope...Envelope inversion(El)is an efficient tool to mitigate the nonlinearity of conventional full waveform inversion(FWI)by utilizing the ultralow-frequency component in the seismic data.However,the performance of envelope inversion depends on the frequency component and initial model to some extent.To improve the convergence ability and avoid the local minima issue,we propose a convolution-based envelope inversion method to update the low-wavenumber component of the velocity model.Besides,the multi-scale inversion strategy(MCEI)is also incorporated to improve the inversion accuracy while guaranteeing the global convergence.The success of this method relies on modifying the original envelope data to expand the overlap region between observed and modeled envelope data,which in turn expands the global minimum basin of misfit function.The accurate low-wavenumber component of the velocity model provided by MCEI can be used as the migration model or an initial model for conventional FWI.The numerical tests on simple layer model and complex BP 2004 model verify that the proposed method is more robust than El even when the initial model is coarse and the frequency component of data is high.展开更多
Using field geological survey,drilling and seismic data,combined with the study of regional tectonic evolution and structural deformation,as well as lithological and sedimentary analysis,we reconstructed the basin fil...Using field geological survey,drilling and seismic data,combined with the study of regional tectonic evolution and structural deformation,as well as lithological and sedimentary analysis,we reconstructed the basin filling process and paleo-geography of north Tarim Basin in Early Cambrian,aiming to analyze the factors controlling the distribution and spatial architecture of the subsalt reservoir and source units and to define the favorable exploration direction.The Late Sinian tectonic activities in the northern Tarim Basin were characterized by different patterns in different areas,which controlled the sedimentary pattern in the Early Cambrian.The boundary faults of Nanhuaian rift basin in the south slope of Tabei uplift and the north slope of Tazhong uplift became reactivated in the Early Cambrian,forming two NEE and EW striking subsidence centers and depocenters,where the predicted thickness of the Yurtusi Formation could reach 250 meters.In the Xiaoerbulake period,the weak rimmed platform was developed in the hanging wall of syndepositional fault.Whereas the Nanhuaian rift system in the Tadong and Manxi areas were uplifted and destroyed in the Late Sinian,and appeared as gently slope transiting toward the subsidence center in the Early Cambrian.The former had the sedimentary features of hybrid facies platform and the latter had the sedimentary features of ramp platform.The black shale of the Yurtus Formation in the footwall of syndepositional fault and the reef bank of Xiaoerbulake Formation platform margin in the hanging wall in Early Cambrian constitute a predicable source-reservoir combination.The activity intensity of syndepositional fault controlled the thickness of black shale and the scale of the reef bank.It is suggested carrying out high accuracy seismic exploration to determine the location of Early Cambrian syndepositional faults,on this basis,to search the reef bank of Xiaoerbulake Formation along the faults westward,and then drill risk exploration wells at sites where traps are shallow in buried depth.展开更多
Through the analysis of logging,field outcrops,cores and geochemical data,and based on the study of the relationships between sea level changes,sequence filling,paleo-geomorphy and lithofacies,the sequence lithofacies...Through the analysis of logging,field outcrops,cores and geochemical data,and based on the study of the relationships between sea level changes,sequence filling,paleo-geomorphy and lithofacies,the sequence lithofacies paleo-geography and evolution process of the Lower Permian Liangshan-Qixia Formation(Qixia Stage for short)in Sichuan Basin and its surrounding areas are restored.The Qixia Stage can be divided into three third-order sequences,in which SQ0,SQ1 and SQ2 are developed in the depression area,and SQ1 and SQ2 are only developed in other areas.The paleo-geomorphy reflected by the thickness of each sequence indicates that before the deposition of the Qixia Stage in the Early Permian,the areas surrounding the Sichuan Basin are characterized by“four uplifts and four depressions”,namely,four paleo-uplifts/paleo-lands of Kangdian,Hannan,Shennongjia and Xuefeng Mountain,and four depressions of Chengdu-Mianyang,Kangdian front,Jiangkou and Yichang;while the interior of the basin is characterized by“secondary uplifts,secondary depressions and alternating convex-concave”.SQ2 is the main shoal forming period of the Qixia Formation,and the high-energy mound shoal facies mainly developed in the highs of sedimentary paleo-geomorphy and the relative slope break zones.The distribution of dolomitic reservoirs(dolomite,limy dolomite and dolomitic limestone)has a good correlation with the sedimentary geomorphic highs and slope break zones.The favorable mound-shoal and dolomitic reservoirs are distributed around depressions at platform-margin and along highs and around sags in the basin.It is pointed out that the platform-margin area in western Sichuan Basin is still the key area for exploration at present;while areas around Chengdu-Mianyang depression and Guangwang secondary depression inside the platform and areas around sags in central Sichuan-southern Sichuan are favorable exploration areas for dolomitic reservoirs of the Qixia Formation in the next step.展开更多
The long-term seepage of hydrocarbons, either as macroseepage or microseepage, can set up near-surface oxidation reduction zones that favor the development of a diverse array of chemical and mineralogical changes. The...The long-term seepage of hydrocarbons, either as macroseepage or microseepage, can set up near-surface oxidation reduction zones that favor the development of a diverse array of chemical and mineralogical changes. The bacterial oxidation of light hydrocarbons can directly or indirectly bring about significant changes in the values of pH and Eh of the surrounding environment, thereby also changing the stability fields of the different mineral species present in that environment. The paper reports the role of hydrocarbon microseepage in surface alterations of trace metal concentrations. In this study trace metal alterations were mapped that appear to be associated with hydrocarbon microseepages in the oil/ gas fields. A total of 50 soil samples were collected near oil and gas fields of the Tatipaka and Pasarlapudi areas of the Krishna Godavari Basin, Andhra Pradesh. The soil samples were collected from a depth of 2-2.5 m. The paper reports the chemical alterations associated with trace metals in soils that are related to hydrocarbon microseepages above some of the major oil and gas fields of this petroliferous region. Trace metals, such as scandium （Sc）, vanadium （V）, chromium （Cr）, cobalt （Co）, nickel （Ni）, copper （Cu）, zinc （Zn）, barium （Ba） and strontium （Sr）, in soil samples were analyzed using inductively coupled plasma- mass spectrometry （ICP-MS）. The concentrations of Sc （8 to 40 mg/kg）, V （197 to 489 mg/kg）, Cr （106 to 287 mg/kg）, Co （31 to 52 mg/kg）, Ni （65 to 110 mg/kg）, Cu （88 to 131 mg/kg）, Zn （88 to 471 mg/kg）, Ba （263 to 3,091 mg/kg） and Sr （119 to 218 mg/kg） were obtained. It was observed that the concentrations of trace elements were tremendously increased when they were compared with their normal concentrations in soils. The analysis of adsorbed soil gas showed the presence of high concentrations of YC2＋ （C2H6, C3H8 and n-C4H10） ranging from 7 to 222 μg/kg respectively. Integrated studies of trace elements over adsorbed light gaseous hydrocarbons （2C2＋） anomalies showed good correlation with the existing oil and gas wells. The carbon isotopic composition of δ13C1 of the samples ranges between -36.6%o to -22.7‰ （Pee Dee Belemnite） values indicate thermogenic origin, which presents convincing evidence that the adsorbed soil gases collected from these sediments are of catagenetic origin. The increase in the concentrations of trace metals near oil/gas producing areas, suggests a soil chemical change to a reducing environment, presumably due to the influence of hydrocarbon microseepage, which could be applied with other geoscientific data to identify areas of future hydrocarbon exploration in frontier areas.展开更多
The Tongnan secondary negative structure in central Sichuan Basin has controls and influences on the structural framework and petroleum geological conditions in the Gaoshiti-Moxi area.To clarify the controls and influ...The Tongnan secondary negative structure in central Sichuan Basin has controls and influences on the structural framework and petroleum geological conditions in the Gaoshiti-Moxi area.To clarify the controls and influences,the deformation characteristics,structural attributes and evolution process of the Tongnan negative structure were investigated through a series of qualitative and quantitative methods such as balanced profile restoration,area-depth-strain(ADS)analysis,and structural geometric forward numerical simulation,after comprehensive structural interpretation of high-precision 3D seismic data.The results are obtained in three aspects.First,above and below the P/AnP(Permian/pre-Permian)unconformity,the Tongnan negative structure demonstrates vertical differential structural deformation.It experiences two stages of structural stacking and reworking:extensional depression(from the Sinian Dengying Formation to the Permian),and compressional syncline deformation(after the Jurassic).The multi-phase trishear deformation of the preexisting deep normal faults dominated the extensional depression.The primary depression episodes occurred in the periods from the end of Late Proterozoic to the deposition of the 1st–2nd members of the Dengying Formation,and from the deposition of Lower Cambrian Longwangmiao Formation–Middle–Upper Cambrian until the Ordovician.Second,the multi-stage evolution process of the Tongnan negative structure controlled the oil and gas migration and adjustment and present-day differential gas and water distribution between the Tongnan negative structure and the Gaoshiti and Moxi-Longnüsi structural highs.Third,the Ordovician,which is limitedly distributed in the Tongnan negative structure and is truncated by the P/AnP unconformity on the top,has basic geological conditions for the formation of weathering karst carbonate reservoirs.It is a new petroleum target deserving attention.展开更多
Metalloporphyrins have devastating effects on the deep processing and efficient utilization of petroleum resources. Various porphyrins are derived from porphin by substitution at the β, meso, or both positions. Herei...Metalloporphyrins have devastating effects on the deep processing and efficient utilization of petroleum resources. Various porphyrins are derived from porphin by substitution at the β, meso, or both positions. Herein, we focus on the interaction between nickel(Ⅱ) and the porphyrin skeleton. Five different density functionals were selected for determining a reliable approach to simulate porphyrin-nickel(Ⅱ) complexes. At D4h symmetry, the 10b1g, 15a1g, 4e1g and 9b2g were dominated by the 3d character of nickel. Compared with pure density functional, the hybrid functional gave larger highest occupied molecular orbital(HOMO) and lowest unoccupied molecular orbital(LUMO) gap, and the higher HOMO-LUMO gap also corresponded to a higher excitation energy. The characteristic Soret and Q bands modeled by timedependent density functional theory(TDDFT) calculations matched well with experimental results, and the absorption was discussed in detail with natural transition orbitals(NTO) analysis. Furthermore, the binding character of nickel(Ⅱ) with various substituted porphyrins were calculated at perdew-burke-ernzerhof(PBE) level, showed that the geometry and complexation behavior of porphyrin-nickel(Ⅱ) complexes could be significantly tuned by different substituents. The symmetry of the complex was reduced and the skeleton ring get distortion when introducing different substituents. Compared with ETIO porphyrin, the introduction of phenyl on the β position can expand the conjugated system and promote the covalent character of Ni-N bond.展开更多
Because of the difference ofoil and gas accumulation condition between the hanging wall and the footwall of a fault, there is a peculiar accumulation mechanism that oil and gas mainly exists in the hanging wall of the...Because of the difference ofoil and gas accumulation condition between the hanging wall and the footwall of a fault, there is a peculiar accumulation mechanism that oil and gas mainly exists in the hanging wall of the basement fault, but in the footwall of the shallow detachment fault in the Nanbaxian pool. The oil and gas of the Nanbaxian pool came from the mature Jurassic hydrocarbon source rock of the Yibei depression located at the south of the Nanbaxian pool. Firstly, the oil and gas accumulated in the traps of the hanging wall of the basement fault by way of the unconformity and the basement faults, and turned into some primary deep pools; and then, the shallow detachment fault that formed in the later tectonic movement broke into the deep primary pools, which caused the oil and gas migration upwards along the basement faults and the shallow detachment faults and the evolvement into some secondary oil and gas pools later. The history of the Nanbaxian oil and gas accumulation can be summarized successively as the syndepositional upheaval controlled by faults; single hydrocarbon source rock; unconformities and faults as migration channels; buoyancy, overpressure and tectonic stress as dynamic forces; multistage migration and accumulation of oil and gas; and finally an overlapped double-floor pattern of oil and gas accumulation. The most important explorative targets in the north of the Qaidam Basin are traps connected with the primary pools in the footwall by shallow detachment faults.展开更多
Based on comprehensive analysis of cores,thin sections,logging and seismic data,the sequence stratigraphy and sedimentary evolution of the third and fourth members of Sinian Dengying Formation(Deng 3 and Deng 4 member...Based on comprehensive analysis of cores,thin sections,logging and seismic data,the sequence stratigraphy and sedimentary evolution of the third and fourth members of Sinian Dengying Formation(Deng 3 and Deng 4 members for short)in the Gaomo area of Sichuan Basin were investigated,and the favorable zones for reservoir development in the Deng 3 Member and Deng 4 Member were predicted.(1)Two Type I and one Type II sequence boundaries are identified in the Deng 3 and Deng 4 members.Based on the identified sequence boundaries,the Deng 3 and Deng 4 members can be divided into two third order sequences SQ3 and SQ4,which are well-developed,isochronal and traceable in this area;the SQ3 thins from west to the east,and the SQ4 thins from northwest to southeast.(2)The sedimentary environment from the depositional period of SQ3 to SQ4 has experienced the evolution from mixed platform to rimmed platform,and the platform rimmed system on the west side is characterized by the development of platform margin microbial mound and grain shoal assemblages.The intraplatform area is a restricted platform facies composed of a variety of dolomites,and there are local micro-geomorphic highlands of different scales and scattered intraplatform mounds and shoals.(3)The Deng 4 Member reservoirs,with obvious facies-controlled characteristic,are mainly distributed in the upper part of high-frequency upward shallow cycle and the high-stand systems tract of the third-order sequence vertically,and are more developed in the platform margin belt than in the intraplatform belt,and more developed in the Gaoshiti platform margin belt than in the west Suining platform margin belt on the plane.(4)Three types of favorable reservoir zones of Deng 4 Member have been finely delineated with 3D seismic data;among them,the mound and shoal facies zones developed in the ancient highlands of the intraplatform are the first choice for the next exploration and development of the Deng 4 Member.展开更多
Surface-related multiples frequently propagate into the subsurface and contain abundant information on small reflection angles.Compared with the conventional migration of primaries,migration of multiples offers comple...Surface-related multiples frequently propagate into the subsurface and contain abundant information on small reflection angles.Compared with the conventional migration of primaries,migration of multiples offers complementary illumination and a higher vertical resolution.However,crosstalk artifacts caused by unrelated multiples during reverse time migration(RTM)using multiples severely degrade the reliability and interpretation of the final migration images.Therefore,we proposed RTM using first-order receiver-side water-bottom-related multiples for eliminating crosstalk artifacts and enhancing vertical resolution.We first backward propagate the first-order receiver-side water-bottom-related multiples using a water-layer model,followed by saving the upper boundary wavefield.Then we produce the source wavefield using a seismic wavelet and the receiver wavefield by back-extrapolating the saved boundary.Finally,the cross-correlation imaging condition is applied to generate the final image.This method transforms the receiver-side multiples into primaries,followed by the conventional migration processing procedures.Numerical examples using synthetic datasets demonstrate that our method significantly enhances the imaging quality by eliminating crosstalk artifacts and improving the resolution.展开更多
The Gaussian beam migration(GBM) is a steady imaging approach, which has high accuracy and efficiency. Its implementation mainly includes the traditional frequency domain and the recent popular space-time domain. Firs...The Gaussian beam migration(GBM) is a steady imaging approach, which has high accuracy and efficiency. Its implementation mainly includes the traditional frequency domain and the recent popular space-time domain. Firstly, we use the upward ray tracing strategy to get the backward wavefields. Then,we use the dominant frequency of the seismic data to simplify the imaginary traveltime calculation of the wavefields, which can cut down the Fourier transform number compared with the traditional GBM in the space-time domain. In addition, we choose an optimized parameter for the take-off angle increment of the up-going and down-going rays. These optimizations help us get an efficient space-time-domain acoustic GBM approach. Typical four examples show that the proposed method can significantly improve the computational efficiency up to one or even two orders of magnitude in different models with different model parameters and produce good imaging results with comparable accuracy and resolution with the traditional GBM in the space-time domain.展开更多
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.展开更多
Acoustic wave velocity has been commonly utilized to predict subsurface geopressure using empirical relations.Acoustic wave velocity is, however, affected by many factors. To estimate pore pressure accurately, we here...Acoustic wave velocity has been commonly utilized to predict subsurface geopressure using empirical relations.Acoustic wave velocity is, however, affected by many factors. To estimate pore pressure accurately, we here propose to use elastic rock physics models to understand and analyze quantitatively the various contributions from these different factors affecting wave velocity. We report a closed-form relationship between the frame flexibility factor(γ) in a rock physics model and differential pressure, which presents the major control of pressure on elastic properties such as bulk modulus and compressional wave velocity. For a gas-bearing shale with abundant micro-cracks and fractures, its bulk modulus is much lower at abnormally high pore pressure(high γ values) where thin cracks and flat pores are open than that at normal hydrostatic pressure(low γ values) where pores are more rounded on average. The developed relations between bulk modulus and differential pressure have been successfully applied to the Upper Ordovician Wufeng and Lower Silurian Longmaxi formations in the Dingshan area of the Sichuan Basin to map the three-dimensional spatial distribution of pore pressure in the shale, integrating core, log and seismic data. The estimated results agree well with field measurements. Pressure coefficient is positively correlated to gas content. The relations and methods reported here could be useful for hydrocarbon exploration, production, and drilling safety in both unconventional and conventional fields.展开更多
基金supported by the National Natural Science Foundation of China under grant no.42374133the Beijing Nova Program under grant no.2022056+1 种基金the Fundamental Research Funds for the Central Universities under grant no.2462020YXZZ006the Young Elite Scientists Sponsorship Program by CAST(YESS)under grant no.2018QNRC001。
文摘(Multichannel)Singular spectrum analysis is considered as one of the most effective methods for seismic incoherent noise suppression.It utilizes the low-rank feature of seismic signal and regards the noise suppression as a low-rank reconstruction problem.However,in some cases the seismic geophones receive some erratic disturbances and the amplitudes are dramatically larger than other receivers.The presence of this kind of noise,called erratic noise,makes singular spectrum analysis(SSA)reconstruction unstable and has undesirable effects on the final results.We robustify the low-rank reconstruction of seismic data by a reweighted damped SSA(RD-SSA)method.It incorporates the damped SSA,an improved version of SSA,into a reweighted framework.The damping operator is used to weaken the artificial disturbance introduced by the low-rank projection of both erratic and random noise.The central idea of the RD-SSA method is to iteratively approximate the observed data with the quadratic norm for the first iteration and the Tukeys bisquare norm for the rest iterations.The RD-SSA method can suppress seismic incoherent noise and keep the reconstruction process robust to the erratic disturbance.The feasibility of RD-SSA is validated via both synthetic and field data examples.
基金The authors thank the funds supported by the China National Nuclear Corporation under Grants Nos.WUQNYC2101 and WUHTLM2101-04National Natural Science Foundation of China(42074132,42274154).
文摘3D eikonal equation is a partial differential equation for the calculation of first-arrival traveltimes and has been widely applied in many scopes such as ray tracing,source localization,reflection migration,seismic monitoring and tomographic imaging.In recent years,many advanced methods have been developed to solve the 3D eikonal equation in heterogeneous media.However,there are still challenges for the stable and accurate calculation of first-arrival traveltimes in 3D strongly inhomogeneous media.In this paper,we propose an adaptive finite-difference(AFD)method to numerically solve the 3D eikonal equation.The novel method makes full use of the advantages of different local operators characterizing different seismic wave types to calculate factors and traveltimes,and then the most accurate factor and traveltime are adaptively selected for the convergent updating based on the Fermat principle.Combined with global fast sweeping describing seismic waves propagating along eight directions in 3D media,our novel method can achieve the robust calculation of first-arrival traveltimes with high precision at grid points either near source point or far away from source point even in a velocity model with large and sharp contrasts.Several numerical examples show the good performance of the AFD method,which will be beneficial to many scientific applications.
基金supported by the Marine S&T Fund of Shandong Province for Pilot National Laboratory for Marine Science and Technology(No.2021QNLM020001)the Major Scientific and Technological Projects of Shandong Energy Group(No.SNKJ2022A06-R23)+2 种基金the Funds of Creative Research Groups of China(No.41821002)National Natural Science Foundation of China Outstanding Youth Science Fund Project(Overseas)(No.ZX20230152)the Major Scientific and Technological Projects of CNPC(No.ZD2019-183-003)。
文摘Compared with the transverse isotropic(TI)medium,the orthorhombic anisotropic medium has both horizontal and vertical symmetry axes and it can be approximated as a set of vertical fissures developed in a group of horizontal strata.Although the full-elastic orthorhombic anisotropic wave equation can accurately simulate seismic wave propagation in the underground media,a huge computational cost is required in seismic modeling,migration,and inversion.The conventional coupled pseudo-acoustic wave equations based on acoustic approximation can be used to significantly reduce the cost of calculation.However,these equations usually suffer from unwanted shear wave artifacts during wave propagation,and the presence of these artifacts can significantly degrade the imaging quality.To solve these problems,we derived a new pure P-wave equation for orthorhombic media that eliminates shear wave artifacts while compromising computational efficiency and accuracy.In addition,the derived equation involves pseudo-differential operators and it must be solved by 3D FFT algorithms.In order to reduce the number of 3D FFT,we utilized the finite difference and pseudo-spectral methods to conduct 3D forward modeling.Furthermore,we simplified the equation by using elliptic approximation and implemented 3D reverse-time migration(RTM).Forward modeling tests on several homogeneous and heterogeneous models confirm that the accuracy of the new equation is better than that of conventional methods.3D RTM imaging tests on three-layer and SEG/EAGE 3D salt models confirm that the ORT media have better imaging quality.
文摘Geo-microbial prospecting for hydrocarbons is an exploration method based on the seepage of light gaseous hydrocarbons from oil/gas reservoirs to the surface and their utilization by hydrocarbon oxidizing bacteria. These bacteria utilize hydrocarbon gases as their sole source of food and are found to be enriched in the near surface soils/sediments above the oil and gas reservoirs. The detection of anomalous populations of n-pentane and n-hexane oxidizing bacteria in the surface soils can help to evaluate the prospects for hydrocarbon exploration. A geo-microbial survey has been carried out in the Bikaner Nagaur basin to investigate the prospects for hydrocarbon exploration. In the present study, bacterial counts for n-pentane utilizing bacteria range between 2.0× 10^2 and 1.26× 10^6 cfu/gm and n-hexane utilizing bacteria range between 2.0× 10^2 and 1.21 × 10^6 cfu/gm. The bacterial concentration distribution maps show four distinct anomalies in the study area. The possibility of discovering oil or gas reservoirs using the microbiological method is emphasized by the tact that the hydrocarbon-oxidizing bacteria range between 103 and 106 cfu/gm in soil/sediment receiving hydrocarbon micro-seepages. In the present study area of the Bikaner Nagaur basin, n-pentane and n-hexane utilizing bacteria are found between 105 and 106 cfu/gm of soil sample, which is significant and thereby substantiates the seepage of lighter hydrocarbon accumulations from oil and gas reservoirs. Geo-microbial prospecting studies suggest that hydrocarbon micro-seepage of subsurface origin is present in the study area and indicate that the area has positive prospects for petroleum exploration.
基金the SWPI group in China University of Petroleum(East China)for financial support and discussions
文摘Full-waveform inversion(FWI)is a powerful tool to reconstruct subsurface geophysical parameters with high resolution.As3 D surveys become widely implemented,corresponding 3 D processing techniques are required to solve complex geological cases,while a large amount of computation is the most challenging problem.We propose an adaptive variable-grid 3 D FWI on graphics processing unit devices to improve computational efficiency without losing accuracy.The irregular-grid discretization strategy is based on a dispersion relation,and the grid size adapts to depth,velocity,and frequency automatically.According to the transformed grid coordinates,we derive a modified acoustic wave equation and apply it to full wavefield simulation.The 3 D variable-grid modeling is conducted on several 3 D models to validate its feasibility,accuracy and efficiency.Then we apply the proposed modeling method to full-waveform inversion for source and residual wavefield propagation.It is demonstrated that the adaptive variable-grid FWI is capable of decreasing computing time and memory requirements.From the inversion results of the 3 D SEG/EAGE overthrust model,our method retains inversion accuracy when recovering both thrust and channels.
基金supported by the National Science Foundation(Grant No.41104069,41274124)National“973 Project”(Grant No.2014CB239006)+1 种基金National Oil and Gas Project(Grant No.2016ZX05014001,2016ZX05002)supported by Tai Shan Science Foundation for the Excellent Youth Scholars.
文摘Envelope inversion(El)is an efficient tool to mitigate the nonlinearity of conventional full waveform inversion(FWI)by utilizing the ultralow-frequency component in the seismic data.However,the performance of envelope inversion depends on the frequency component and initial model to some extent.To improve the convergence ability and avoid the local minima issue,we propose a convolution-based envelope inversion method to update the low-wavenumber component of the velocity model.Besides,the multi-scale inversion strategy(MCEI)is also incorporated to improve the inversion accuracy while guaranteeing the global convergence.The success of this method relies on modifying the original envelope data to expand the overlap region between observed and modeled envelope data,which in turn expands the global minimum basin of misfit function.The accurate low-wavenumber component of the velocity model provided by MCEI can be used as the migration model or an initial model for conventional FWI.The numerical tests on simple layer model and complex BP 2004 model verify that the proposed method is more robust than El even when the initial model is coarse and the frequency component of data is high.
基金Supported by the National Key Research and Development Program of China(2017YFC0603101)Strategic Priority Research Program of the Chinese Academy of Sciences(XDA14010101)Scientific Research and Technological Development Project of PetroChina(2018A-01).
文摘Using field geological survey,drilling and seismic data,combined with the study of regional tectonic evolution and structural deformation,as well as lithological and sedimentary analysis,we reconstructed the basin filling process and paleo-geography of north Tarim Basin in Early Cambrian,aiming to analyze the factors controlling the distribution and spatial architecture of the subsalt reservoir and source units and to define the favorable exploration direction.The Late Sinian tectonic activities in the northern Tarim Basin were characterized by different patterns in different areas,which controlled the sedimentary pattern in the Early Cambrian.The boundary faults of Nanhuaian rift basin in the south slope of Tabei uplift and the north slope of Tazhong uplift became reactivated in the Early Cambrian,forming two NEE and EW striking subsidence centers and depocenters,where the predicted thickness of the Yurtusi Formation could reach 250 meters.In the Xiaoerbulake period,the weak rimmed platform was developed in the hanging wall of syndepositional fault.Whereas the Nanhuaian rift system in the Tadong and Manxi areas were uplifted and destroyed in the Late Sinian,and appeared as gently slope transiting toward the subsidence center in the Early Cambrian.The former had the sedimentary features of hybrid facies platform and the latter had the sedimentary features of ramp platform.The black shale of the Yurtus Formation in the footwall of syndepositional fault and the reef bank of Xiaoerbulake Formation platform margin in the hanging wall in Early Cambrian constitute a predicable source-reservoir combination.The activity intensity of syndepositional fault controlled the thickness of black shale and the scale of the reef bank.It is suggested carrying out high accuracy seismic exploration to determine the location of Early Cambrian syndepositional faults,on this basis,to search the reef bank of Xiaoerbulake Formation along the faults westward,and then drill risk exploration wells at sites where traps are shallow in buried depth.
基金Supported by the PetroChina and Southwest Petroleum University Innovation Consortium Science and Technology Cooperation Project(2020CX010000)Basic Forward-Looking Project in Upstream Field of CNPC(2021DJ0501)General Program of NSFC(42172166).
文摘Through the analysis of logging,field outcrops,cores and geochemical data,and based on the study of the relationships between sea level changes,sequence filling,paleo-geomorphy and lithofacies,the sequence lithofacies paleo-geography and evolution process of the Lower Permian Liangshan-Qixia Formation(Qixia Stage for short)in Sichuan Basin and its surrounding areas are restored.The Qixia Stage can be divided into three third-order sequences,in which SQ0,SQ1 and SQ2 are developed in the depression area,and SQ1 and SQ2 are only developed in other areas.The paleo-geomorphy reflected by the thickness of each sequence indicates that before the deposition of the Qixia Stage in the Early Permian,the areas surrounding the Sichuan Basin are characterized by“four uplifts and four depressions”,namely,four paleo-uplifts/paleo-lands of Kangdian,Hannan,Shennongjia and Xuefeng Mountain,and four depressions of Chengdu-Mianyang,Kangdian front,Jiangkou and Yichang;while the interior of the basin is characterized by“secondary uplifts,secondary depressions and alternating convex-concave”.SQ2 is the main shoal forming period of the Qixia Formation,and the high-energy mound shoal facies mainly developed in the highs of sedimentary paleo-geomorphy and the relative slope break zones.The distribution of dolomitic reservoirs(dolomite,limy dolomite and dolomitic limestone)has a good correlation with the sedimentary geomorphic highs and slope break zones.The favorable mound-shoal and dolomitic reservoirs are distributed around depressions at platform-margin and along highs and around sags in the basin.It is pointed out that the platform-margin area in western Sichuan Basin is still the key area for exploration at present;while areas around Chengdu-Mianyang depression and Guangwang secondary depression inside the platform and areas around sags in central Sichuan-southern Sichuan are favorable exploration areas for dolomitic reservoirs of the Qixia Formation in the next step.
文摘The long-term seepage of hydrocarbons, either as macroseepage or microseepage, can set up near-surface oxidation reduction zones that favor the development of a diverse array of chemical and mineralogical changes. The bacterial oxidation of light hydrocarbons can directly or indirectly bring about significant changes in the values of pH and Eh of the surrounding environment, thereby also changing the stability fields of the different mineral species present in that environment. The paper reports the role of hydrocarbon microseepage in surface alterations of trace metal concentrations. In this study trace metal alterations were mapped that appear to be associated with hydrocarbon microseepages in the oil/ gas fields. A total of 50 soil samples were collected near oil and gas fields of the Tatipaka and Pasarlapudi areas of the Krishna Godavari Basin, Andhra Pradesh. The soil samples were collected from a depth of 2-2.5 m. The paper reports the chemical alterations associated with trace metals in soils that are related to hydrocarbon microseepages above some of the major oil and gas fields of this petroliferous region. Trace metals, such as scandium （Sc）, vanadium （V）, chromium （Cr）, cobalt （Co）, nickel （Ni）, copper （Cu）, zinc （Zn）, barium （Ba） and strontium （Sr）, in soil samples were analyzed using inductively coupled plasma- mass spectrometry （ICP-MS）. The concentrations of Sc （8 to 40 mg/kg）, V （197 to 489 mg/kg）, Cr （106 to 287 mg/kg）, Co （31 to 52 mg/kg）, Ni （65 to 110 mg/kg）, Cu （88 to 131 mg/kg）, Zn （88 to 471 mg/kg）, Ba （263 to 3,091 mg/kg） and Sr （119 to 218 mg/kg） were obtained. It was observed that the concentrations of trace elements were tremendously increased when they were compared with their normal concentrations in soils. The analysis of adsorbed soil gas showed the presence of high concentrations of YC2＋ （C2H6, C3H8 and n-C4H10） ranging from 7 to 222 μg/kg respectively. Integrated studies of trace elements over adsorbed light gaseous hydrocarbons （2C2＋） anomalies showed good correlation with the existing oil and gas wells. The carbon isotopic composition of δ13C1 of the samples ranges between -36.6%o to -22.7‰ （Pee Dee Belemnite） values indicate thermogenic origin, which presents convincing evidence that the adsorbed soil gases collected from these sediments are of catagenetic origin. The increase in the concentrations of trace metals near oil/gas producing areas, suggests a soil chemical change to a reducing environment, presumably due to the influence of hydrocarbon microseepage, which could be applied with other geoscientific data to identify areas of future hydrocarbon exploration in frontier areas.
基金Supported by the National Natural Science Foundation of China(U19B6003-01).
文摘The Tongnan secondary negative structure in central Sichuan Basin has controls and influences on the structural framework and petroleum geological conditions in the Gaoshiti-Moxi area.To clarify the controls and influences,the deformation characteristics,structural attributes and evolution process of the Tongnan negative structure were investigated through a series of qualitative and quantitative methods such as balanced profile restoration,area-depth-strain(ADS)analysis,and structural geometric forward numerical simulation,after comprehensive structural interpretation of high-precision 3D seismic data.The results are obtained in three aspects.First,above and below the P/AnP(Permian/pre-Permian)unconformity,the Tongnan negative structure demonstrates vertical differential structural deformation.It experiences two stages of structural stacking and reworking:extensional depression(from the Sinian Dengying Formation to the Permian),and compressional syncline deformation(after the Jurassic).The multi-phase trishear deformation of the preexisting deep normal faults dominated the extensional depression.The primary depression episodes occurred in the periods from the end of Late Proterozoic to the deposition of the 1st–2nd members of the Dengying Formation,and from the deposition of Lower Cambrian Longwangmiao Formation–Middle–Upper Cambrian until the Ordovician.Second,the multi-stage evolution process of the Tongnan negative structure controlled the oil and gas migration and adjustment and present-day differential gas and water distribution between the Tongnan negative structure and the Gaoshiti and Moxi-Longnüsi structural highs.Third,the Ordovician,which is limitedly distributed in the Tongnan negative structure and is truncated by the P/AnP unconformity on the top,has basic geological conditions for the formation of weathering karst carbonate reservoirs.It is a new petroleum target deserving attention.
基金supported by National Key R & D program of China (No.2021YFA1501204)Contract Projects of China Petroleum & Chemical Corporation (SINOPEC Corp.) (grant number 123017grant number 121043-2)。
文摘Metalloporphyrins have devastating effects on the deep processing and efficient utilization of petroleum resources. Various porphyrins are derived from porphin by substitution at the β, meso, or both positions. Herein, we focus on the interaction between nickel(Ⅱ) and the porphyrin skeleton. Five different density functionals were selected for determining a reliable approach to simulate porphyrin-nickel(Ⅱ) complexes. At D4h symmetry, the 10b1g, 15a1g, 4e1g and 9b2g were dominated by the 3d character of nickel. Compared with pure density functional, the hybrid functional gave larger highest occupied molecular orbital(HOMO) and lowest unoccupied molecular orbital(LUMO) gap, and the higher HOMO-LUMO gap also corresponded to a higher excitation energy. The characteristic Soret and Q bands modeled by timedependent density functional theory(TDDFT) calculations matched well with experimental results, and the absorption was discussed in detail with natural transition orbitals(NTO) analysis. Furthermore, the binding character of nickel(Ⅱ) with various substituted porphyrins were calculated at perdew-burke-ernzerhof(PBE) level, showed that the geometry and complexation behavior of porphyrin-nickel(Ⅱ) complexes could be significantly tuned by different substituents. The symmetry of the complex was reduced and the skeleton ring get distortion when introducing different substituents. Compared with ETIO porphyrin, the introduction of phenyl on the β position can expand the conjugated system and promote the covalent character of Ni-N bond.
文摘Because of the difference ofoil and gas accumulation condition between the hanging wall and the footwall of a fault, there is a peculiar accumulation mechanism that oil and gas mainly exists in the hanging wall of the basement fault, but in the footwall of the shallow detachment fault in the Nanbaxian pool. The oil and gas of the Nanbaxian pool came from the mature Jurassic hydrocarbon source rock of the Yibei depression located at the south of the Nanbaxian pool. Firstly, the oil and gas accumulated in the traps of the hanging wall of the basement fault by way of the unconformity and the basement faults, and turned into some primary deep pools; and then, the shallow detachment fault that formed in the later tectonic movement broke into the deep primary pools, which caused the oil and gas migration upwards along the basement faults and the shallow detachment faults and the evolvement into some secondary oil and gas pools later. The history of the Nanbaxian oil and gas accumulation can be summarized successively as the syndepositional upheaval controlled by faults; single hydrocarbon source rock; unconformities and faults as migration channels; buoyancy, overpressure and tectonic stress as dynamic forces; multistage migration and accumulation of oil and gas; and finally an overlapped double-floor pattern of oil and gas accumulation. The most important explorative targets in the north of the Qaidam Basin are traps connected with the primary pools in the footwall by shallow detachment faults.
基金Petro China-Southwest Petroleum University Innovation Consortium Science and Technology Cooperation Project(2020CX010000)。
文摘Based on comprehensive analysis of cores,thin sections,logging and seismic data,the sequence stratigraphy and sedimentary evolution of the third and fourth members of Sinian Dengying Formation(Deng 3 and Deng 4 members for short)in the Gaomo area of Sichuan Basin were investigated,and the favorable zones for reservoir development in the Deng 3 Member and Deng 4 Member were predicted.(1)Two Type I and one Type II sequence boundaries are identified in the Deng 3 and Deng 4 members.Based on the identified sequence boundaries,the Deng 3 and Deng 4 members can be divided into two third order sequences SQ3 and SQ4,which are well-developed,isochronal and traceable in this area;the SQ3 thins from west to the east,and the SQ4 thins from northwest to southeast.(2)The sedimentary environment from the depositional period of SQ3 to SQ4 has experienced the evolution from mixed platform to rimmed platform,and the platform rimmed system on the west side is characterized by the development of platform margin microbial mound and grain shoal assemblages.The intraplatform area is a restricted platform facies composed of a variety of dolomites,and there are local micro-geomorphic highlands of different scales and scattered intraplatform mounds and shoals.(3)The Deng 4 Member reservoirs,with obvious facies-controlled characteristic,are mainly distributed in the upper part of high-frequency upward shallow cycle and the high-stand systems tract of the third-order sequence vertically,and are more developed in the platform margin belt than in the intraplatform belt,and more developed in the Gaoshiti platform margin belt than in the west Suining platform margin belt on the plane.(4)Three types of favorable reservoir zones of Deng 4 Member have been finely delineated with 3D seismic data;among them,the mound and shoal facies zones developed in the ancient highlands of the intraplatform are the first choice for the next exploration and development of the Deng 4 Member.
基金partially funded by the National Natural Science Foundation of China(Grant No.41730425)the Special Fund of the Institute of Geophysics,China Earthquake Administration(Grant No.DQJB20K42)the Institute of Geology and Geophysics,Chinese Academy of Sciences Project(Grant No.IGGCAS-2019031)。
文摘Surface-related multiples frequently propagate into the subsurface and contain abundant information on small reflection angles.Compared with the conventional migration of primaries,migration of multiples offers complementary illumination and a higher vertical resolution.However,crosstalk artifacts caused by unrelated multiples during reverse time migration(RTM)using multiples severely degrade the reliability and interpretation of the final migration images.Therefore,we proposed RTM using first-order receiver-side water-bottom-related multiples for eliminating crosstalk artifacts and enhancing vertical resolution.We first backward propagate the first-order receiver-side water-bottom-related multiples using a water-layer model,followed by saving the upper boundary wavefield.Then we produce the source wavefield using a seismic wavelet and the receiver wavefield by back-extrapolating the saved boundary.Finally,the cross-correlation imaging condition is applied to generate the final image.This method transforms the receiver-side multiples into primaries,followed by the conventional migration processing procedures.Numerical examples using synthetic datasets demonstrate that our method significantly enhances the imaging quality by eliminating crosstalk artifacts and improving the resolution.
基金jointly supported by the National Key Research and Development Program of China (2019YFC0605503)the National Natural Science Foundation of China (41821002, 41922028,41874149)+1 种基金the Strategic Priority Research Program of the Chinese Academy of Sciences (XDA14010303)the Major Scientific and Technological Projects of CNPC (ZD2019-183-003)。
文摘The Gaussian beam migration(GBM) is a steady imaging approach, which has high accuracy and efficiency. Its implementation mainly includes the traditional frequency domain and the recent popular space-time domain. Firstly, we use the upward ray tracing strategy to get the backward wavefields. Then,we use the dominant frequency of the seismic data to simplify the imaginary traveltime calculation of the wavefields, which can cut down the Fourier transform number compared with the traditional GBM in the space-time domain. In addition, we choose an optimized parameter for the take-off angle increment of the up-going and down-going rays. These optimizations help us get an efficient space-time-domain acoustic GBM approach. Typical four examples show that the proposed method can significantly improve the computational efficiency up to one or even two orders of magnitude in different models with different model parameters and produce good imaging results with comparable accuracy and resolution with the traditional GBM in the space-time domain.
基金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.
文摘Acoustic wave velocity has been commonly utilized to predict subsurface geopressure using empirical relations.Acoustic wave velocity is, however, affected by many factors. To estimate pore pressure accurately, we here propose to use elastic rock physics models to understand and analyze quantitatively the various contributions from these different factors affecting wave velocity. We report a closed-form relationship between the frame flexibility factor(γ) in a rock physics model and differential pressure, which presents the major control of pressure on elastic properties such as bulk modulus and compressional wave velocity. For a gas-bearing shale with abundant micro-cracks and fractures, its bulk modulus is much lower at abnormally high pore pressure(high γ values) where thin cracks and flat pores are open than that at normal hydrostatic pressure(low γ values) where pores are more rounded on average. The developed relations between bulk modulus and differential pressure have been successfully applied to the Upper Ordovician Wufeng and Lower Silurian Longmaxi formations in the Dingshan area of the Sichuan Basin to map the three-dimensional spatial distribution of pore pressure in the shale, integrating core, log and seismic data. The estimated results agree well with field measurements. Pressure coefficient is positively correlated to gas content. The relations and methods reported here could be useful for hydrocarbon exploration, production, and drilling safety in both unconventional and conventional fields.
