Multi-scale seismic full waveform inversion in the frequency-domain with a multi-grid method

Although full waveform inversion in the frequency domain can overcome the local minima problem in the time direction, such problem still exists in the space direction because of the media subsurface complexity. Based on the optimal steep descent methods, we present an algorithm which combines the preconditioned bi-conjugated gradient stable method and the multi-grid method to compute the wave propagation and the gradient space. The multiple scale prosperity of the waveform inversion and the multi-grid method can overcome the inverse problems local minima defect and accelerate convergence. The local inhomogeneous three-hole model simulated results and the Marmousi model certify the algorithm effectiveness.

Harmonic balance method with alternating frequency/time domain technique for nonlinear dynamical system with fractional exponential

Comparisons of the common methods for obtaining the periodic responses show that the harmonic balance method with alternating frequency/time （HB-AFT） do- main technique has some advantages in dealing with nonlinear problems of fractional exponential models. By the HB-AFT method, a rigid rotor supported by ball bearings with nonlinearity of Hertz contact and ball passage vibrations is considered. With the aid of the Floquet theory, the movement characteristics of interval stability are deeply studied. Besides, a simple strategy to determine the monodromy matrix is proposed for the stability analysis.

Application of pseudo-random frequency domain electromagnetic method in mining areas with strong interferences

Due to the strong electromagnetic interferences and human interference,traditional electromagnetic methods cannot obtain high quality resistivity data of mineral deposits in Chinese mines.The wide field electromagnetic method(WFEM),in which the pseudo-random signal is taken as the transmitter source,can extract high quality resistivity data in areas with sever interference by only measuring the electric field component.We use the WFEM to extract the resistivity information of the Dongguashan mine in southeast China.Compared with the audio magnetotelluric(AMT)method,and the controlled source audio-frequency magnetotelluric(CSAMT) method,the WFEM can obtain data with higher quality and simpler operations.The inversion results indicate that the WFEM can accurately identify the location of the main ore-body,which can be used for deep mine exploration in areas with strong interference.

Acoustic reflection well logging modeling using the frequency-domain finite-element method with a hybrid PML

In this paper, we propose a hybrid PML （H-PML） combining the normal absorption factor of convolutional PML （C-PML） with tangential absorption factor of Mutiaxial PML （M-PML）. The H-PML boundary conditions can better suppress the numerical instability in some extreme models, and the computational speed of finite-element method and the dynamic range are greatly increased using this HPML. We use the finite-element method with a hybrid PML to model the acoustic reflection of the interface when wireline and well logging while drilling （LWD）, in a formation with a reflector outside the borehole. The simulation results suggests that the PS- and SP- reflected waves arrive at the same time when the inclination between the well and the outer interface is zero, and the difference in arrival times increases with increasing dip angle. When there are fractures outside the well, the reflection signal is clearer in the subsequent reflection waves and may be used to identify the fractured zone. The difference between the dominant wavelength and the model scale shows that LWD reflection logging data are of higher resolution and quality than wireline acoustic reflection logging.

Modified slow-fast analysis method for slow-fast dynamical systems with two scales in frequency domain

A modified slow-fast analysis method is presented for the periodically excited non-autonomous dynamical system with an order gap between the exciting frequency and the natural frequency.By regarding the exciting term as a slow-varying parameter,a generalized autonomous fast subsystem can be defined,the equilibrium branches as well as the bifurcations of which can be employed to account for the mechanism of the bursting oscillations by combining the transformed phase portrait introduced.As an example,a typical periodically excited Hartley model is used to demonstrate the validness of the method,in which the exciting frequency is far less than the natural frequency.The equilibrium branches and their bifurcations of the fast subsystem with the variation of the slow-varying parameter are presented.Bursting oscillations for two typical cases are considered,which reveals that,fold bifurcation may cause the the trajectory to jump between different equilibrium branches,while Hopf bifurcation may cause the trajectory to oscillate around the stable limit cycle.

Comparison of delay compensation methods for real-time hybrid simulation using frequency-domain evaluation index

The delay compensation method plays an essential role in maintaining the stability and achieving accurate real-time hybrid simulation results. The effectiveness of various compensation methods in different test scenarios, however, needs to be quantitatively evaluated. In this study, four compensation methods （i.e., the polynomial extrapolation, the linear acceleration extrapolation, the inverse compensation and the adaptive inverse compensation） are selected and compared experimentally using a frequency evaluation index （FEI） method. The effectiveness of the FEI method is first verified through comparison with the discrete transfer fimction approach for compensation methods assuming constant delay. Incomparable advantage is further demonstrated for the FEI method when applied to adaptive compensation methods, where the discrete transfer function approach is difficult to implement. Both numerical simulation and laboratory tests with predefined displacements are conducted using sinusoidal signals and random signals as inputs. Findings from numerical simulation and experimental results demonstrate that the FEI method is an efficient and effective approach to compare the performance of different compensation methods, especially for those requiring adaptation of compensation parameters.

Goal-oriented error estimation applied to direct solution of steady-state analysis with frequency-domain finite element method

Based on the concept of the constitutive relation error along with the residuals of both the origin and the dual problems, a goal-oriented error estimation method with extended degrees of freedom is developed. It leads to the high quality locM error bounds in the problem of the direct-solution steady-state dynamic analysis with a frequency-domain finite element, which involves the enrichments with plural variable basis functions. The solution of the steady-state dynamic procedure calculates the harmonic response directly in terms of the physical degrees of freedom in the model, which uses the mass, damping, and stiffness matrices of the system. A three-dimensional finite element example is carried out to illustrate the computational procedures.

A Krylov Space-Based Finite Element Time Domain Method for Broadband Frequency Domain Solutions

A Krylov space based time domain method for wave propagation problems is introduced. The proposed method uses the Arnoldi algorithm to obtain broad-band frequency domain solutions. This method is especially advantageous in cases where slow convergence is observed when using traditional time domain methods. The efficiency of the method is examined in several test cases to show its fast convergence in such problems.

Acoustic viscoelastic modeling by frequency-domain boundary element method

Earth medium is not completely elastic, with its viscosity resulting in attenuation and dispersion of seismic waves. Most viscoelastic numerical simulations are based on the finite-difference and finite-element methods. Targeted at viscoelastic numerical modeling for multilayered media, the constant-Q acoustic wave equation is transformed into the corresponding wave integral representation with its Green＇s function accounting for viscoelastic coefficients. An efficient alternative for full-waveform solution to the integral equation is proposed in this article by extending conventional frequency-domain boundary element methods to viscoelastic media. The viscoelastic boundary element method enjoys a distinct characteristic of the explicit use of boundary continuity conditions of displacement and traction, leading to a semi-analytical solution with sufficient accuracy for simulating the viscoelastic effect across irregular interfaces. Numerical experiments to study the viscoelastic absorption of different Q values demonstrate the accuracy and applicability of the method.

Frequency-time Domain Shielding Effectiveness of Materials and Their Strong Electromagnetic Field Test

To describe the shielding ability of materials accurately and comprehensively,the frequency-domain and time-domain shielding effectiveness(SE) of material is investigated.The relevance between them was analyzed based on the minimum phase method,and the time-domain SE can be derived from frequency-domain SE.The SE of an energy selective surface(ESS) made of a novel material is investigated,and the relationship between SE and radiation field intensity are analyzed.The results show that not only material,but also the intensity of radiation electric field shows influence on SE in its frequency; for some materials,the dependence of SE on radiation electric field intensity needs to be considered.Therefore,it is necessary to research on the SE of shielding material in high-intensity electromagnetic environment.

Characteristic signal extracted from a continuous time signal on the aspect of frequency domain

Extracting characteristic signal from a continuous signal can effectively reduce the difficulty of analyzing the running states of a single-variable nonlinear system.Whether the extracted characteristic signal can accurately reflect the running states of the system is very important.In this paper, a method called automatic sampling method(ASM) for extracting characteristic signals is investigated.The complete definition is described, the effectiveness is proved theoretically, and the general formulas of the extracted characteristic signals are derived for the first time.Furthermore, typical Chua's circuit is used to accomplish a lot of experimental research on the aspect of frequency domain.The experimental results show that ASM is feasible and practical, and can automatically generate a characteristic signal with the change of the original signal.

Error Analysis in Frequency Domain for Linear Multipass Algorithms

Error analysis methods in frequency domain are developed in this paper for determining the characteristic root and transfer function errors when the linear multipass algorithms are used to solve linear differential equations. The relation between the local truncation error in time domain and the error in frequency domain is established, which is the basis for developing the error estimation methods. The error estimation methods for the digital simulation model constructed by using the Runge-Kutta algorithms and the linear multistep predictor-corrector algorithms are also given.

THE STUDY ON MEASURING TECHNIQUE OF PARTIAL DISCHARGE IN GAS INSULATED SWITCHGEAR USING ULTRA HIGH FREQUENCY METHOD WITH EXTERNAL SENSORS

Objective In order to find early latent faults and prevent catastrophic failures, diagnosis of insulation condition by measuring technique of partial discharge(PD) in gas insulated switchgear (GIS) is applied in this paper, which is one of the most basic ways for diagnosis of insulation condition. Methods Ultra high frequency(UHF) PD detection method by using internal sensors has been proved efficient, because it may avoid the disturbance of corona, but the sensor installation of this method will be limited by the structure and operation condition of GIS. There are some of electromagnetic (E-M) waves leak from the place of insulation spacer, therefore, the external sensors UHF measuring PD technique is applied, which isn't limited by the operation condition of GIS. Results This paper analyzes propagated electromagnetic (E-M) waves of partial discharge pulse excited by using the finite-difference time-domain (FDTD) method. The signal collected at the outer point is more complex than that of the inner point, and the signals' amplitude of outer is about half of the inner, because it propagates through spacer and insulation slot. Set up UHF PD measuring system. The typical PD in 252kV GIS bus bar was measured using PD detection UHF technique with external sensors. Finally, compare the results of UHF measuring technique using external sensors with the results of FDTD method simulation and the traditional IEC60270 method detection. Conclusion The results of experiment shows that the UHF technique can realize the diagnosis of insulation condition, the results of FDTD method simulation and the result UHF method detection can demonstrate each other, which gives references to further researches and application for UHF PD measuring technique.

Applying a Step Approach Method in Solving the Multi-Frequency Radiation From a Complex Obstacle

In this paper,a step approach method in the time domain is developed to calculate the radiated waves from an arbitrary obstacle pulsating with multiple frequencies.The computing scheme is based on the Boundary Integral Equation and derived in the time domain;thus,the time-harmonic Neumann boundary condition can be imposed.By the present method,the values of the initial conditions are set to zero,and the approach process is carried forward in a loop from the first time step to the last.At each time step,the radiated pressure on each element is updated.After several loops,the correct radiated pressures can be obtained.A sphere pulsating with a monopole frequency in an infinite acoustic domain is calculated first.This result is compared with the analytical solution,and both of them are in good agreement.Then,a complex-shaped radiator is taken as the studied case.The pulsating frequency of this case is multiple,and the waves propagate in half space.It is shown that the present method can treat multiple-frequency pulsation well,even when the radiator is a complex shape,and a robust convergence can be attained quickly.

Efficient solution of large-scale matrix of acoustic wave equations in 3D frequency domain

In 3D frequency domain seismic forward and inversion calculation,the huge amount of calculation and storage is one of the main factors that restrict the processing speed and calculation efficiency.The frequency domain finite-difference forward simulation algorithm based on the acoustic wave equation establishes a large bandwidth complex matrix according to the discretized acoustic wave equation,and then the frequency domain wave field value is obtained by solving the matrix equation.In this study,the predecessor's optimized five-point method is extended to a 3D seven-point finite-difference scheme,and then a perfectly matched layer absorbing boundary condition(PML)is added to establish the corresponding matrix equation.In order to solve the complex matrix,we transform it to the equivalent real number domain to expand the solvable range of the matrix,and establish two objective functions to transform the matrix solving problem into an optimization problem that can be solved using gradient methods,and then use conjugate gradient algorithm to solve the problem.Previous studies have shown that in the conjugate gradient algorithm,the product of the matrix and the vector is the main factor that affects the calculation efficiency.Therefore,this study proposes a method that transform bandwidth matrix and vector product problem into some equivalent vector and vector product algorithm,thereby reducing the amount of calculation and storage.

Hydroelasticity Analysis in Frequency Domain and Time Domain

Hydroelasticity has been introduced in ship seakeeping assessment for more than three decades, and it finally becomes an essential tool in marine industry for design of some types of ship. In the 35 years of evolution, hydroelasticity methods applied in industry of marine and offshore energy grown up from two dimensional to three dimensional and now has analysis models of linear model in frequency domain and nonlinear model in time domain. In this paper, we present the three dimensional hydroelasticity theory model in frequency domain and time domain, show the difference in the approach, and discuss their applications in wave-structure interaction.

Finite difference time domain method forward simulation of complex geoelectricity ground penetrating radar model

The ground penetrating radar(GPR) forward simulation all aims at the singular and regular models, such as sandwich model, round cavity, square cavity, and so on, which are comparably simple. But as to the forward of curl interface underground or “v” figure complex model, it is difficult to realize. So it is important to forward the complex geoelectricity model. This paper takes two Maxwell’s vorticity equations as departure point, makes use of the principles of Yee’s space grid model theory and the basic principle finite difference time domain method, and deduces a GPR forward system of equation of two dimensional spaces. The Mur super absorbed boundary condition is adopted to solve the super strong reflection on the interceptive boundary when there is the forward simulation. And a self-made program is used to process forward simulation to two typical geoelectricity model.

超临界循环流化床锅炉深度调峰时水冷壁流动不稳定特性计算分析

针对深度调峰时超临界循环流化床锅炉水冷壁管道内的汽液两相流动不稳定问题,提出了一种适用于不同工况的频域法数学模型。通过对质量、能量、动量方程进行小扰动线性化后经拉普拉斯变换得到了用于描述管道内汽液流体流动稳定性的传递函数,通过图解奈奎斯特图的方法判断管内工质流动的稳定性。然后,通过该模型计算了某350 MW超临界循环流化床锅炉水冷壁管段的不稳定边界并研究了不同参数对流动不稳定特性的影响。计算结果表明:所选回路在20%和50%锅炉最大连续蒸发量负荷发生流动不稳定时,热负荷分别为76.09 kW·m^(-2)和113.52 kW·m^(-2),水冷壁管道流动稳定且安全;入口过冷度对临界热流密度的影响呈非单值性;增大质量流速使流体的进口和出口密度差减小,有利于流动的稳定;增大入口节流系数,可以抑制入口处流量的脉动,有利于流动的稳定;在不同的工况条件下,改变管段倾斜角度会对流动稳定性产生不同的影响。

非水平海底地形下浅海全波导声场波动成分分析

为分析浅海非水平海底地形下甚低频声波激发海底地震波场波动成分,阐明海底地形变化对水中声源所激发海底地震波场的影响规律,开展了理论建模与数值仿真研究。基于时域有限元法,在柱坐标系下推导并建立了一种浅海水体/海底全波导下的声场计算模型,结合仿真算例讨论了上坡海底、下坡海底、海底山三类典型非水平海底地形对海底地震波场成分的影响规律及机理。研究结果表明,利用时域有限元方法可以准确计算浅海各类非水平海底地形下的地震波场;水中声源所激发海底地震波场包含水中声波、海底界面波、海底纵波、海底横波4个波系;与理想水平海底相比,上坡海底会加速声能量的衰减,且坡度越大,声波远距离传播时向海底泄漏的越多,声能量衰减越大。而下坡海底的波场变化则相反,海底山的存在使得部分声波发生反射及透射到海底的现象,声能量衰减程度增加。