Reflection coefficients of P-SV waves in weak anisotropic media

We introduce the Thomsen anisotropic parameters into the approximate linear reflection coefficient equation for P-SV wave in weakly anisotropic HTI media. From this we get a new, more effective, and practical reflection coefficient equation. We performed forward modeling to AVO attributes, obtaining excellent results. The combined AVO attribute analysis of PP and PS reflection data can greatly reduce ambiguity, obtain better petrophysical parameters, and improve parameter accuracy.

Scattering of Water Waves by Dual Symmetric Inclined Floating Porous Barriers Using the DBEM

The scattering of normally incident water waves by two surface-piercing inclined perforated barriers in water with a uniform finite depth is investigated within the framework of linear water wave theory.Considering that thin barriers are zero-thickness,a novel numerical method involving the the coupling of the dual boundary element method(DBEM)with damping layers is applied.In order to effectively damp out the reflected waves,two damping layers,instead of pseudoboundaries are implemented near the two side boundaries of the computational domain.Thus,the modified linearized free surface boundary conditions are formulated and used for solving both the ordinary boundary integral equation as well as the hypersingular boundary integral equation for degenerate boundaries.The newly developed numerical method is validated against analytical methods using the matched eigenfunction expansion method for the special case of two vertical barriers or the inclined angle to the vertical being zero.The influence of the length of the two damping layers has been discussed.Moreover,these findings are also validated against previous results for several cases.After validation,the numerical results for the reflection coefficient,transmission coefficient and dissipation coefficient are obtained by varying the inclination angle and porosity-effect parameter.The effects of both the inclination angle and the porosity on the amplitudes of wave forces acting on both the front and rear barriers are also investigated.It is found that the effect of the inclination angle mainly shifts the location of the extremal values of the reflection and the transmission coefficients.Additionally,a moderate value of the porosity-parameter is quite effective at dissipating wave energy and mitigating the wave loads on dual barriers.

不均匀海底地形影响下柔性板与斜向来浪的相互作用

The present work analyzes the interaction of oblique waves by a porous flexible breakwater in the presence of a step-type bottom.The physical models for both scattering and trapping cases are considered and developed within the framework of small amplitude water-wave theory.Darcy’s law is used to model the wave interaction with the porous medium.It is assumed that the varying bottom extends over a finite interval,connected by a finite length of uniform bottom near an impermeable wall,and a semi-infinite length of bottom in the open water region.The boundary value problem is solved using the eigenfunction expansion method in the uniform bottom regions,while a modified mild-slope equation(MMSE)is used for the region with the varying bottom.Additionally,a mass-conserving jump condition is employed to handle the solution at slope discontinuities in the bottom.A system of equations is derived by matching the solutions at interfaces.The reflection coefficient and force on the breakwater and impermeable wall are plotted and analyzed for various parameters,such as the length of the varying bottom,depth ratio,angle of incidence,and flexural rigidity.It is observed that moderate values of flexural rigidity and depth ratio significantly contribute to an optimum reflection coefficient and reduce the wave force on the wall and breakwater.Remarkably,the outcomes of this study are assumed to be applicable in the construction of this type of breakwater in coastal regions.

Wave Reflection by Rectangular Breakwaters for Coastal Protection

In this study,we focus on the numerical modelling of the interaction between waves and submerged structures in the presence of a uniform flow current.Both the same and opposite senses of wave propagation are considered.The main objective is an understanding of the effect of the current and various geometrical parameters on the reflection coefficient.The wave used in the study is based on potential theory,and the submerged structures consist of two rectangular breakwaters positioned at a fixed distance from each other and attached to the bottom of a wave flume.The numerical modeling approach employed in this work relies on the Boundary Element Method(BEM).The results are compared with experimental data to validate the approach.The findings of the study demonstrate that the double rectangular breakwater configuration exhibits superior wave attenuation abilities if compared to a single rectangular breakwater,particularly at low wavenumbers.Furthermore,the study reveals that wave mitigation is more pronounced when the current and wave propagation are coplanar,whereas it is less effective in the case of opposing current.

On the eigenvalues and eigendisplacement of the critical mode in horizontally layered media

Wave propagation in horizontally layered media is a classical problem in seismic-wave theory.In semi-infinite space,a nondispersive Rayleigh wave mode exists,and the eigendisplacement decays exponentially with depth.In a layered model with increasing layer velocity,the phase velocity of the Rayleigh wave varies between the S-wave velocity of the bottom half-space and that of the classical Rayleigh wave propagated in a supposed half-space formed by the parameters of the top layer.If the phase velocity is the same as the P-or S-wave velocity of the layer,which is called the critical mode or critical phase velocity of surface waves,the general solution of the wave equation is not a homogeneous(expressed by trigonometric functions)or inhomogeneous(expressed by exponential functions)plane wave,but one whose amplitude changes linearly with depth(expressed by a linear function).Theories based on a general solution containing only trigonometric or exponential functions do not apply to the critical mode,owing to the singularity at the critical phase velocity.In this study,based on the classical framework of generalized reflection and transmission coefficients,the propagation of surface waves in horizontally layered media was studied by introducing a solution for the linear function at the critical phase velocity.Therefore,the eigenvalues and eigenfunctions of the critical mode can be calculated by solving a singular problem.The eigendisplacement characteristics associated with the critical phase velocity were investigated for different layered models.In contrast to the normal mode,the eigendisplacement associated with the critical phase velocity exhibits different characteristics.If the phase velocity is equal to the S-wave velocity in the bottom half-space,the eigendisplacement remains constant with increasing depth.

孔隙率对矩形沟槽上垂直多孔屏障的波浪散射效应影响

The effect of porosity on surface wave scattering by a vertical porous barrier over a rectangular trench is studied here under the assumption of linearized theory of water waves.The fluid region is divided into four subregions depending on the position of the barrier and the trench.Using the Havelock’s expansion of water wave potential in different regions along with suitable matching conditions at the interface of different regions,the problem is formulated in terms of three integral equations.Considering the edge conditions at the submerged end of the barrier and at the edges of the trench,these integral equations are solved using multi-term Galerkin approximation technique taking orthogonal Chebyshev’s polynomials and ultra-spherical Gegenbauer polynomial as its basis function and also simple polynomial as basis function.Using the solutions of the integral equations,the reflection coefficient,transmission coefficient,energy dissipation coefficient and horizontal wave force are determined and depicted graphically.It was observed that the rate of convergence of the Galerkin method in computing the reflection coefficient,considering special functions as basis function is more than the simple polynomial as basis function.The change of porous parameter of the barrier and variation of trench width and height significantly contribute to the change in the scattering coefficients and the hydrodynamic force.The present results are likely to play a crucial role in the analysis of surface wave propagation in oceans involving porous barrier over submarine trench.

Seismic reflection and transmission coefficients of a single layer sandwiched between two dissimilar poroelastic solids

The seismic reflection and transmission characteristics of a single layer sandwiched between two dissimilar poroelastic solids saturated with two immiscible viscous fluids are investigated. The sandwiched layer is modeled as a porous solid with finite thickness. The propagation of waves is represented with potential functions. The displacements of particles in different phases of the aggregate are defined in terms of these potential functions. Due to the presence of viscosity in pore fluids, the reflected and transmitted waves are inhomogeneous in nature, i.e., with different directions of propagation and attenuation. The closed-form analytical expressions for reflection and transmission coefficients are derived theoretically for appropriate boundary conditions. These expressions are calculated as a non-singular system of linear algebraic equations and depend on the various parameters involved in this non-singular system. Hence,numerical examples are studied to determine the effects of various properties of the sandwich layer on reflection and transmission coefficients. The essential features of layer thickness, incident direction, wave frequency, liquidsaturation and capillary pressure of the porous layer on reflection and transmission coefficients are depicted graphically and discussed. The analysis shows that reflection and transmission coefficients are strongly associated with incident direction and various properties of the porous layer.

Analysis of Directional Spectra and Reflection Coefficients in Incident and Reflected Wave Field

In this paper, the modified Bayesian method for the analysis of directional wave spectra and reflection coefficients is verified by numerical and physical simulation of waves. The results show that the method can basically separate the incident and reflected directional spectra. In addition, the effect of the type of wave gage arrays, the number of measured wave properties, and the distance between the wave gage array and the reflection line on the resolution of the method are investigated. Some suggestions are proposed for practical application.

Monitoring the change in horizontal stress with multi-wave time-lapse seismic response based on nonlinear elasticity theory

Monitoring the change in horizontal stress from the geophysical data is a tough challenge, and it has a crucial impact on broad practical scenarios which involve reservoir exploration and development, carbon dioxide (CO_(2)) injection and storage, shallow surface prospecting and deep-earth structure description. The change in in-situ stress induced by hydrocarbon production and localized tectonic movements causes the changes in rock mechanic properties (e.g. wave velocities, density and anisotropy) and further causes the changes in seismic amplitudes, phases and travel times. In this study, the nonlinear elasticity theory that regards the rock skeleton (solid phase) and pore fluid as an effective whole is used to characterize the effect of horizontal principal stress on rock overall elastic properties and the stress-dependent anisotropy parameters are therefore formulated. Then the approximate P-wave, SV-wave and SH-wave angle-dependent reflection coefficient equations for the horizontal-stress-induced anisotropic media are proposed. It is shown that, on the different reflectors, the stress-induced relative changes in reflectivities (i.e., relative difference) of elastic parameters (i.e., P- and S-wave velocities and density) are much less than the changes in contrasts of anisotropy parameters. Therefore, the effects of stress change on the reflectivities of three elastic parameters are reasonably neglected to further propose an AVO inversion approach incorporating P-, SH- and SV-wave information to estimate the change in horizontal principal stress from the corresponding time-lapse seismic data. Compared with the existing methods, our method eliminates the need for man-made rock-physical or fitting parameters, providing more stable predictive power. 1D test illustrates that the estimated result from time-lapse P-wave reflection data shows the most reasonable agreement with the real model, while the estimated result from SH-wave reflection data shows the largest bias. 2D test illustrates the feasibility of the proposed inversion method for estimating the change in horizontal stress from P-wave time-lapse seismic data.

Deghosting towed streamer data in τ/p domain based on rough sea surface reflectivity

Currently, the deghosting of towed streamer seismic data assumes a flat sea level and a sea-surface reflection coefficient of-1; this decreases the precision of deghosting. A new method that considers the rough sea surface is proposed to suppress ghost reflections. The proposed deghosting method obtains the rough sea surface reflection coefficient using Gaussian statistics, and calculates the optimized deghosting operator in the r/p domain. The proposed method is closer to the actual sea conditions, offers an improved deghosting operator, removes the ghost reflections from marine towed seismic data, widens the bandwidth and restores the low-frequency information, and finally improves the signal-to- noise ratio and resolution of the seismic data.

Estimating primaries by sparse inversion of the 3D Curvelet transform and the L1-norm constraint

In this paper, we built upon the estimating primaries by sparse inversion （EPSI） method. We use the 3D curvelet transform and modify the EPSI method to the sparse inversion of the biconvex optimization and Ll-norm regularization, and use alternating optimization to directly estimate the primary reflection coefficients and source wavelet. The 3D curvelet transform is used as a sparseness constraint when inverting the primary reflection coefficients, which results in avoiding the prediction subtraction process in the surface-related multiples elimination （SRME） method. The proposed method not only reduces the damage to the effective waves but also improves the elimination of multiples. It is also a wave equation- based method for elimination of surface multiple reflections, which effectively removes surface multiples under complex submarine conditions.

Application of BDRM Theory in the Two-AXIS Underwater Channel

Two-axis underwater channel often exists in deep ocean. Because of the coupling between surface channel and SOFAR channel, sound propagation in the two-axis underwater channel is complex and so its calculations of acoustic fields are difficult. The beam-displacement ray-mode （BDRM） theory is a normal mode method for propagation modeling in the common horizontally stratified shallow water. We improve the theory, proposing a new method for computing the upper boundary reflection coefficient, and apply it to calculate the acoustic fields of two-axis underwater channel. Transmission losses in the two-axis underwater channel are calculated by the BDRM theory. The results are in good agreement with the KRAKEN code and the computational speed excels those of the other methods.

AVO approximation for PS-wave and its application in PP/PS joint inversion

Multi-component exploration has many advantages over ordinary P-wave exploration. PP/PS joint AVO analysis and inversion are useful and powerful methods to discriminate between reservoir and non-productive lithology. In this paper, we derive a new PS-wave reflection coefficient approximation equation which is more accurate at larger incidence angles. The equation is simplified for small incidence angles, which makes AVO analysis clearer and easier for angles less than 30 degrees. Based on this approximation, a PP/PS joint inversion is introduced. A real data example shows that oil sands, brine sands and shales can be differentiated based on the P- to S-wave velocity ratio from the PP/PS joint inversion. Fluid factors and Poisson＇s ratio also indicate an anomaly in the target zone at the oil well location.

Reflection and transmission of bottom simulating reflectors in gas hydrate-bearing sediments: Two-phase media models

The bottom simulating reflector （BSR） in gas hydrate-bearing sediments is a physical interface which is composed of solid, gas, and liquid and is influenced by temperature and pressure. Deep sea floor sediment is a porous, unconsolidated, fluid saturated media. Therefore, the reflection and transmission coefficients computed by the Zoeppritz equation based on elastic media do not match reality. In this paper, a two-phase media model is applied to study the reflection and transmission at the bottom simulating reflector in order to find an accurate wave propagation energy distribution and the relationship between reflection and transmission and fluid saturation on the BSR. The numerical experiments show that the type I compressional （fast） and shear waves are not sensitive to frequency variation and the velocities change slowly over the whole frequency range. However, type II compressional （slow） waves are more sensitive to frequency variation and the velocities change over a large range. We find that reflection and transmission coefficients change with the amount of hydrate and free gas. Frequency, pore fluid saturation, and incident angle have different impacts on the reflection and transmission coefficients. We can use these characteristics to estimate gas hydrate saturation or detect lithological variations in the gas hydrate-bearing sediments.

AVO inversion based on common shot migration

Conventional AVO inversion utilizes the trace amplitudes of CMP gathers. There are three main factors affecting the accuracy of the inversion. First, CMP gathers are based on the hypothesis of horizontal layers but most real layers are not horizontal. Greater layer dip results in a greater difference between the observed CMP gathers and their real location. Second, conventional processing flows such as NMO, DMO, and deconvolution will distort amplitudes. Third, the formulation of reflection coefficient is related to incidence angles and it is difficult to get the relationship between amplitude and incidence angle. Wave equation prestack depth migration has the ability of imaging complex media and steeply dipping layers. It can reduce the errors of conventional processing and move amplitudes back to their real location. With true amplitude migration, common angle gathers abstraction, and AVO inversion, we suggest a method of AVO inversion from common shot gathers in order to reduce the effect of the above factors and improve the accuracy of AVO inversion.

Improving seismic interpretation: a high-contrast approximation to the reflection coefficient of a plane longitudinal wave

Linearized approximations of reflection and transmission coefficients set a foundation for amplitude versus offset（AVO） analysis and inversion in exploration geophysics.However,the weak properties contrast hypothesis of those linearized approximate equations leads to big errors when the two media across the interface vary dramatically.To extend the application of AVO analysis and inversion to high contrast between the properties of the two layers,we derive a novel nonlinearized high-contrast approximation of the PP-wave reflection coefficient,which establishes the direct relationship between PPwave reflection coefficient and P-wave velocities,S-wave velocities and densities across the interface.（A PP wave is a reflected compressional wave from an incident compressional wave（P-wave）.） This novel approximation is derived from the exact reflection coefficient equation with Taylor expansion for the incident angle.Model tests demonstrate that,compared with the reflection coefficients of the linearized approximations,the reflection coefficients of the novel nonlinearized approximate equation agree with those of the exact PP equation better for a high contrast interface with a moderate incident angle.Furthermore,we introduce a nonlinear direct inversion method utilizing the novel reflection coefficient equation as forward solver,to implement the direct inversion for the six parameters including P-wave velocities,S-wave velocities,and densities in the upper and lower layers across the interface.This nonlinear inversion algorithm is able to estimate the inverse of the nonlinear function in terms of model parameters directly rather than in a conventional optimization way.Three examples verified the feasibility and suitability of this novel approximation for a high contrast interface,and we still could estimate the six parameters across the interface reasonably when the parameters in both media across the interface vary about 50%.

Ultrasonic characterization of modified Cr_2O_3 coatings by reflection coefficient spectroscopy

Pores,microcracks and density of plasma sprayed Cr2O3 coatings before and after high-intensity pulsed ion beam(HIPIB) irradiation were investigated using the ultrasonic reflection coefficient spectroscopy(URCS).The URCS was analyzed based on an acoustic transmission model for the multi-layered structure.The longitudinal velocity in the coatings was calculated from the experimental URCS,and the attenuation coefficient expression was deduced by comparing the experimental and numerical fitting amplitude spectral lines.The longitudinal velocity of as-sprayed Cr2O3 coating is 2 002 m/s,and increases to 2 099 and 2 148 m/s after being irradiated by HIPIB with 1 and 5 shots.Correspondingly,the factor A changes from 0.046 to 0.026 and 0.020 and n from 1.702 to 1.658 and 1.649 in the attenuation coefficient expression of α=Af n.It is observed that the surface morphology of Cr2O3 coatings changes from rough and porous to smooth and uniform with the increase of shot number,which accords with the ultrasonic analyses reasonably.The URCS seems to provide a convenient and nondestructive method to characterize surface modification of the plasma sprayed coatings.

Penetration characteristics of VLF wave from atmosphere into lower ionosphere

The factors affecting the reflection and transmission coefficient of the ionosphere have been analyzed.These factors include wave frequency,incident angle,geomagnetic inclination,electron density and collision frequency in the ionosphere.The ionosphere refractive index is also analyzed.The ionosphere above 70 km is considered to be homogeneous and anisotropic,and the reflection and transmission coefficient matrix is calculated using matrix method.Simultaneously the Booker quartic equation is solved to get the refractive index in the ionosphere.The results show that when the wave frequency is higher,it is easier to penetrate into the ionosphere from its bottom boundary and the propagation attenuation in the ionosphere is smaller.TE（traverse electric） wave and TM（traverse magnetic） wave can both penetrate into the ionosphere with a small incident angle,while TE wave can hardly transmit into the ionosphere when the incident angle is large.The transmission coefficient decreases as the geomagnetic inclination increases.TE and TM wave cannot penetrate into the ionosphere at magnetic equator.When the electron collision frequency is higher,it is easier for VLF wave to penetrate into the ionosphere and the attenuation of ordinary wave is weaker,which may be caused by the energy transportation between the waves and the particles.The ordinary（O） wave experiences severer attenuation than extraordinary（X） wave,and X wave is a penetration mode whereas O wave is a non-penetration mode in the ionosphere.All the results indicate that VLF wave with higher frequency is easier to penetrate into the ionosphere and to be recorded by the satellites at high latitude.It is hard for ULF and the lower frequency VLF wave to transmit into the ionosphere directly for the severe reflection and attenuation.It may transmit into the ionosphere with a small incident angle due to the nonlinear effect,for example,the interaction between the waves and the particles or cross modulation,and then propagate along the whistle duct with small attenuation.This work may be a preliminary theoretical exploration for the future calculation on the response of ground based VLF artificial transmitter in the ionosphere and further study on the seismic ionosphere coupling model.

Analysis of the ambiguity of log-constrained seismic impedance inversion

Although the ambiguity of seismic inversion is widely recognized in both theory and practice, so far as a concrete inversion example is concerned, there is not any objective, controllable method or any standard for how to evaluate and determine its ambiguity and reliability, especially for the high frequency components beyond the effective seismic frequency band. Taking log-constrained impedance inversion as an example, a new appraisal method is proposed on the basis of analyzing a simple geological model. Firstly, the inverted impedance model is transformed to a reflection coefficient series. Secondly, the maximum effective frequency of the real seismic data is chosen as a cutoff point and the reflection coefficient series is decomposed into two components by low-pass and high-pass filters. Thirdly, the geometrical reflection characteristics of the high-frequency components and that of the real seismic data are compared and analyzed. Then, the reliability of the inverted impedance model is appraised according to the similarity of geometrical characteristics between the high-frequency components and the real seismic data. The new method avoids some subjectivity in appraising the inverted result, and helps to enhance the reliability of reservoir prediction by impedance inversion technology.

Hilbert Transferin Applied to Separation of Waves

The analytical method (AM) for separation of composite waves is presented based on the Hilbert transform. It is applicable to both regular and irregular trains; of waves. The wave data series measured with two wave gauges in the experiments are separated into two series of incident and reflected waves. Then, the reflection coefficient can be easily obtained. The arrival of reflected waves can also be detected for improvement of the accuracy of the reflection coefficient. The reflection performance of the physical model can be estimated exactly without calculation of wave height and phase difference. Numerical samples developed to test the method are proved to be accurate. Physical experiments are conducted and compared with Goda's method and satisfactory results are obtained.