P-and SV-wave dispersion and attenuation have been extensively investigated in saturated poroelastic media with aligned fractures.However,there are few existing models that incorporate the multiple wave attenuation me...P-and SV-wave dispersion and attenuation have been extensively investigated in saturated poroelastic media with aligned fractures.However,there are few existing models that incorporate the multiple wave attenuation mechanisms from the microscopic scale to the macroscopic scale.Hence,in this work,we developed a unified model to incorporate the wave attenuation mechanisms at different scales,which includes the microscopic squirt flow between the microcracks and pores,the mesoscopic wave-induced fluid flow between fractures and background(FB-WIFF),and the macroscopic Biot's global flow and elastic scattering(ES)from the fractures.Using Tang's modified Biot's theory and the mixed-boundary conditions,we derived the exact frequency-dependent solutions of the scattering problem for a single penny-shaped fracture with oblique incident P-and SV-waves.We then developed theoretical models for a set of aligned fractures and randomly oriented fractures using the Foldy approximation.The results indicated that microcrack squirt flow considerably influences the dispersion and attenuation of P-and SV-wave velocities.The coupling effects of microcrack squirt flow with the FB-WIFF and ES of fractures cause much higher velocity dispersion and attenuation for P waves than for SV waves.Randomly oriented fractures substantially reduce the attenuation caused by the FB-WIFF and ES,particularly for the ES attenuation of SV waves.Through a comparison with existing models in the limiting cases and previous experimental measurements,we validated our model.展开更多
In marine seismic exploration, ocean-bottom cable techniques accurately record the multicomponent seismic wavefield; however, the seismic wave propagation in fluid–solid media cannot be simulated by a single wave equ...In marine seismic exploration, ocean-bottom cable techniques accurately record the multicomponent seismic wavefield; however, the seismic wave propagation in fluid–solid media cannot be simulated by a single wave equation. In addition, when the seabed interface is irregular, traditional finite-difference schemes cannot simulate the seismic wave propagation across the irregular seabed interface. Therefore, an acoustic–elastic forward modeling and vector-based P-and S-wave separation method is proposed. In this method, we divide the fluid–solid elastic media with irregular interface into orthogonal grids and map the irregular interface in the Cartesian coordinates system into a horizontal interface in the curvilinear coordinates system of the computational domain using coordinates transformation. The acoustic and elastic wave equations in the curvilinear coordinates system are applied to the fluid and solid medium, respectively. At the irregular interface, the two equations are combined into an acoustic–elastic equation in the curvilinear coordinates system. We next introduce a full staggered-grid scheme to improve the stability of the numerical simulation. Thus, separate P-and S-wave equations in the curvilinear coordinates system are derived to realize the P-and S-wave separation method.展开更多
Joint PP–PS inversion offers better accuracy and resolution than conventional P-wave inversion. P-and S-wave elastic moduli determined through data inversions are key parameters for reservoir evaluation and fluid cha...Joint PP–PS inversion offers better accuracy and resolution than conventional P-wave inversion. P-and S-wave elastic moduli determined through data inversions are key parameters for reservoir evaluation and fluid characterization. In this paper, starting with the exact Zoeppritz equation that relates P-and S-wave moduli, a coefficient that describes the reflections of P-and converted waves is established. This method effectively avoids error introduced by approximations or indirect calculations, thus improving the accuracy of the inversion results. Considering that the inversion problem is ill-posed and that the forward operator is nonlinear, prior constraints on the model parameters and modified low-frequency constraints are also introduced to the objective function to make the problem more tractable. This modified objective function is solved over many iterations to continuously optimize the background values of the velocity ratio, which increases the stability of the inversion process. Tests of various models show that the method effectively improves the accuracy and stability of extracting P and S-wave moduli from underdetermined data. This method can be applied to provide inferences for reservoir exploration and fluid extraction.展开更多
Elastic reverse time migration(RTM)uses the elastic wave equation to extrapolate multicomponent seismic data to the subsurface and separate the elastic wavefield into P-and S-waves.P-and S-wave separation is a necessa...Elastic reverse time migration(RTM)uses the elastic wave equation to extrapolate multicomponent seismic data to the subsurface and separate the elastic wavefield into P-and S-waves.P-and S-wave separation is a necessary step in elastic RTM to avoid crosstalk between coupled wavefields.However,the current curl-divergence operator-based separation method has a polarity reversal problem in PS imaging,and vector separation methods often have separation artifacts at the interface,which affects the quality of the imaging stack.We propose a non-artifact P-and S-wave separation method based on the first-order velocity-strain equation.This equation is used for wavefield extrapolation and separation in the first-order staggered-grid finite-difference scheme,and the storage and calculation amounts are consistent with the classical first-order velocity-stress equation.The separation equation does not calculate the partial derivatives of the elastic parameters,and thus,there is no artifact in the separated Pand S-waves.During wavefield extrapolation,the dynamic characteristics of the reflected wave undergo some changes,but the transmitted wavefield is accurate;therefore,it does not affect the dynamic characteristics of the final migration imaging.Through numerical examples of 2 D simple models,part SEAM model,BP model,and 3 D 4-layer model,different wavefield separation methods and corresponding elastic RTM imaging results are analyzed.We found that the velocity-strain based elastic RTM can image subsurface structures well,without spike artifacts caused by separation artifacts,and without polarity reversal phenomenon of the PS imaging.展开更多
West Java in the western part of the Sunda Arc has a relatively high seismicity due to subduction activity and faults.In this study,double-difference tomography was used to obtain the 3D velocity tomograms of P and S ...West Java in the western part of the Sunda Arc has a relatively high seismicity due to subduction activity and faults.In this study,double-difference tomography was used to obtain the 3D velocity tomograms of P and S waves beneath the western part of Java.To infer the geometry of the structure beneath the study area,precise earthquake hypo・center determination was first performed before tomographic imaging.For this,earthquake waveform data were extracted from the regional Meteorological,Climatological,Geophysical Agency(BMKG)network of Indonesia from South Sumatra to Central Java.The P and S arrival times for about 1,000 events in the period April 2009 to July 2016 were selected,the key features being events of magnitude>3,azimuthal gap<210°and number of phases>8.A nonlinear method using the oct-tree sampling algorithm from the NonLinLoc program was employed to determine the earthquake hypocenters.The hypocenter locations were then relocated using double-difference tomography(tomoDD).A significant reduction of travel-time(root mean square basis)and a better clustering of earthquakes were achieved which correlated well with the geological structure in West Java.Double-difference tomography was found to give a clear velocity structure,especially beneath the volcanic arc area,i.e.,under Mt Anak Krakatau,Mt Salak and the mountains complex in the southern part of West Java.Low velocity anomalies for the P and S waves as well as the vp/vs ratio below the volcanoes indicated possible partial melting of the upper mantle which ascended from the subducted slab beneath the volcanic arc.展开更多
Based on Hudson’s theoretical hypothesis of equivalent fracture model,inserting aligned round chips in solid model can simulate fractured media. The effect of fractures on the propagation of P and S waves can be obse...Based on Hudson’s theoretical hypothesis of equivalent fracture model,inserting aligned round chips in solid model can simulate fractured media. The effect of fractures on the propagation of P and S waves can be observed by changing the fracture thickness. The base model is made of epoxy resin,and the material of fractures is a kind of low-velocity mixture containing silicon rubber. With constant diameter and number of fractures in each model,one group of models can be formed through changing the thickness of fracture. These models have the same fracture density. By using the ultrasonic pulse transmission method,the experiment records time and waveform of P and S waves in the direction parallel and perpendicular to the fracture orientation. The result shows that,with the same fracture density,changing fracture aperture will affect both velocity and amplitude of P and S waves,and the effect on P-wave amplitude is much greater than that on the velocity. Moreover,the variation in velocity of S wave is more obvious in the slow shear wave (S2),while the variation in amplitude is more obvious in the fast shear wave (S1). These properties of wave propagation are useful for seismic data processing and interpretation.展开更多
The partial derivative equations of Zoeppritz equations are established and the derivatives of each matrix entry with respect to wave vectors are derived in this paper.By solving the partial derivative equations we ob...The partial derivative equations of Zoeppritz equations are established and the derivatives of each matrix entry with respect to wave vectors are derived in this paper.By solving the partial derivative equations we obtained the partial derivatives of seismic wave reflection coefficients with respect to wave vectors,and computed the Goos-Hnchen shift for reflected P-and VS-waves.By plotting the curves of Goos-Hnchen shift,we gained some new insight into the lateral shift of seismic reflection wave.The lateral shifts are very large for glancing wave or the wave of the incidence angle near the critical angle,meaning that the seismic wave propagates a long distance along the reflection interface before returning to the first medium.For the reflection waves of incidence angles away from the critical angle,the lateral shift is in the same order of magnitude as the wavelength.The lateral shift varies significantly with different reflection interfaces.For example,the reflected P-wave has a negative shift at the reflection interface between mudstone and sandstone.The reflected VS-wave has a large lateral shift at or near the critical angle.The lateral shift of the reflected VS-wave tends to be zero when the incidence angle approaches 90°.These observations suggest that Goos-Hnchen effect has a great influence on the reflection wave of wide-angles.The correction for the error caused by Goos-Hnchen effect,therefore,should be made before seismic data processing,such as the depth migration and the normal-moveout correction.With the theoretical foundation established in this paper,we can further study the correction of Goos-Hnchen effect for the reflection wave of large incidence angle.展开更多
基金This work was supported by the Laoshan National Laboratory Science and Technology Innovation Project(No.LSKJ202203407)the National Natural Science Foundation of China(Grant Nos.42174145,41821002,42274146)+1 种基金Guangdong Provincial Key Laboratory of Geophysical High-resolution Imaging Technology(2022B1212010002)Shenzhen Stable Support Plan Program for Higher Education Institutions(20220815110144003).
文摘P-and SV-wave dispersion and attenuation have been extensively investigated in saturated poroelastic media with aligned fractures.However,there are few existing models that incorporate the multiple wave attenuation mechanisms from the microscopic scale to the macroscopic scale.Hence,in this work,we developed a unified model to incorporate the wave attenuation mechanisms at different scales,which includes the microscopic squirt flow between the microcracks and pores,the mesoscopic wave-induced fluid flow between fractures and background(FB-WIFF),and the macroscopic Biot's global flow and elastic scattering(ES)from the fractures.Using Tang's modified Biot's theory and the mixed-boundary conditions,we derived the exact frequency-dependent solutions of the scattering problem for a single penny-shaped fracture with oblique incident P-and SV-waves.We then developed theoretical models for a set of aligned fractures and randomly oriented fractures using the Foldy approximation.The results indicated that microcrack squirt flow considerably influences the dispersion and attenuation of P-and SV-wave velocities.The coupling effects of microcrack squirt flow with the FB-WIFF and ES of fractures cause much higher velocity dispersion and attenuation for P waves than for SV waves.Randomly oriented fractures substantially reduce the attenuation caused by the FB-WIFF and ES,particularly for the ES attenuation of SV waves.Through a comparison with existing models in the limiting cases and previous experimental measurements,we validated our model.
基金financially supported by the Natural Science Foundation of China(No.41774133)the Open Funds of SINOPEC Key Laboratory of Geophysics(No.wtyjy-wx2017-01-04)National Science and Technology Major Project of the Ministry of Science and Technology of China(No.2016ZX05024-003-011)
文摘In marine seismic exploration, ocean-bottom cable techniques accurately record the multicomponent seismic wavefield; however, the seismic wave propagation in fluid–solid media cannot be simulated by a single wave equation. In addition, when the seabed interface is irregular, traditional finite-difference schemes cannot simulate the seismic wave propagation across the irregular seabed interface. Therefore, an acoustic–elastic forward modeling and vector-based P-and S-wave separation method is proposed. In this method, we divide the fluid–solid elastic media with irregular interface into orthogonal grids and map the irregular interface in the Cartesian coordinates system into a horizontal interface in the curvilinear coordinates system of the computational domain using coordinates transformation. The acoustic and elastic wave equations in the curvilinear coordinates system are applied to the fluid and solid medium, respectively. At the irregular interface, the two equations are combined into an acoustic–elastic equation in the curvilinear coordinates system. We next introduce a full staggered-grid scheme to improve the stability of the numerical simulation. Thus, separate P-and S-wave equations in the curvilinear coordinates system are derived to realize the P-and S-wave separation method.
基金supported by the National Science and Technology Major Project(No.2016ZX05047-002-001)
文摘Joint PP–PS inversion offers better accuracy and resolution than conventional P-wave inversion. P-and S-wave elastic moduli determined through data inversions are key parameters for reservoir evaluation and fluid characterization. In this paper, starting with the exact Zoeppritz equation that relates P-and S-wave moduli, a coefficient that describes the reflections of P-and converted waves is established. This method effectively avoids error introduced by approximations or indirect calculations, thus improving the accuracy of the inversion results. Considering that the inversion problem is ill-posed and that the forward operator is nonlinear, prior constraints on the model parameters and modified low-frequency constraints are also introduced to the objective function to make the problem more tractable. This modified objective function is solved over many iterations to continuously optimize the background values of the velocity ratio, which increases the stability of the inversion process. Tests of various models show that the method effectively improves the accuracy and stability of extracting P and S-wave moduli from underdetermined data. This method can be applied to provide inferences for reservoir exploration and fluid extraction.
基金supported by the National Natural Science Foundation of China,Grant No.41774142
文摘Elastic reverse time migration(RTM)uses the elastic wave equation to extrapolate multicomponent seismic data to the subsurface and separate the elastic wavefield into P-and S-waves.P-and S-wave separation is a necessary step in elastic RTM to avoid crosstalk between coupled wavefields.However,the current curl-divergence operator-based separation method has a polarity reversal problem in PS imaging,and vector separation methods often have separation artifacts at the interface,which affects the quality of the imaging stack.We propose a non-artifact P-and S-wave separation method based on the first-order velocity-strain equation.This equation is used for wavefield extrapolation and separation in the first-order staggered-grid finite-difference scheme,and the storage and calculation amounts are consistent with the classical first-order velocity-stress equation.The separation equation does not calculate the partial derivatives of the elastic parameters,and thus,there is no artifact in the separated Pand S-waves.During wavefield extrapolation,the dynamic characteristics of the reflected wave undergo some changes,but the transmitted wavefield is accurate;therefore,it does not affect the dynamic characteristics of the final migration imaging.Through numerical examples of 2 D simple models,part SEAM model,BP model,and 3 D 4-layer model,different wavefield separation methods and corresponding elastic RTM imaging results are analyzed.We found that the velocity-strain based elastic RTM can image subsurface structures well,without spike artifacts caused by separation artifacts,and without polarity reversal phenomenon of the PS imaging.
基金the Directorate General of Resources for Science Technologythe Higher Education of the Republic of Indonesia for granting a PMDSU scholarship to SR
文摘West Java in the western part of the Sunda Arc has a relatively high seismicity due to subduction activity and faults.In this study,double-difference tomography was used to obtain the 3D velocity tomograms of P and S waves beneath the western part of Java.To infer the geometry of the structure beneath the study area,precise earthquake hypo・center determination was first performed before tomographic imaging.For this,earthquake waveform data were extracted from the regional Meteorological,Climatological,Geophysical Agency(BMKG)network of Indonesia from South Sumatra to Central Java.The P and S arrival times for about 1,000 events in the period April 2009 to July 2016 were selected,the key features being events of magnitude>3,azimuthal gap<210°and number of phases>8.A nonlinear method using the oct-tree sampling algorithm from the NonLinLoc program was employed to determine the earthquake hypocenters.The hypocenter locations were then relocated using double-difference tomography(tomoDD).A significant reduction of travel-time(root mean square basis)and a better clustering of earthquakes were achieved which correlated well with the geological structure in West Java.Double-difference tomography was found to give a clear velocity structure,especially beneath the volcanic arc area,i.e.,under Mt Anak Krakatau,Mt Salak and the mountains complex in the southern part of West Java.Low velocity anomalies for the P and S waves as well as the vp/vs ratio below the volcanoes indicated possible partial melting of the upper mantle which ascended from the subducted slab beneath the volcanic arc.
基金Supported by the Major State Basic Research Development Program of China (973 Project) (Grant No.2006CB202306)the Applied Basic Project of China National Petroleum Corporation (Grant No.200510101)
文摘Based on Hudson’s theoretical hypothesis of equivalent fracture model,inserting aligned round chips in solid model can simulate fractured media. The effect of fractures on the propagation of P and S waves can be observed by changing the fracture thickness. The base model is made of epoxy resin,and the material of fractures is a kind of low-velocity mixture containing silicon rubber. With constant diameter and number of fractures in each model,one group of models can be formed through changing the thickness of fracture. These models have the same fracture density. By using the ultrasonic pulse transmission method,the experiment records time and waveform of P and S waves in the direction parallel and perpendicular to the fracture orientation. The result shows that,with the same fracture density,changing fracture aperture will affect both velocity and amplitude of P and S waves,and the effect on P-wave amplitude is much greater than that on the velocity. Moreover,the variation in velocity of S wave is more obvious in the slow shear wave (S2),while the variation in amplitude is more obvious in the fast shear wave (S1). These properties of wave propagation are useful for seismic data processing and interpretation.
基金supported by Funding Project for Academic Human Resources Development in Institutions of Higher Learning(Grant No.PHR201107145)
文摘The partial derivative equations of Zoeppritz equations are established and the derivatives of each matrix entry with respect to wave vectors are derived in this paper.By solving the partial derivative equations we obtained the partial derivatives of seismic wave reflection coefficients with respect to wave vectors,and computed the Goos-Hnchen shift for reflected P-and VS-waves.By plotting the curves of Goos-Hnchen shift,we gained some new insight into the lateral shift of seismic reflection wave.The lateral shifts are very large for glancing wave or the wave of the incidence angle near the critical angle,meaning that the seismic wave propagates a long distance along the reflection interface before returning to the first medium.For the reflection waves of incidence angles away from the critical angle,the lateral shift is in the same order of magnitude as the wavelength.The lateral shift varies significantly with different reflection interfaces.For example,the reflected P-wave has a negative shift at the reflection interface between mudstone and sandstone.The reflected VS-wave has a large lateral shift at or near the critical angle.The lateral shift of the reflected VS-wave tends to be zero when the incidence angle approaches 90°.These observations suggest that Goos-Hnchen effect has a great influence on the reflection wave of wide-angles.The correction for the error caused by Goos-Hnchen effect,therefore,should be made before seismic data processing,such as the depth migration and the normal-moveout correction.With the theoretical foundation established in this paper,we can further study the correction of Goos-Hnchen effect for the reflection wave of large incidence angle.
