Danping Cao

China University of Petroleum, East China

Author Statistics

Papers

Citation

H-Index

i-10 index

Research Trends

Author Order

Document Type

Co-Authors

Xingyao Yin

China University of Petroleum, East China

Bo Zhang

Aerospace Information Research Institute

Chao Jin

Khalifa University of Science and Technology

Zhiyu Hou

China University of Petroleum, East China

Yihuai Lou

Zhejiang University

Zhaolin Zhu

Sanya University

Kai Liang

China University of Petroleum, East China

Guochen Wu

China University of Petroleum, East China

Qiang Liu

China University of Petroleum, East China

Yitao Pu

Geological Survey of Alabama

Cooperative Institutions

China University of Petroleum, East China

China University of Petroleum, Beijing

Chinese Academy of Sciences

Qingdao National Laboratory for Marine Science and Technology

University of Chinese Academy of Sciences

Zhejiang University

China National Petroleum Corporation (China)

Institut de physique du globe de Paris

Sinopec (China)

Chengdu University of Technology

Author Statistics

Papers

Citation

H-Index

i-10 index

Research Field

Shear wave velocity estimation based on the rock physics model with the multi-variables information constraint

International Geophysical Conference, Qingdao, China, 17-20 April 2017 (2017)

Qijie Zhou Xingyao Yin Danping Cao Wenguo Sun Huixin Wang

Summary S-wave velocity plays important roles in pre-stack seismic inversion, AVO processing, also multicomponent seismic interpretation. To build a rock physics model is an alternative way to get the S-wave when lack of S-sonic logging data. Considering the Xu-White model which includes some parameters (such as mineral modulus, porosity, shale content, pore structure, pore fluid, etc.) that are inavailable. We proposed a method that retrieve the parameters from multi-variables logging data, like SP, GR, RT, etc.. Nevertheless, there is difference between interpretations from other single logging data. To deal with the contradiction, we give an inversion-based strategy. Firstly, specifying a range of the content of reservoir physical parameters. Then, constraining these physical property parameters through the P-wave velocity to improve the prediction accuracy of S-wave velocity and elastic parameters

Wave velocity

10.1190/igc2017-276

Cite

Citations (1)

Viscoelastic wave equations based on different rheological model in time domain modeling

Danping Cao

There are lots of viscoelastic wave equations expressed by different formalisms in time domain for the modeling of seismic wave with attenuation. In general, most viscoelastic wave equations are built based on two generalized rheological models, which are expressed by different relaxation parameters, attenuation variables and elastic modulus. With the combination of different attenuation variables (including the memory variables, the anelastic functions and the material-independent anelastic functions) and the elastic modulus (including the unrelaxed modulus and the relaxed modulus), 6 types of viscoelastic wave equations are built. Consequently, 24 kinds of formalisms can be expressed according to the displacement, displacement-stress, displacement-velocity-stress and velocity-stress wave equations. Several formalisms have been used in the literature by different authors. With the equivalent relations of the generalized rheological models, all 24 viscoelastic wave equations based on different rheological models are identical and can be translated into others easily. This work is beneficial to understand the different formulations of the viscoelastic wave equations and give authors more choice to use their favorite viscoelastic wave equations when modeling the seismic wave characteristics in viscoelastic medium.

Stress relaxation

10.1190/igcbeijing2014-172

Cite

Citations (1)

Joint viscoacoustic waveform inversion with the separated upgoing and downgoing wavefields in zero‐offset vertical seismic profile data

Geophysical Prospecting (2022)

Chao Jin Danping Cao Bing Zhou Yue Li Yuxi Wang

Abstract The attenuation of seismic waves propagating in reservoirs can be obtained accurately from the data analysis of vertical seismic profile in terms of the quality‐factor Q . The common methods usually use the downgoing wavefields in vertical seismic profile data. However, the downgoing wavefields consist of more than 90% energy of the spectrum of the vertical seismic profile data, making it difficult to estimate the viscoacoustic parameters accurately. Thus, a joint viscoacoustic waveform inversion of velocity and quality‐factor is proposed based on the multi‐objective functions and analysis of the difference between the results inverted from the separated upgoing and downgoing wavefields. A simple separating step is accomplished by the reflectivity method to obtain the individual wavefields in vertical seismic profile data, and then a joint inversion is carried out to make full use of the information of the individual wavefields and improve the convergence of viscoacoustic full‐waveform inversion. The sensitivity analysis of the different wavefields to the velocity and quality‐factor shows that the upgoing and downgoing wavefields contribute differently to the viscoacoustic parameters. A numerical example validates our method can improve the accuracy of viscoacoustic parameters compared with the direct inversion using full wavefield and the separate inversion using upgoing or downgoing wavefield. The application on real field data indicates our method can recover a reliable viscoacoustic model, which helps reservoir appraisal.

Vertical seismic profile

Anelastic attenuation factor

Seismic migration

10.1111/1365-2478.13281

Cite

Citations (2)

Pure qP-wave equation in VTI media based on approximate method of completing the square

Kai Liang Shangrao Sun Danping Cao Xingyao Yin

The elastic wave propagation is coupled in VTI media. Different approximations are used to obtain single wave equation, such as acoustic approximation and weak anisotropic approximation. In this paper, we start from the exact elastic wave phase velocities in VTI media, and propose an approximate qP-wave phase velocity based on the approximate method of completing the square. Then we derive the corresponding pure qP-wave equation in VTI media, and choosing hybrid method as the solver. The accuracy of corresponding approximate phase velocity is tested. The results show a good accuracy in weak and strong anisotropy. Numerical examples show that our proposed equation is stable and free of S wave. The wavefront of pure qP-wave has good similarity with that of conventional elastic wave.

Solver

10.1190/segam2021-3579895.1

Cite

Citations (0)

Impedance inversion method constrained with crosswell seismic data

Danping Cao Xingyao Yin Fanchang Zhang

Impedance inversion of surface seismic data is an inherently non‐unique problem. The conventional seismic inversion resolution is limited to characterize the subsurface reservoirs because of the narrow frequency bandwidth. The high frequency components of crosswell seismic data are rich enough to connect the surface seismic data and the logging curve. Based on the sparse spike inversion assumption, the crosswell seismic data is introduced into the surface seismic impedance inversion flow. The objective function of impedance inversion method is built constrained with the crosswell seismic data and well logging data. The field data inversion result shows that the crosswell seismic data constrained inversion method can improve the resolution of the conventional well‐logging constrained inversion method. It is a new way to solve the geology problems with kinds of geophysical data.

10.1190/1.3603786

Cite

Citations (2)

Non-orthogonal beam coordinate system wave propagation and reverse time migration

Journal of Geophysics and Engineering (2019)

Zhaolin Zhu Danping Cao Bangyu Wu Xingyao Yin Ying Wang

Abstract Grid size has a significant influence on the computation efficiency and accuracy of finite-difference seismic modeling and can change the workload of reverse time migration (RTM) remarkably. This paper proposes a non-orthogonal analytical coordinate system, beam coordinate system (BCS), for the solution of seismic wave propagation and RTM. Starting with an optical Gaussian beam width equation, we expand the representation on vertically variable velocity media, which is the most common scenario in seismic exploration. The BCS based on this representation can be used to implement an irregular-grid increment finite-difference that improves the efficiency of RTM. Based on the Laplacian expression in Riemannian space, we derive the wave equation in the BCS. The new coordinate system can generate an irregular grid with increment increasing vertically along depth. Through paraxial ray tracing, it can be extended to non-analytical beam coordinate system (NBCS). Experiments for the RTM on the Marmousi model with the BCS demonstrate that the proposed method improves the efficiency about 52% while maintaining good image quality.

Seismic migration

Gaussian beam

Finite difference

Representation

Coordinate Descent

10.1093/jge/gxz067

Cite

Citations (4)

Q estimation with the upgoing and downgoing waves for VSP

Danping Cao Yue Li Chao Jin

Quality factor Q is a key parameter describing seismic wave attenuation. Estimating Q accurately is important to seismic data processing and reservoir characterization. VSP records have the information of upgoing and downgoing waves, and can show the changes of seismic wavelet at different depths, which has certain advantages in Q value estimation. Downgoing and upgoing waves are coming from above and below the receiver containing the strata information, respectively. Downgoing wave has always been the key in Q estimation, while the upgoing wave is seldom used. Thus, the information contained in the upgoing wave is not fully utilized, leading to problems of multi-solution and low accuracy. To fix these problems, this study proposes a joint estimation method of Q, which uses both upgoing and downgoing waves. The nonlinear relation among the upgoing wave, downgoing wave, and Q value is jointly used to construct an objective function, which measures the energy difference between the actual and calculated amplitude spectra. As the energy between upgoing and downgoing waves is not balanced, a weight coefficient is used to balance the weight of the two kinds of waves in Q estimation. The optimal Q can be estimated by minimizing the objective function. Numerical simulation results indicate that the proposed method can provide more accurate results. The application on field data demonstrates that good results can be achieved, and the estimated Q has great potential in reservoir characterization. Presentation Date: Tuesday, October 13, 2020 Session Start Time: 8:30 AM Presentation Time: 10:35 AM Location: 360D Presentation Type: Oral

10.1190/segam2020-3427591.1

Cite

Citations (0)

The small-scale forward modeling method for large models based on the wave equation

Journal of Geophysics and Engineering (2015)

Yitao Pu Xingyao Yin Danping Cao Chong Zhang

The small-scale modeling of an underground geological model can avoid the equivalent large-scale assumptions and show the seismic reflection characteristics more precisely. Considering the large number of perfectly matched layer (PML) boundaries in small-scale modeling of large models, we analyze the absorption effect of the PML boundary in different spatial spacing, and apply variable-grid technology to small-scale seismic modeling. Using small grids to compute the model region and large grids to compute the PML boundary, we successfully reduce the amount of PML boundary layers in small-scale seismic modeling. This method can reduce data storage and improve computational efficiency dramatically.

Perfectly Matched Layer

Reflection

Scale model

10.1088/1742-2132/12/6/988

Cite

Citations (1)

An Adjoint-Based Hybrid Computational Method for Crosswell Seismic Inversion

Computing in Science & Engineering (2014)

Danping Cao Wenyuan Liao

Waveform tomography, or full waveform inversion (FWI), is in general better than ray tomography; however, a reliable initial model is usually required to ensure success. Here, the authors designed a cascade-like hybrid tomography technology to solve the crosswell seismic inversion problem. The new method is a combination of several widely used tomography technologies, implemented in a sequence, where the cheaper method is used as the initial model feeder for the more expensive method. They start with the Radon-transform-based back-projection method; then the Linear Iterative Reconstruction method, such as selective internal radiation therapy, is adopted to update this initial model; which is further improved by the nonlinear, Gradient-based Eikonal Equation Tomography method. Finally, the velocity model reconstructed from the previous multiple ray tomography methods is used as the initial model for the FWI.

Radon transform

Seismic Tomography

10.1109/mcse.2014.1

Cite

Citations (2)

Velocity and Q estimation from the separated upgoing and downgoing wavefields in VSP data

Chao Jin Danping Cao Bing Zhou

Summary Quality factor (Q-factor) evaluates the attenuation of seismic wave propagation, playing a fundamental role of reservoir characterization, which can be obtained accurately from Vertical Seismic Profile (VSP). The common methods usually use the downgoing wavefields in VSP data. However, the downgoing wavefields consist of more than 90% energy of the spectrum of the VSP data due to the energy fraction of the upgoing and downgoing wavefields, which makes difficult to estimate the viscoacoustic parameters accurately. Thus, a joint viscoacoustic waveform inversion of velocity and Q-factor is proposed to measure the difference between the separated upgoing and downgoing wavefields in VSP data based on the multi-objective functions. A simple separating step is accomplished by the reflectivity method to obtain the pure individual wavefields in VSP data, and then a joint inversion step is carried out to make full use of the information of the individual wavefields and improve the convergence of viscoacoustic waveform inversion. The sensitivity analysis about the velocity and Q-factor shows that the upgoing and downgoing wavefields contribute differently to the viscoacoustic parameters. Numerical examples and a field test indicate the accuracy and efficiency of the proposed method.

10.3997/2214-4609.202113294

Cite

Citations (0)