Optimized finite-difference operator for broadband seismic wave modeling

Geophysics ◽  
2013 ◽  
Vol 78 (1) ◽  
pp. A13-A18 ◽  
Author(s):  
Jin-Hai Zhang ◽  
Zhen-Xing Yao

High-resolution image and waveform inversion of small-scale targets requires the handling of high-frequency seismic wavefields. However, conventional finite-difference (FD) methods have strong numerical dispersions in the presence of high-frequency components. To reduce these numerical dispersions, we optimized the constant coefficients of the FD operator by maximizing the wavenumber coverage within a given error limitation. We set up three general criteria to enhance the convergence of the algorithm and reduce the optimization effort. We selected the error limitation to be 0.0001, this being the smallest in the literature, which led to perfect agreement between theoretical analyses and numerical experiments. The accuracy of our optimized FD methods can even reach that of much higher order unoptimized FD methods, which means great savings of computational efforts and memory demand. These advantages become even more apparent with 3D modeling, especially for saving memory demand.

2016 ◽  
Vol 2016 ◽  
pp. 1-15 ◽  
Author(s):  
Manuel Kindelan ◽  
Miguel Moscoso ◽  
Pedro Gonzalez-Rodriguez

We present a method to obtain optimal finite difference formulas which maximize their frequency range of validity. The optimization is based on the idea of keeping an error of interest (dispersion, phase, or group velocities errors) below a given threshold for a wavenumber interval as large as possible. To find the weights of these optimal finite difference formulas we solve a system of nonlinear equations. Furthermore, we give compact formulas for the optimal weights as function of the error bound. Several numerical experiments illustrate the performance of the obtained finite difference formulas compared to the standard ones.


1991 ◽  
Vol 81 (6) ◽  
pp. 2351-2370
Author(s):  
Zoltan A. Der ◽  
Robert H. Shumway ◽  
Michael R. Hirano

Abstract Waveform modeling in the time domain matches the various frequency components of seismic signals unevenly; the agreement is better at low frequencies and becomes progressively worse towards higher frequencies. The net effect of this kind of time-domain modeling is that the resolution in the spatial details of the source is less than optimal since the high-frequency components of the signal with their short wavelengths to resolve finer details do not fit the data. These problems are demonstrated by numerical simulations and by the reanalysis of some aspects of the El Golfo earthquake in using a new seismic imaging technique based on a generalization of an f-k algorithm. This procedure computes a statistic that can be used to derive confidence limits of the parameters sought in the inversion, thus providing a quantitative measure of the uncertainties in the results.


2018 ◽  
Vol 26 (02) ◽  
pp. 1850005 ◽  
Author(s):  
Jian Wang ◽  
Xiaohong Meng ◽  
Hong Liu ◽  
Wanqiu Zheng ◽  
Zhiwei Liu

Staggered-grid finite-difference forward modeling in the time domain has been widely used in reverse time migration and full waveform inversion because of its low memory cost and ease to implementation on GPU, however, high dominant frequency of wavelet and big grid interval could result in significant numerical dispersion. To suppress numerical dispersion, in this paper, we first derive a new weighted binomial window function (WBWF) for staggered-grid finite-difference, and two new parameters are included in this new window function. Then we analyze different characteristics of the main and side lobes of the amplitude response under different parameters and accuracy of the numerical solution between the WBWF method and some other optimum methods which denotes our new method can drive a better finite difference operator. Finally, we perform elastic wave numerical forward modeling which denotes that our method is more efficient than other optimum methods without extra computing costs.


1983 ◽  
Vol 55 (4) ◽  
pp. 1120-1127 ◽  
Author(s):  
G. Charbonneau ◽  
J. L. Racineux ◽  
M. Sudraud ◽  
E. Tuchais

An accurate recording system was set up and used for analyzing normal and asthmatic breath-sound features. Breath sounds are recorded at the trachea simultaneously with the airflow signal at 0.5- and 1-1/s levels. The study was carried out in the frequency domain using a fast-fourier transform (FFT). FFTs are taken on 1,024-sample blocks (one block = 200 ms) over a duration of about 20 s. Different characteristics of the spectra are calculated in the range 60-1,260 Hz for 11 normal and 10 asthmatic subjects. This allows the computation of an index that discriminates (P less than 0.0005) asthma cases from normal cases. Spectral features strongly depend on the flow rate both for normal and asthmatic subjects. Increasing the flow rate raises the high-frequency components of the spectra.


Author(s):  
G. Y. Fan ◽  
J. M. Cowley

It is well known that the structure information on the specimen is not always faithfully transferred through the electron microscope. Firstly, the spatial frequency spectrum is modulated by the transfer function (TF) at the focal plane. Secondly, the spectrum suffers high frequency cut-off by the aperture (or effectively damping terms such as chromatic aberration). While these do not have essential effect on imaging crystal periodicity as long as the low order Bragg spots are inside the aperture, although the contrast may be reversed, they may change the appearance of images of amorphous materials completely. Because the spectrum of amorphous materials is continuous, modulation of it emphasizes some components while weakening others. Especially the cut-off of high frequency components, which contribute to amorphous image just as strongly as low frequency components can have a fundamental effect. This can be illustrated through computer simulation. Imaging of a whitenoise object with an electron microscope without TF limitation gives Fig. 1a, which is obtained by Fourier transformation of a constant amplitude combined with random phases generated by computer.


2019 ◽  
Vol 14 (7) ◽  
pp. 658-666
Author(s):  
Kai-jian Xia ◽  
Jian-qiang Wang ◽  
Jian Cai

Background: Lung cancer is one of the common malignant tumors. The successful diagnosis of lung cancer depends on the accuracy of the image obtained from medical imaging modalities. Objective: The fusion of CT and PET is combining the complimentary and redundant information both images and can increase the ease of perception. Since the existing fusion method sare not perfect enough, and the fusion effect remains to be improved, the paper proposes a novel method called adaptive PET/CT fusion for lung cancer in Piella framework. Methods: This algorithm firstly adopted the DTCWT to decompose the PET and CT images into different components, respectively. In accordance with the characteristics of low-frequency and high-frequency components and the features of PET and CT image, 5 membership functions are used as a combination method so as to determine the fusion weight for low-frequency components. In order to fuse different high-frequency components, we select the energy difference of decomposition coefficients as the match measure, and the local energy as the activity measure; in addition, the decision factor is also determined for the high-frequency components. Results: The proposed method is compared with some of the pixel-level spatial domain image fusion algorithms. The experimental results show that our proposed algorithm is feasible and effective. Conclusion: Our proposed algorithm can better retain and protrude the lesions edge information and the texture information of lesions in the image fusion.


Author(s):  
Priya R. Kamath ◽  
Kedarnath Senapati ◽  
P. Jidesh

Speckles are inherent to SAR. They hide and undermine several relevant information contained in the SAR images. In this paper, a despeckling algorithm using the shrinkage of two-dimensional discrete orthonormal S-transform (2D-DOST) coefficients in the transform domain along with shock filter is proposed. Also, an attempt has been made as a post-processing step to preserve the edges and other details while removing the speckle. The proposed strategy involves decomposing the SAR image into low and high-frequency components and processing them separately. A shock filter is used to smooth out the small variations in low-frequency components, and the high-frequency components are treated with a shrinkage of 2D-DOST coefficients. The edges, for enhancement, are detected using a ratio-based edge detection algorithm. The proposed method is tested, verified, and compared with some well-known models on C-band and X-band SAR images. A detailed experimental analysis is illustrated.


Author(s):  
E. Faghand ◽  
S. Karimian ◽  
E. Mehrshahi ◽  
N. Karimian

Abstract A new simple computational tool is proposed for the synthesis of multi-section coupled-line filters based on combined modified planar circuit method (MPCM) and transmission line method (TLM) analysis, referred to as MPCM-TLM. Due to its fundamentally simple architecture, the presented tool offers significantly faster optimization of coupled-line filters – for exactly the same initial simulation set-up – than other costly commercially-available tools, giving equally reliable results. Validity and accuracy of the proposed tool have been verified through the design of 3rd, 5th, and 7th order coupled-line filters and comparative analysis between results obtained from the proposed approach and the high-frequency structure simulator. A remarkable 99% time reduction in the analysis is recorded in the case of 7th order filter using the proposed tool, for almost identical results to HFSS. Therefore, it can be confidently claimed that the proposed technique can be used as a reliable alternative to existing complex, costly, processor-intensive CAD tools.


2021 ◽  
Vol 11 (9) ◽  
pp. 4100
Author(s):  
Rasa Supankanok ◽  
Sukanpirom Sriwong ◽  
Phisan Ponpo ◽  
Wei Wu ◽  
Walairat Chandra-ambhorn ◽  
...  

Evacuated-tube solar collector (ETSC) is developed to achieve high heating medium temperature. Heat transfer fluid contained inside a copper heat pipe directly affects the heating medium temperature. A 10 mol% of ethylene-glycol in water is the heat transfer fluid in this system. The purpose of this study is to modify inner structure of the evacuated tube for promoting heat transfer through aluminum fin to the copper heat pipe by inserting stainless-steel scrubbers in the evacuated tube to increase heat conduction surface area. The experiment is set up to measure the temperature of heat transfer fluid at a heat pipe tip which is a heat exchange area between heat transfer fluid and heating medium. The vapor/ liquid equilibrium (VLE) theory is applied to investigate phase change behavior of the heat transfer fluid. Mathematical model validated with 6 experimental results is set up to investigate the performance of ETSC system and evaluate the feasibility of applying the modified ETSC in small-scale industries. The results indicate that the average temperature of heat transfer fluid in a modified tube increased to 160.32 °C which is higher than a standard tube by approximately 22 °C leading to the increase in its efficiency by 34.96%.


Sign in / Sign up

Export Citation Format

Share Document