Simultaneous-source algorithms can increase the efficiency of full waveform inversion (FWI) dramatically through reducing the number of times of wavefield simulations. However, the multiple sources ...induce the crosstalk artifacts, which severely contaminate the inversion results. To solve this problem, we provide a crosstalk-free simultaneous-source algorithm with the normalized seismic data. Sine harmonic functions with arbitrary phases are used as different wavelets in a super shot and as the encoding operator. Based on this algorithm, the crosstalk artifacts are eliminated by deblending the simultaneous-source wavefield with little additional computation. Moreover, the estimation or inversion of source wavelets at each iteration for simultaneous-source data, which is crucial for successful FWI but would severely reduce the efficiency of simultaneous-source algorithms, are avoided by normalizing the seismic data with deconvolution. Since the simultaneous-source data are deblended, the proposed algorithm is naturally applicable to the marine mobile streamer seismic data. Furthermore, it is convenient to select the reference traces for deblended data, instead of simultaneous-source data, to eliminate or unify the wavelet information by deconvolution. Finally, we verify the proposed algorithm with the synthetic data.
This article is concerned with the convergence properties of the Strang splitting method for the Degasperis-Procesi equation, which models shallow water dynamics. The challenges of analyzing ...splitting methods for this equation lie in the fact that the involved suboperators are both nonlinear. In this paper, instead of building the second order convergence in L2 for the proposed method directly, we first show that the Strang splitting method has first order convergence in H2. In the analysis, the Lie derivative bounds for the local errors are crucial. The obtained first order convergence result provides the H2 boundedness of the approximate solutions, thereby enabling us to subsequently establish the second order convergence in L2 for the Strang splitting method.
Celotno besedilo
Dostopno za:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, UILJ, UKNU, UL, UM, UPUK
Elastic full-waveform inversion (FWI) estimates the subsurface velocity model by inverting multi-component seismic data. Most of the current elastic FWI solves the elastic wave equation by a temporal ...second-order finite-difference time-domain (FDTD) method. However, the traditional temporal second-order FDTD introduces large temporal dispersion in the case of a large time-stepping size. To suppress the temporal dispersion, a relatively small time step should be used, which increases the computational cost of FWI. We introduce a temporal fourth-order and spatial arbitrary even-order FDTD method for elastic wave modeling and inversion. The temporal fourth-order scheme can simulate a highly accurate wavefield, even when a large time step is used for temporal extrapolation. Compared with the traditional temporal second-order FDTD scheme, the temporal fourth-order FDTD enables a larger time step for its more relaxed stability condition and smaller dispersion. As a result, FWI with temporal fourth-order FDTD scheme can achieve high accuracy inversion without reducing much computational efficiency. We run our elastic FWI codes on a Graphic Processing Unit (GPU) computing platform, and to accelerate the GPU computing, we propose to use the GPU shared memory to reduce the data access delay. By transferring the data frequently accessed for the FDTD approximations, from the global memory to the shared memory, we realize the FDTD calculation in a more efficient way. Forward modeling and FWI examples are carried out to confirm the effectiveness of our elastic FWI.
•Incorporating temporal fourth-order time-domain FD method into elastic FWI.•Developing a GPU shared memory access strategy to improve FWI efficiency.•Revealing the effects of temporal dispersion on FWI results.
A nonlinear dynamical model for the plankton population in a fixed sea area under the influence of asymmetric multiple factors, including atmospheric CO2 concentration, atmospheric temperature, ...nutrient concentration, seawater temperature, light intensity, and predator density is proposed to address the survival of the plankton population due to global warming. The model’s accuracy is confirmed by comparison with actual data, and numerical simulations are carried out to justify the relevant findings. The results suggest that increasing plankton’s ability to absorb atmospheric CO2 or regulate atmospheric temperature can help to mitigate global warming. Furthermore, if the population density of fish, the primary predator of plankton, falls within a certain range, the increase in atmospheric temperature will be mitigated. Additionally, the stability conditions for the suggested model are obtained, along with the equilibrium point of the system. Overall, this paper considers the effects of asymmetric multifactor interaction on plankton population density and establishes a mathematical connection between environmental ecosystems and plankton that might aid in addressing the challenges posed by global warming and preserving the plankton population.
Full-waveform inversion (FWI) serves as a useful tool to quantitatively investigate the properties of the media. One of the greatest challenges is its enormous computational cost, especially for the ...three-dimensional (3D) case. The multisource encoding strategy is a crucial method to speed up FWI. However, current 3D elastic FWI algorithms with the simultaneous sources, such as random time series and random polarity methods, may be subjected to the crosstalk noise, which severely degrades inversion quality. We present a novel 3D elastic crosstalk-free multisource full-waveform inversion method (MFWI). The encoded sources with a series of single-frequency components are excited to perform the seismic wavefield simulation. Extracting the wavefield corresponding to each single-frequency harmonic source by phase-sensitive detection (PSD) algorithm, the multisource wavefield can be decomposed. Therefore, summing the gradients with respect to the model parameters of all decomposed single frequencies can effectively avoid the crosstalk in our MFWI algorithm. Furthermore, to increase the practicability and universality of the algorithm, we adopt the graphic processing unit (GPU) based on the domain decomposition to improve the MFWI algorithm. The multi-scale strategy with a flexible encoding number in each scale helps to enhance the robustness of our algorithm. The numerical examples demonstrate that our method can achieve domain decomposition simulation and obtain 3D crosstalk-free multisource elastic wave inversion results.
•Proposed a crosstalk-free multisource 3D elastic FWI.•Develop an efficient FWI algorithm suitable for distributed GPU clusters.•Incorporate a GPUs-based domain decomposition scheme into the proposed method.
High-precision elastic reverse-time migration (ERTM) imaging has always been one of the trends in the development of geophysics. However, current wavefield simulations using time-domain ...finite-difference (FD) approaches in ERTM have second-order temporal accuracy, resulting in travel time changes and waveform distortion in wavefield propagation with large time steps, i.e., temporal dispersion. Errors caused by the temporal dispersion can lead to erroneous imaging locations and out-of-focus diffraction events. A new ERTM and its workflow are established here using temporal and spatial high-order FD accuracy wavefields and the vector-based imaging condition. Our method computes elastic vector-based wavefields by solving a P- and S-wave decomposition form of a quasi-stress–velocity equation. An advanced finite-difference scheme is employed in the wavefield solution to achieve simulation with temporal fourth-order accuracy and spatial arbitrary even-order accuracy. The normalized dot-product imaging condition of the source and receiver P/S wavefields is then applied to generate high-quality images. The elastic wavefield simulation and ERTM numerical examples presented here reveal that the anti-dispersion workflow can improve modeling and imaging accuracy. In addition, the field data application shows that our method can achieve reasonable and reliable ERTM images. This method can integrate the most advanced imaging techniques into this computational framework to improve imaging accuracy.
Hydrodynamic lubrication is widely used between two relatively moving objects, and the effect of fluid flow state and temperature distribution on lubrication performance in wedge-shaped gaps is a ...popular topic to study. In this paper, the incompressible double-distribution lattice Boltzmann method (LBM) is applied to study the effect of micro convex surface texture on micro lubrication and heat transfer in wedge-shaped channels. By comparing this model with the analytical solution of an infinitely wide wedge slider, the maximum pressure calculated by LBM is 0.1081 MPa, and the maximum pressure calculated by the Reynolds equation is 0.1079 MPa. The error of the maximum pressure is 1.11%, and the Reynolds equation result is slightly smaller. The reason is that the Reynolds equation ignores the influence of fluid inertia force on oil film pressure. The results indicate that the application of LBM can be used to study lubrication problems. Compared with the Reynolds equation, LBM can calculate the velocity field and pressure field in the film thickness direction, and can also observe precise flow field details such as vortices. Three micro convex texture shapes were established to study the effects of different convex textures on micro lubrication and oil film temperature distribution, and the velocity distribution, temperature distribution and oil film pressure along the oil film thickness direction were given. Under the same conditions, comparing the oil film pressure with and without surface texture, the results show that the maximum oil film pressure with surface texture 3 is increased by about 4.34% compared with that without surface texture. The slightly convex texture can increase the hydrodynamic lubrication effect and obtain greater load-bearing capacity, helping to reduce the possibility of contact friction. The results show that the convex surface texture can improve the hydrodynamic lubrication performance, increase the load carrying capacity and reduce the possibility of contact friction, and the convex surface texture can influence the temperature distribution of the oil film. At 3.6 mm in the slider length direction and 7.5 μm in the oil film thickness direction, the temperature of surface texture 1 is 402.64 K, the temperature of surface texture 2 is 403.31 K, and the temperature of surface texture 3 is 403.99 K. The presence of vortices is captured at a high convergence ratio.
Full waveform inversion (FWI) is regarded as an effective method to obtain high accuracy subsurface parameters like velocity and density. When it is applied to real seismic data, some challenges may ...appear. Real seismic data contain complex wavefields including compressional waves, shear waves and their interconversions, however, acoustic FWI only matches compressional waves but not including the converted compressional waves and treats others as noise. Though elastic FWI can take all types of wave modes into consideration, the inversion results will be compromised by the lack of multi-component input since only the vertical component is generally acquired in practice. It is interesting to investigate the effects of a single-component input generated by an elastic model on acoustic and elastic FWI. The investigation is carried out using a modified Marmousi2 model. Numerical examples illustrate that elastic FWI using a single vertical component observation data has an advantage in obtaining higher resolution inversion results than acoustic FWI. However, the ability of elastic FWI to recover more accurate velocities, especially for S-wave velocity, is limited to a degree if the horizontal component is missing. This situation will worsen in the presence of noise. It is also found that the inversion results of acoustic FWI from the elastic data generated by an explosive source are acceptable, whereas those of acoustic FWI from the elastic data created by a vertical-directional normal stress source are of very low accuracy.
•A convolutional objective function is used to acoustic and elastic FWI.•Reveal the effects of the elastic data excited by two typesof sources on elastic single component FWI and acoustic FWI.•Following the point 2, we also test effects on FWI when the observed data are corrupted by noise.
Three-dimensional (3D) elastic reverse-time migration (ERTM) can image the subsurface 3D seismic structures, and it is an important tool for the Earth’s interior imaging. A common simulation kernel ...used in 3D ERTM is the current staggered-grid finite-difference (SGFD) method of the first-order elastic wave equation. However, the mere second-order accuracy in time of the current SGFD method can bring non-negligible time dispersion, which reduces the simulation accuracy and further leads to the distortion of the imaging results. This paper proposes a vector-based 3D ERTM using the high-order accuracy SGFD method in time to obtain high-accuracy images. This approach is a new high-resolution ERTM workflow that improves the imaging accuracy of conventional ERTM from numerical simulation. The proposed ERTM workflow is established on a quasi-stress–velocity wave equation and its vector wavefield decomposition form. Advanced SGFD schemes and their corresponding coefficients with fourth-order temporal accuracy solve the quasi-linear wave equation system. The normalized dot product imaging condition produces high-quality images using high-accuracy vector wavefields solved using the SGFD method. Through the numerical examples, we test the simulation efficiency and analyze how temporal accuracy in numerical simulations affects migration imaging quality. We include that the proposed method obtains highly accurate images.
Delay differential equations have been used to model numerous phenomena in nature. We extend the previous work of one of the authors to analyze the stability properties of the explicit exponential ...Rosenbrock methods for stiff differential equations with constant delay. We first derive sufficient conditions so that the exponential Rosenbrock methods satisfy the desired stability property. We accomplish this without relying on some extreme constraints, which are usually necessary in stability analysis. Then, with the aid of the integral form of the method coefficients, we provide a simple stability criterion that can be easily verified. We also present a theorem on the order barrier for the proposed methods, stating that there is no method of order five or higher that satisfies the simple criterion. Numerical tests are carried out to validate the theoretical results.
•Stable explicit exponential Rosenbrock methods for stiff delay differential equations are derived without extreme constraints.•A simple stability criterion is presented using integral representation of method coefficients.•Order barrier theorem: No method of order five or higher satisfies the stability criterion.