A joint absorbing boundary for the multiple-relaxation-time lattice Boltzmann method in seismic acoustic wavefield modeling

Conventional seismic wave forward simulation generally uses mathematical means to solve the macroscopic wave equation,and then obtains the corresponding seismic wavefield.Usually,when the subsurface structure is finely constructed and the continuity of media is poor,this strategy is difficult to meet the requirements of accurate wavefield calculation.This paper uses the multiple-relaxation-time lattice Boltzmann method(MRT-LBM)to conduct the seismic acoustic wavefield simulation and verify its computational accuracy.To cope with the problem of severe reflections at the truncated boundaries,we analogize the viscous absorbing boundary and perfectly matched layer(PML)absorbing boundary based on the single-relaxation-time lattice Boltzmann(SRT-LB)equation to the MRT-LB equation,and further,propose a joint absorbing boundary through comparative analysis.We give the specific forms of the modified MRT-LB equation loaded with the joint absorbing boundary in the two-dimensional(2D)and three-dimensional(3D)cases,respectively.Then,we verify the effects of this absorbing boundary scheme on a 2D homogeneous model,2D modified British Petroleum(BP)gas-cloud model,and 3D homogeneous model,respectively.The results reveal that by comparing with the viscous absorbing boundary and PML absorbing boundary,the joint absorbing boundary has the best absorption performance,although it is a little bit complicated.Therefore,this joint absorbing boundary better solves the problem of truncated boundary reflections of MRT-LBM in simulating seismic acoustic wavefields,which is pivotal to its wide application in the field of exploration seismology.

FINITE ELEMENT ANALYSIS OF A VERTICAL RECTANGULAR PLATE COUPLED WITH AN UNBOUNDED FLUID DOMAIN ON ONE SIDE USING A TRUNCATED FAR BOUNDARY

The dynamic pressure distribution on a rectangular plate attached to a rigid wall and supporting an infinitely large extent of fluid subjected to a harmonic ground excitation is evaluated in the time domain. Governing equations for the fluid domain are set considering the compressibility of the fluid with negligibly small change in density and a linearized free surface. A far boundary condition for the three-dimensional fluid domain is developed so that the far boundary is truncated at a closer proximity to the structure. The coupled problem is solved independently for the structure and the fluid domain by transferring the acceleration of the plate to the fluid and pressure of the fluid to the plate in sequence. Helmholtz equation for the three-dimensional fluid domain and Mindlin＇s theory for the two-dimensional plate are used for the solution of the interacting domains. Finite element technique is adopted for the solution of this problem with pressure as nodal variable for the fluid domain and displacement for the plate. The time dependent equations are solved in each of the interacting domain using Newmark-fl method. The effectiveness of the technique is demonstrated and the influences of surface wave, exciting frequency and flexibility of the plate on dynamic pressure are investigated.