For the simulation of seismic wavefields in the ocean, it is very important to consider the influences of the irregular seabed and the rough sea surface. To numerically model seismic waves in the ocean, the fluid seawater and the solid seabed must be treated separately, and the air–fluid and fluid–solid boundary conditions must be considered simultaneously. A fluid–solid configuration is created with the influences of a rough sea surface and an irregular seabed to carry out wavefield simulation. Based on the boundary conformal grids, the 2D acoustic wave equations (in the fluid media), the 2D elastic wave equations (in the solid media) and two boundary conditions (at the interface between air and fluid and the interface between fluid and solid) are converted from the Cartesian coordinate system into a curvilinear coordinate system. In curvilinear coordinates, the seismic wave equations and boundary conditions are discretized with finite-difference operators. At the left, right and bottom boundaries of the model, the convolutional perfectly matched layer boundary condition is used to absorb artificial reflections. The numerical simulation results show that the snapshots and shot gathers of the rough sea surface model are more complex and contain stronger reflections than those of the flat model. The rough sea surface and the pattern of the irregular seabed have a great influence on the propagation of seismic waves in the ocean.

Amplitude variation with incidence angle (AVA) analysis is an essential tool for discriminating lithology in hydrocarbon reservoirs. Compared to traditional AVA inversion using only compressional wave (P-wave) information, joint AVA inversion using PP and PS seismic data provides better estimation of rock properties (e.g., density, P- and shear wave [S-wave] velocities). Currently, the most used AVA inversions depend on the approximations of the Zoeppritz equations (e.g., the Shuey and Aki-Richards approximations), which are not suitable for formations with strong contrast interfaces and seismic data with large incidence angles. Based on the previous derivation of the accurate Jacobian matrix, we have found that the sign of each partial derivative of reflection coefficient with respect to the P-, S-wave velocities and density changes across the interface represents a good indicator for the reflection interfaces. Accordingly, we adopt an adaptive stratified joint PP and PS AVA inversion using the accurate Jacobian matrix that can automatically obtain the layer information and can be further used as a constraint in the inversion of in-layer rock properties (density and P- and S-wave velocities). Due to the use of the exact Zoeppritz equations and accurate Jacobian matrix, this proposed inversion method is more accurate than traditional AVA inversion methods, has higher computational efficiency, and can be applied to seismic wide-angle reflection data or seismic data acquired for formations with strong contrast interfaces. The model study shows that this proposed inversion method works better than the classic Shuey and Aki-Richards approximations at estimating reflection interfaces and in-layer rock properties. It also works well in handling a part of the complex Marmousi 2 model and real seismic data.