Acoustic reflection imaging in deep water wells is a new application scope for offshore hydrocarbon exploration.Two-dimensional(2 D)geological structure images can be obtained away from a one-dimensional(1 D)borehole ...Acoustic reflection imaging in deep water wells is a new application scope for offshore hydrocarbon exploration.Two-dimensional(2 D)geological structure images can be obtained away from a one-dimensional(1 D)borehole using single-well acoustic reflection imaging.Based on the directivity of dipole source and four-component dipole data,one can achieve the azimuth detection and the three-dimensional(3 D)structural information around the wellbore can be obtained.We first perform matrix rotation on the field fourcomponent data.Then,a series of processing steps are applied to the rotated dipole data to obtain the reflector image.According to the above dipole shear-wave imaging principle,we used four-component cross-dipole logging data from a deviated well in the South China Sea to image geological structures within 50 m of a deviated well,which can delineate the structural configuration and determine its orientation.The configuration of near-borehole bedding boundaries and fault structures from shear-wave imaging results agrees with those from the Inline and Xline seismic profiles of the study area.In addition,the configuration and orientation of the fault structure images are consistent with regional stress maps and the results of the borehole stress anisotropy analysis.Furthermore,the dip azimuth of the bedding boundary images was determined using borehole wall resistivity data.Results of this study indicate that integrating borehole acoustic reflection with seismic imaging not only fills the gap between the two measurement scales but also accurately delineates geological structures in the borehole vicinity.展开更多
The technology of simultaneous-source acquisition of seismic data excited by several sources can significantly improve the data collection efficiency. However, direct imaging of simultaneous-source data or blended dat...The technology of simultaneous-source acquisition of seismic data excited by several sources can significantly improve the data collection efficiency. However, direct imaging of simultaneous-source data or blended data may introduce crosstalk noise and affect the imaging quality. To address this problem, we introduce a structure-oriented filtering operator as preconditioner into the multisource least-squares reverse-time migration (LSRTM). The structure-oriented filtering operator is a nonstationary filter along structural trends that suppresses crosstalk noise while maintaining structural information. The proposed method uses the conjugate-gradient method to minimize the mismatch between predicted and observed data, while effectively attenuating the interference noise caused by exciting several sources simultaneously. Numerical experiments using synthetic data suggest that the proposed method can suppress the crosstalk noise and produce highly accurate images.展开更多
The distributions of local velocity and local phase holdup along the radial direction of pipes are complicated because of gravity differentiation,and the distribution of fluid velocity fi eld changes along the gravity...The distributions of local velocity and local phase holdup along the radial direction of pipes are complicated because of gravity differentiation,and the distribution of fluid velocity fi eld changes along the gravity direction in horizontal wells.Therefore,measuring the mixture flow and water holdup is difficult,resulting in poor interpretation accuracy of the production logging output profile.In this paper,oil–water two-phase flow dynamic simulation logging experiments in horizontal oil–water two-phase fl ow simulation wells were conducted using the Multiple Array Production Suite,which comprises a capacitance array tool(CAT)and a spinner array tool(SAT),and then the response characteristics of SAT and CAT in diff erent fl ow rates and water cut production conditions were studied.According to the response characteristics of CAT in diff erent water holdup ranges,interpolation imaging along the wellbore section determines the water holdup distribution,and then,the oil–water two-phase velocity fi eld in the fl ow section was reconstructed on the basis of the fl ow section water holdup distribution and the logging value of SAT and combined with the rheological equation of viscous fl uid,and the calculation method of the oil–water partial phase fl ow rate in the fl ow section was proposed.This new approach was applied in the experiment data calculations,and the results are basically consistent with the experimental data.The total fl ow rate and water holdup from the calculation are in agreement with the set values in the experiment,suggesting that the method has high accuracy.展开更多
Least-squares reverse time migration(LSRTM)can eliminate imaging artifacts in an iterative way based on the concept of inversion,and it can restore imaging amplitude step by step.LSRTM can provide a high-resolution mi...Least-squares reverse time migration(LSRTM)can eliminate imaging artifacts in an iterative way based on the concept of inversion,and it can restore imaging amplitude step by step.LSRTM can provide a high-resolution migration section and can be applied to irregular and poor-quality seismic data and achieve good results.Steeply dipping refl ectors and complex faults are imaged by using wavefi eld extrapolation based on a two-way wave equation.However,the high computational cost limits the method’s application in practice.A fast approach to realize LSRTM in the imaging domain is provided in this paper to reduce the computational cost signifi cantly and enhance its computational effi ciency.The method uses the Kronecker decomposition algorithm to estimate the Hessian matrix.A low-rank matrix can be used to calculate the Kronecker factor,which involves the calculation of Green’s function at the source and receiver point.The approach also avoids the direct construction of the whole Hessian matrix.Factorization-based LSRTM calculates the production of low-rank matrices instead of repeatedly calculating migration and demigration.Unlike traditional LSRTM,factorization-based LSRTM can reduce calculation costs considerably while maintaining comparable imaging quality.While having the same imaging eff ect,factorization-based LSRTM consumes half the running time of conventional LSRTM.In this regard,the application of factorization-based LSRTM has a promising advantage in reducing the computational cost.Ambient noise caused by this method can be removed by applying a commonly used fi ltering method without signifi cantly degrading the imaging quality.展开更多
In marine seismic exploration,especially in deep-water and hard ocean-bottom cases,free-surface multiples are strongly developed.Compared with primary waves,the wider illumination aperture of the multiples is benefici...In marine seismic exploration,especially in deep-water and hard ocean-bottom cases,free-surface multiples are strongly developed.Compared with primary waves,the wider illumination aperture of the multiples is beneficial for high-resolution seismic imaging.In this study,by introducing a new compound source composed of primaries and free-surface multiples and by ignoring internal multiples,we derive a new linearized forward problem(free-surface-multiple prediction model)under a weak-scattering assumption(i.e.,first-order Born approximation).On the basis of the new linearized problem,we propose a joint inversion-imaging method by simultaneously using the primaries and free-surface multiples under the general framework of least square inversion.To eliminate the crosstalk artifacts introduced by the cross-correlation of multiples with different orders,we prove that the crosstalk artifacts can be gradually eliminated during the inversion if a proper step length is selected.Synthetic-andfield-data tests demonstrate the effectiveness of the proposed method.展开更多
基金supported by the National Natural Science Foundation of China(Nos.41804124,41774138,41804121,41604109)China Academy of Sciences Strategic Leading Science and Technology Project(Grant Nos.XDA14020304,XDA14020302)+2 种基金Shandong Provincial Natural Science Foundation,China(No.ZR2019BD039)Shandong Province Postdoctoral Innovation Project(No.201901011)China Postdoctoral Science Foundation(Grant Nos.2019T120615,2018M632745)
文摘Acoustic reflection imaging in deep water wells is a new application scope for offshore hydrocarbon exploration.Two-dimensional(2 D)geological structure images can be obtained away from a one-dimensional(1 D)borehole using single-well acoustic reflection imaging.Based on the directivity of dipole source and four-component dipole data,one can achieve the azimuth detection and the three-dimensional(3 D)structural information around the wellbore can be obtained.We first perform matrix rotation on the field fourcomponent data.Then,a series of processing steps are applied to the rotated dipole data to obtain the reflector image.According to the above dipole shear-wave imaging principle,we used four-component cross-dipole logging data from a deviated well in the South China Sea to image geological structures within 50 m of a deviated well,which can delineate the structural configuration and determine its orientation.The configuration of near-borehole bedding boundaries and fault structures from shear-wave imaging results agrees with those from the Inline and Xline seismic profiles of the study area.In addition,the configuration and orientation of the fault structure images are consistent with regional stress maps and the results of the borehole stress anisotropy analysis.Furthermore,the dip azimuth of the bedding boundary images was determined using borehole wall resistivity data.Results of this study indicate that integrating borehole acoustic reflection with seismic imaging not only fills the gap between the two measurement scales but also accurately delineates geological structures in the borehole vicinity.
基金supported by the National Natural Science Foundation of China(Nos.41374122 and 41504100)
文摘The technology of simultaneous-source acquisition of seismic data excited by several sources can significantly improve the data collection efficiency. However, direct imaging of simultaneous-source data or blended data may introduce crosstalk noise and affect the imaging quality. To address this problem, we introduce a structure-oriented filtering operator as preconditioner into the multisource least-squares reverse-time migration (LSRTM). The structure-oriented filtering operator is a nonstationary filter along structural trends that suppresses crosstalk noise while maintaining structural information. The proposed method uses the conjugate-gradient method to minimize the mismatch between predicted and observed data, while effectively attenuating the interference noise caused by exciting several sources simultaneously. Numerical experiments using synthetic data suggest that the proposed method can suppress the crosstalk noise and produce highly accurate images.
基金supported by National Natural Science Foundation of China(41474115,42174155)Open Fund of Key Laboratory of Exploration Technologies for Oil and Gas Resources(Yangtze University)Ministry of Education of China(No K2018-02)。
文摘The distributions of local velocity and local phase holdup along the radial direction of pipes are complicated because of gravity differentiation,and the distribution of fluid velocity fi eld changes along the gravity direction in horizontal wells.Therefore,measuring the mixture flow and water holdup is difficult,resulting in poor interpretation accuracy of the production logging output profile.In this paper,oil–water two-phase flow dynamic simulation logging experiments in horizontal oil–water two-phase fl ow simulation wells were conducted using the Multiple Array Production Suite,which comprises a capacitance array tool(CAT)and a spinner array tool(SAT),and then the response characteristics of SAT and CAT in diff erent fl ow rates and water cut production conditions were studied.According to the response characteristics of CAT in diff erent water holdup ranges,interpolation imaging along the wellbore section determines the water holdup distribution,and then,the oil–water two-phase velocity fi eld in the fl ow section was reconstructed on the basis of the fl ow section water holdup distribution and the logging value of SAT and combined with the rheological equation of viscous fl uid,and the calculation method of the oil–water partial phase fl ow rate in the fl ow section was proposed.This new approach was applied in the experiment data calculations,and the results are basically consistent with the experimental data.The total fl ow rate and water holdup from the calculation are in agreement with the set values in the experiment,suggesting that the method has high accuracy.
基金funded by the National Natural Science Foundation of China (No.41574098&41630964)the Fundamental Research Funds for the Central Universities (No.18CX02059A)+3 种基金the Development Fund of Key Laboratory of Deep Oil&Gas (No. 20CX02111A)SINOPEC Key Laboratory of Geophysics open fund (No. wtyjy-wx2018-01-07)Shandong Natural Science Foundation of China(No. ZR2020MD048)the Major Scientific and Technological Projects of CNPC (No. ZD2019-183-003)
文摘Least-squares reverse time migration(LSRTM)can eliminate imaging artifacts in an iterative way based on the concept of inversion,and it can restore imaging amplitude step by step.LSRTM can provide a high-resolution migration section and can be applied to irregular and poor-quality seismic data and achieve good results.Steeply dipping refl ectors and complex faults are imaged by using wavefi eld extrapolation based on a two-way wave equation.However,the high computational cost limits the method’s application in practice.A fast approach to realize LSRTM in the imaging domain is provided in this paper to reduce the computational cost signifi cantly and enhance its computational effi ciency.The method uses the Kronecker decomposition algorithm to estimate the Hessian matrix.A low-rank matrix can be used to calculate the Kronecker factor,which involves the calculation of Green’s function at the source and receiver point.The approach also avoids the direct construction of the whole Hessian matrix.Factorization-based LSRTM calculates the production of low-rank matrices instead of repeatedly calculating migration and demigration.Unlike traditional LSRTM,factorization-based LSRTM can reduce calculation costs considerably while maintaining comparable imaging quality.While having the same imaging eff ect,factorization-based LSRTM consumes half the running time of conventional LSRTM.In this regard,the application of factorization-based LSRTM has a promising advantage in reducing the computational cost.Ambient noise caused by this method can be removed by applying a commonly used fi ltering method without signifi cantly degrading the imaging quality.
基金the sponsors of the WPI group for their financial supportfinancially supported by the National Key R&D Program of China (Grant Number: 2018YFA0702503, 2019YFC0312004)+2 种基金National Natural Science Foundation of China (Grant Number: 41774126)Southern Marine Science and Engineering Guangdong Laboratory (Zhanjiang) (ZJW-2019-04)National Science and Technology Major Project of China (Grant Number: 2016ZX05024-001, 2016ZX05006-002)。
文摘In marine seismic exploration,especially in deep-water and hard ocean-bottom cases,free-surface multiples are strongly developed.Compared with primary waves,the wider illumination aperture of the multiples is beneficial for high-resolution seismic imaging.In this study,by introducing a new compound source composed of primaries and free-surface multiples and by ignoring internal multiples,we derive a new linearized forward problem(free-surface-multiple prediction model)under a weak-scattering assumption(i.e.,first-order Born approximation).On the basis of the new linearized problem,we propose a joint inversion-imaging method by simultaneously using the primaries and free-surface multiples under the general framework of least square inversion.To eliminate the crosstalk artifacts introduced by the cross-correlation of multiples with different orders,we prove that the crosstalk artifacts can be gradually eliminated during the inversion if a proper step length is selected.Synthetic-andfield-data tests demonstrate the effectiveness of the proposed method.