Angle-domain common-image gathers （ADCIGs） transformed from the shot- domain common-offset gathers are input to migration velocity analysis （MVA） and prestack inversion. ADCIGs are non-illusion prestack inversion ...Angle-domain common-image gathers （ADCIGs） transformed from the shot- domain common-offset gathers are input to migration velocity analysis （MVA） and prestack inversion. ADCIGs are non-illusion prestack inversion gathers, and thus, accurate. We studied the extraction of elastic-wave ADCIGs based on amplitude-preserving elastic-wave reverse- time migration for calculating the incidence angle of P- and S-waves at each image point and for different source locations. The P- and S-waves share the same incident angle, namely the incident angle of the source P-waves. The angle of incidence of the source P-wavefield was the difference between the source P-wave propagation angle and the reflector dips. The propagation angle of the source P-waves was obtained from the polarization vector of the decomposed P-waves. The reflectors＇ normal direction angle was obtained using the complex wavenumber of the stacked reverse-time migration （RTM） images. The ADCIGs of P- and S-waves were obtained by rearranging the common-shot migration gathers based on the incident angle. We used a horizontally layered model, the graben medium model, and part of the Marmousi-II elastic model and field data to test the proposed algorithm. The results suggested that the proposed method can efficiently extract the P- and S-wave ADCIGs of the elastic-wave reverse-time migration, the P- and S-wave incident angle, and the angle-gather amplitude fidelity, and improve the MVA and prestack inversion.展开更多
Conventional full-waveform inversion is computationally intensive because it considers all shots in each iteration. To tackle this, we establish the number of shots needed and propose multiscale inversion in the frequ...Conventional full-waveform inversion is computationally intensive because it considers all shots in each iteration. To tackle this, we establish the number of shots needed and propose multiscale inversion in the frequency domain while using only the shots that are positively correlated with frequency. When using low-frequency data, the method considers only a small number of shots and raw data. More shots are used with increasing frequency. The random-in-group subsampling method is used to rotate the shots between iterations and avoid the loss of shot information. By reducing the number of shots in the inversion, we decrease the computational cost. There is no crosstalk between shots, no noise addition, and no observational limits. Numerical modeling suggests that the proposed method reduces the computing time, is more robust to noise, and produces better velocity models when using data with noise.展开更多
Accuracy of angle-domain common-image gathers(ADCIGs)is the key to multiwave AVA inversion and migration velocity analysis,and of which Poynting vectors of pure P-and S-wave are the decisive factors in obtaining multi...Accuracy of angle-domain common-image gathers(ADCIGs)is the key to multiwave AVA inversion and migration velocity analysis,and of which Poynting vectors of pure P-and S-wave are the decisive factors in obtaining multi-component seismic data ADCIGs.A Poynting vector can be obtained from conventional velocity-stress elastic wave equations,but it focused on the propagation direction of mixed P-and S-wave fields,and neither on the propagation direction of the P-wave nor the direction of the S-wave.The Poynting vectors of pure P-or pure S-wave can be calculated from first-order velocity-dilatation-rotation equations.This study presents a method of extracting ADCIGs based on first order velocitydilatation-rotation elastic wave equations reverse-time migration algorithm.The method is as follows:calculating the pure P-wave Poynting vector of source and receiver wavefields by multiplication of P-wave particle-velocity vector and dilatation scalar,calculating the pure S-wave Poynting vector by vector multiplying S-wave particle-velocity vector and rotation vector,selecting the Poynting vector at the time of maximum P-wave energy of source wavefield as the propagation direction of incident P-wave,and obtaining the reflected P-wave(or converted S-wave)propagation direction of the receiver wavefield by the Poynting vector at the time of maximum P-(S-)wave energy in each grid point.Then,the P-wave incident angle is computed by the two propagation directions.Thus,the P-and S-wave ADGICs can obtained Numerical tests show that the proposed method can accurately compute the propagation direction and incident angle of the source and receiver wavefields,thereby achieving high-precision extraction of P-and S-wave ADGICs.展开更多
An air-gun source is the most commonly used excitation method in off shore seismic exploration.The excitation characteristics of an air-gun source aff ect seismic data quality.Far-field wavelet simulation is an import...An air-gun source is the most commonly used excitation method in off shore seismic exploration.The excitation characteristics of an air-gun source aff ect seismic data quality.Far-field wavelet simulation is an important approach to study these characteristics.Compared to the measured wavelet,far-field wavelet simulation based on a traditional bubble-motion equation and ideal gas wavelet model has some disadvantages,such as a greater amplitude and smaller pulse attenuation velocity.Here,we start from the linear acoustic wave equation in the spherical coordinate system to deduce an improved,simpler bubble-motion equation and develop a Van der Waals gas wavelet model based on this equation.Unlike the existing methods,our method considers the high-pressure environment during actual excitation,heat exchange between the bubble and outside water,and change in the air fl ow at the muzzle.The results show that the far-fi eld wavelet simulated using this model is closer to the measured wavelet than that of the ideal gas wavelet model.At the same time,our method has a more succinct equation and a higher calculation effi ciency.展开更多
基金supported by Financially Supported by Qingdao National Laboratory for Marine Science and Technology(No.QNLM2016ORP0206)National Science and Technology Major Project(No.2016ZX05027-002)+6 种基金China Postdoctoral Science Foundation(No.2017M612219)National Key R&D Plan(Nos.2017YFC0306706 and 2017YFC0307400)Financially Supported by Qingdao National Laboratory for Marine Science and Technology(No.QNLM201708)Natural Science Foundation of Shandong Province(No.ZR2016DB10)National Natural Science Foundation of China(Nos.41674118,41504109,and 41506084)Key Laboratory of Submarine Geosciences Foundation of SOA(No.KLSG1603)Qingdao Municipal Applied Research Projects(No.2016238)
文摘Angle-domain common-image gathers （ADCIGs） transformed from the shot- domain common-offset gathers are input to migration velocity analysis （MVA） and prestack inversion. ADCIGs are non-illusion prestack inversion gathers, and thus, accurate. We studied the extraction of elastic-wave ADCIGs based on amplitude-preserving elastic-wave reverse- time migration for calculating the incidence angle of P- and S-waves at each image point and for different source locations. The P- and S-waves share the same incident angle, namely the incident angle of the source P-waves. The angle of incidence of the source P-wavefield was the difference between the source P-wave propagation angle and the reflector dips. The propagation angle of the source P-waves was obtained from the polarization vector of the decomposed P-waves. The reflectors＇ normal direction angle was obtained using the complex wavenumber of the stacked reverse-time migration （RTM） images. The ADCIGs of P- and S-waves were obtained by rearranging the common-shot migration gathers based on the incident angle. We used a horizontally layered model, the graben medium model, and part of the Marmousi-II elastic model and field data to test the proposed algorithm. The results suggested that the proposed method can efficiently extract the P- and S-wave ADCIGs of the elastic-wave reverse-time migration, the P- and S-wave incident angle, and the angle-gather amplitude fidelity, and improve the MVA and prestack inversion.
基金financially supported by the Fundamental Research Funds for the Central Universities(No.201822011)the National Natural Science Foundation of China(No.41674118)the National Science and Technology Major Project(No.2016ZX05027002)
文摘Conventional full-waveform inversion is computationally intensive because it considers all shots in each iteration. To tackle this, we establish the number of shots needed and propose multiscale inversion in the frequency domain while using only the shots that are positively correlated with frequency. When using low-frequency data, the method considers only a small number of shots and raw data. More shots are used with increasing frequency. The random-in-group subsampling method is used to rotate the shots between iterations and avoid the loss of shot information. By reducing the number of shots in the inversion, we decrease the computational cost. There is no crosstalk between shots, no noise addition, and no observational limits. Numerical modeling suggests that the proposed method reduces the computing time, is more robust to noise, and produces better velocity models when using data with noise.
基金financially supported by the Fundamental Research Funds for the Central Universities(No.201822011)the National Key R&D Program of China(No.2018YFC1405900)+1 种基金the National Natural Science Foundation of China(Nos.41674118 and 41574105)the National Science and Technology Major Project(No.2016ZX05027002)。
文摘Accuracy of angle-domain common-image gathers(ADCIGs)is the key to multiwave AVA inversion and migration velocity analysis,and of which Poynting vectors of pure P-and S-wave are the decisive factors in obtaining multi-component seismic data ADCIGs.A Poynting vector can be obtained from conventional velocity-stress elastic wave equations,but it focused on the propagation direction of mixed P-and S-wave fields,and neither on the propagation direction of the P-wave nor the direction of the S-wave.The Poynting vectors of pure P-or pure S-wave can be calculated from first-order velocity-dilatation-rotation equations.This study presents a method of extracting ADCIGs based on first order velocitydilatation-rotation elastic wave equations reverse-time migration algorithm.The method is as follows:calculating the pure P-wave Poynting vector of source and receiver wavefields by multiplication of P-wave particle-velocity vector and dilatation scalar,calculating the pure S-wave Poynting vector by vector multiplying S-wave particle-velocity vector and rotation vector,selecting the Poynting vector at the time of maximum P-wave energy of source wavefield as the propagation direction of incident P-wave,and obtaining the reflected P-wave(or converted S-wave)propagation direction of the receiver wavefield by the Poynting vector at the time of maximum P-(S-)wave energy in each grid point.Then,the P-wave incident angle is computed by the two propagation directions.Thus,the P-and S-wave ADGICs can obtained Numerical tests show that the proposed method can accurately compute the propagation direction and incident angle of the source and receiver wavefields,thereby achieving high-precision extraction of P-and S-wave ADGICs.
基金supported by National Natural Science Foundation of China (No. 41674118)
文摘An air-gun source is the most commonly used excitation method in off shore seismic exploration.The excitation characteristics of an air-gun source aff ect seismic data quality.Far-field wavelet simulation is an important approach to study these characteristics.Compared to the measured wavelet,far-field wavelet simulation based on a traditional bubble-motion equation and ideal gas wavelet model has some disadvantages,such as a greater amplitude and smaller pulse attenuation velocity.Here,we start from the linear acoustic wave equation in the spherical coordinate system to deduce an improved,simpler bubble-motion equation and develop a Van der Waals gas wavelet model based on this equation.Unlike the existing methods,our method considers the high-pressure environment during actual excitation,heat exchange between the bubble and outside water,and change in the air fl ow at the muzzle.The results show that the far-fi eld wavelet simulated using this model is closer to the measured wavelet than that of the ideal gas wavelet model.At the same time,our method has a more succinct equation and a higher calculation effi ciency.
