When seismic exploration is conducted in a special geological environment such as a tunnel space,the traditional imaging method in the Cartesian coordinate system cannot accurately discretize the air column in that en...When seismic exploration is conducted in a special geological environment such as a tunnel space,the traditional imaging method in the Cartesian coordinate system cannot accurately discretize the air column in that environment.Thus,obtaining Thus,obtaining highquality imaging results is diffi cult.Therefore,an elastic-wave reverse-time migration method based on the polar coordinate system is proposed.In this method,three boundary conditions exist:outer,inner,and corner boundaries.In the outer boundary,the polar-coordinated absorbing boundary in the radial direction is used to suppress the artifi cial-boundary refl ection.The free-surface boundary condition is adopted in the tunnel space at the inner boundary.In the angular boundaries,we use two diff erent boundary conditions for two cases.The air column in the tunnel space is usually not an irregular circle.Therefore,the irregular tunnelspace geological body in the polar coordinate system is meshed into curvilinear grids and transformed into a regular one in an auxiliary polar coordinate system using the mapping method.Finally,elastic reverse-time migration technology is applied into the auxiliary polar coordinate system.In the numerical examples,two typical models are used to test the proposed method,which verify that the proposed method can obtain accurate images from the datasets in the tunnel space.展开更多
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.展开更多
基金financially supported by the National Natural Science Foundation of China (grant Nos. 41904101 and 41774133)Natural Science Foundation of Shandong Province (grant No. ZR2019QD004)+1 种基金Fundamental Research Funds for the Central Universities (grant No. 19CX02010A)the Open Funds of SINOPEC Key Laboratory of Geophysics (grant No. wtyjy-wx2019-01-03)。
文摘When seismic exploration is conducted in a special geological environment such as a tunnel space,the traditional imaging method in the Cartesian coordinate system cannot accurately discretize the air column in that environment.Thus,obtaining Thus,obtaining highquality imaging results is diffi cult.Therefore,an elastic-wave reverse-time migration method based on the polar coordinate system is proposed.In this method,three boundary conditions exist:outer,inner,and corner boundaries.In the outer boundary,the polar-coordinated absorbing boundary in the radial direction is used to suppress the artifi cial-boundary refl ection.The free-surface boundary condition is adopted in the tunnel space at the inner boundary.In the angular boundaries,we use two diff erent boundary conditions for two cases.The air column in the tunnel space is usually not an irregular circle.Therefore,the irregular tunnelspace geological body in the polar coordinate system is meshed into curvilinear grids and transformed into a regular one in an auxiliary polar coordinate system using the mapping method.Finally,elastic reverse-time migration technology is applied into the auxiliary polar coordinate system.In the numerical examples,two typical models are used to test the proposed method,which verify that the proposed method can obtain accurate images from the datasets in the tunnel space.
基金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.
