Viscoacoustic prestack reverse time migration based onthe optimal time-space domain high-order finite-difference method

Prestack reverse time migration （RTM） is an accurate imaging method ofsubsurface media. The viscoacoustic prestack RTM is of practical significance because itconsiders the viscosity of the subsurface media. One of the steps of RTM is solving thewave equation and extrapolating the wave field forward and backward; therefore, solvingaccurately and efficiently the wave equation affects the imaging results and the efficiencyof RTM. In this study, we use the optimal time-space domain dispersion high-order finite-difference （FD） method to solve the viscoacoustic wave equation. Dispersion analysis andnumerical simulations show that the optimal time-space domain FD method is more accurateand suppresses the numerical dispersion. We use hybrid absorbing boundary conditions tohandle the boundary reflection. We also use source-normalized cross-correlation imagingconditions for migration and apply Laplace filtering to remove the low-frequency noise.Numerical modeling suggests that the viscoacoustic wave equation RTM has higher imagingresolution than the acoustic wave equation RTM when the viscosity of the subsurface isconsidered. In addition, for the wave field extrapolation, we use the adaptive variable-lengthFD operator to calculate the spatial derivatives and improve the computational efficiencywithout compromising the accuracy of the numerical solution.

Benchmarking of two three-dimensional numerical models in time/space domain to predict railway-induced ground vibrations

In the last 30 years,the scientific community has developed and proposed different models and numerical approaches for the study of vibrations induced by railway traffic.Most of them are formulated in the frequency/wave number domain and with a 2.5D approach.Three-dimensional numerical models formulated in the time/space domain are less frequently used,mainly due to their high computational cost.Notwithstanding,these models present very attractive characteristics,such as the possibility of considering nonlinear behaviors or the modelling of excess pore pressure and non-homogeneous and non-periodic geometries in the longitudinal direction of the track.In this study,two 3D numerical approaches formulated in the time/space domain are compared and experimentally validated.The first one consists of a finite element approach and the second one of a finite difference approach.The experimental validation in an actual case situated in Carregado(Portugal)shows an acceptable fitting between the numerical results and the actual measurements for both models.However,there are some differences among them.This study therefore includes some recommendations for their use in practical soil dynamics and geotechnical engineering.

3D Gravity Inversion with Correlation Image in Space Domain and Application to the Northern Sinai Peninsula

We present a 3D inversion method to recover density distribution from gravity data in space domain.Our method firstly employs 3D correlation image of the vertical gradient of gravity data as a starting model to generate a higher resolution image for inversion.The 3D density distribution is then obtained by inverting the correlation image of gravity data to fit the observed data based on classical inversion method of the steepest descent method.We also perform the effective equivalent storage and subdomain techniques in the starting model calculation,the forward modeling and the inversion procedures,which allow fast computation in space domain with reducing memory consumption but maintaining accuracy.The efficiency and stability of our method is demonstrated on two sets of synthetic data and one set of the Northern Sinai Peninsula gravity data.The inverted 3D density distributions show that high density bodies beneath Risan Aniza and low density bodies exist to the southeast of Risan Aniza at depths between 1~10 and 20 km,which may be originated from hot anomalies in the lower crust.The results show that our inversion method is useful for 3D quantitative interpretation.

Modeling of Surface Waves in a Fluid Saturated Poro-Elastic Medium under Initial Stress Using Time-Space Domain Higher Order Finite Difference Method

In this present context, mathematical modeling of the propagation of surface waves in a fluid saturated poro-elastic medium under the influence of initial stress has been considered using time dependent higher order finite difference method (FDM). We have proved that the accuracy of this finite-difference scheme is 2M when we use 2nd order time domain finite-difference and 2M-th order space domain finite-difference. It also has been shown that the dispersion curves of Love waves are less dispersed for higher order FDM than of lower order FDM. The effect of initial stress, porosity and anisotropy of the layer in the propagation of Love waves has been studied here. The numerical results have been shown graphically. As a particular case, the phase velocity in a non porous elastic solid layer derived in this paper is in perfect agreement with that of Liu et al. (2009).

Multilinear Calder′on Zygmund operators on variableexponent Morrey spaces over domains

The boundedness of multilinear singular integrals of Calder′on-Zygmund type onproduct of variable exponent Lebesgue spaces over both bounded and unbounded domains areobtained. Further more, the boundedness for this type multilinear operators on product ofvariable exponent Morrey spaces over domains is shown in the paper.

Solution Domain Boundary Analysis Method and Its Application in Parameter Spaces of Nonlinear Gear System

Mastering the influence laws of parameters on the solution structure of nonlinear systems is the basis of carrying out vibration isolation and control.Many researches on solution structure and bifurcation phenomenon in parameter spaces are carried out broadly in many fields,and the research on nonlinear gear systems has attracted the attention of many scholars.But there is little study on the solution domain boundary of nonlinear gear systems.For a periodic non-autonomous nonlinear dynamic system with several control parameters,a solution domain boundary analysis method of nonlinear systems in parameter spaces is proposed,which combines the cell mapping method based on Poincaré point mapping in phase spaces with the domain decomposition technique of parameter spaces.The cell mapping is known as a global analysis method to analyze the global behavior of a nonlinear dynamic system with finite dimensions,and the basic idea of domain decomposition techniques is to divide and rule.The method is applied to analyze the solution domain boundaries in parameter spaces of a nonlinear gear system.The distribution of different period domains,chaos domain and the domain boundaries between different period domains and chaotic domain are obtained in control parameter spaces constituted by meshing damping ratio with excitation frequency,fluctuation coefficient of meshing stiffness and average exciting force respectively by calculation.The calculation results show that as the meshing damping increases,the responses of the system change towards a single motion,while the variations of the excitation frequency,meshing stiffness and exciting force make the solution domain presenting diversity.The proposed research contribution provides evidence for vibration control and parameter design of the gear system,and confirms the validity of the solution domain boundary analysis method.

基于CiteSpace的用户画像知识图谱分析

为了明确用户画像研究热点和趋势,文章运用CiteSpace软件定量化分析知网1496篇用户画像文献。研究结果表明,国内对用户画像研究呈现增长态势,总体经历了研究萌芽、快速发展和质量提升三个阶段,虽然研究机构和作者群体稳定,但尚未形成核心作者群和跨专业领域的合作。研究趋势方面,大数据、数据挖掘、图书馆、机器学习、人工智能、精准营销、知识服务、协同过滤等是主要研究热点;前沿研究聚焦于信息服务、融媒体、数据采集、媒体融合和聚类。

Three-dimensional numerical modeling of gravity anomalies based on Poisson equation in spacewavenumber mixed domain

In gravity-anomaly-based prospecting, the computational and memory requirements for practical numerical modeling are potentially enormous. Achieving an efficient and precise inversion for gravity anomaly imaging over large-scale and complex terrain requires additional methods. To this end, we have proposed a new topography-capable By performing a two-dimensional Fourier transform in the horizontal directions, threedimensional partial differential equations in the spatial domain were transformed into a group of independent, one-dimensional differential equations engaged with different wave numbers. These independent differential equations are highly parallel across different wave numbers. differential equations with different wave numbers, and the efficiency of solving fixedbandwidth linear equations was further improved by a chasing method. In a synthetic test, a prism model was used to verify the accuracy and reliability of the proposed algorithm by comparing the numerical solution with the analytical solution. We studied the computational precision and efficiency with and without topography using different Fourier transform methods. The results showed that the Guass-FFT method has higher numerical precision, while the standard FFT method is superior, in terms of computation time, for inversion and quantitative interpretation under complicated terrain.

THE EQUIVALENT CHARACTERIZATION OF F(p,q,s)SPACE ON BOUNDED SYMMETRIC DOMAINS OF Cn

LetΩ be a bounded symmetric domain in Cn. The purpose of this article is to define and characterize the general function space F（p, q, s） on Ω. Characterizing functions in the F（p, q, s） space is a work of considerable interest nowadays. In this article, the authors give several equivalent descriptions of the functions in the F（p, q, s） space on Ω in terms of fractional differential operators. At the same time, the authors give the relationship between F（p, q, s） space and Bloch type space on Ω too.

Compactness of Composition Operators from the p-Bloch Space to the q-Bloch Space on the Classical Bounded Symmetric Domains

In this paper, we introduce the weighted Bloch spaces on the first type of classical bounded symmetric domains , and prove the equivalence of the norms and . Furthermore, we study the compactness of composition operator from to , and obtain a sufficient and necessary condition for to be compact.

Three-dimensional numerical simulation of geothermal fi eld in space-wavenumber domain

Large-scale,fine,and efficient numerical simulation of a geothermal field plays an important role in geothermal energy development.Confronted with the problem of large computation and high storage requirements for complex underground models in a three-dimensional(3-D)numerical simulation of a geothermal fi eld,a mixed space-wavenumber domain 3-D numerical simulation algorithm is proposed in this paper.According to the superposition principle of temperature field,the geothermal field is decomposed into background and abnormal temperature fi elds for calculation.The uniform layered model is used to solve the background field.When the abnormal field is solved,the horizontal two-dimensional(2-D)Fourier transform is used to transform the 3-D diff erential equation satisfi ed by an abnormal field into a series of one-dimensional ordinary differential equations with diff erent wavenumbers,which greatly reduces the calculation and storage.The unit division of an ordinary diff erential equation is fl exible,and the calculation amount is small.The algorithm fully takes advantage of the effi ciency of the Fourier transform and the quickness of the catch-up method to solve linear equations with a fixed bandwidth,which effectively improves the computational efficiency.Compared with the COMSOL Multiphysics professional simulation finite element software,the time consumption and memory requirements of the algorithm proposed in this paper are reduced by multiple orders of magnitude in terms of ensuring accuracy and the same mesh division.The more the number of calculated nodes is,the more obvious is the advantage.We design models to study the thermal conductivity,heat fl ux boundary,regional tectonic morphology,and topographic relief of the geothermal fi eld distribution.A 3-D geophysical model is developed based on topographic elevation data,geothermal geology,and geophysical exploration data in the Qiabuqia area of Gonghe Basin,Qinghai Province,China.Numerical simulation of the geothermal fi eld in this area is realized,which shows that the algorithm is suitable for precise and effi cient simulation of an arbitrary complex terrain and geological conditions.

Direct data domain approach to space-time adaptive processing

In non-homogeneous environment, traditional space-time adaptive processing doesn＇t effectively suppress interference and detect target, because the secondary data don＇ t exactly reflect the statistical characteristic of the range cell under test. A ravel methodology utilizing the direct data domain approach to space-time adaptive processing （ STAP ） in airbome radar non-homogeneous environments is presented. The deterministic least squares adaptive signal processing technique operates on a ＂snapshot-by-snapshot＂ basis to dethrone the adaptive adaptive weights for nulling interferences and estimating signal of interest （SOI）. Furthermore, this approach eliminates the requirement for estimating the covariance through the data of neighboring range cell, which eliminates calculating the inverse of covariance, and can be implemented to operate in real-time. Simulation results illustrate the efficiency of interference suppression in non-homogeneous environment.

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.

Weighted Composition Operators from the Bloch Space to Weighted Banach Spaces on Bounded Homogeneous Domains

We study the bounded and the compact weighted composition operators from the Bloch space into the weighted Banach spaces of holomorphic functions on bounded homogeneous domains, with particular attention to the unit polydisk. For bounded homogeneous domains, we characterize the bounded weighted composition operators and determine the operator norm. In addition, we provide sufficient conditions for compactness. For the unit polydisk, we completely characterize the compact weighted composition operators, as well as provide ＂computable＂ estimates on the operator norm.

LARGE TIME ASYMPTOTIC BEHAVIOR OF THE COMPRESSIBLE NAVIER-STOKES EQUATIONS IN PARTIAL SPACE-PERIODIC DOMAINS

In this article, we study the large time behavior of the 3-D isentropic compressible Navier-Stokes equation in the partial space-periodic domains, and simultaneously show that the related profile systems can be described by like Navier-Stokes equations with suitable ＂pressure＂ functions in lower dimensions. Our proofs are based on the energy methods together with some delicate analysis on the corresponding linearized problems.

EQUIVALENCES BETWEEN BLOCH TYPE SPACES AND BMO SPACES ON THE HYPERBOLIC DOmAINS IN C

For an analytic function f on the hyperbolic domain Ω in C,the following conclusions are obtained: (i)f∈B(Ω)=BMOA(Ω,m)if and only if Ref∈B(?)(Ω)=BMOH(Ω,m).(ii)QB_h(Ω)=B_h(Ω) (BMOH,(Ω,m)=BMOH(Ω,m)if and only if C(Ω)=inf{Z_o(z)·δ_o(z)·z≡Ω}>0,Also some applica- lions to automorphic function are considered.

Toeplitz Operators on the Bergman Space of Multi-connected Domains

In this paper, we investigate the Toeplitz operators with positive measure symbols on the Bergman spaces of bounded multi-connected domains and show that a Toeplitz operator is bounded or compact if and only if the symbol measure is a Carleson or vanishing Carleson measure respectively.

POTENTIAL SPACE ESTIMATES IN LOCAL HARDY SPACES FOR GREEN POTENTIALS IN CONVEX DOMAINS

Let Ω be a bounded convex domain in Rn(n≥3) and G(x,y) be the Green function of the Laplace operator -△on Ω Let hPT(Ω) = {f∈D'(Ω) : (?)F ∈ hP(Rn), s.t. F|Ω = f}, by the atom characterization of Local Hardy spaces in a bounded Lipschitz domain, the bound of f Ω(?)2(Gf) for every f ∈hPr(Ω) is obtained, where n/(n + 1) <p≤1.

任意起伏地形下重力异常三维正演及并行计算

为了进一步提高空间-波数域三维重力异常正演算法的适用范围和计算效率,本文采用任意傅里叶变换算法实现了空间-波数域三维重力异常正演,且在NVIDIA CUDA平台上进行CPU-GPU并行加速.任意傅里叶变换算法的基本思想是将二维傅里叶变换转化为两个一维傅里叶变换,一维傅里叶变换积分离散为多个单元积分累加和,离散单元中原函数采用二次插值形函数拟合,求出单元积分的解析表达式.相比现有的傅里叶变换算法,新方法具有采样灵活、积分精度高、计算速度快和傅里叶变换的截断效应小等优势.利用空间-波数域算法的高度并行性,采用CPU并行求解常微分方程,GPU并行计算任意傅里叶变换,实现了CPU-GPU并行加速方案,进一步提升了本文算法效率.利用常密度模型,对比数值解和解析解,结果表明本文算法正确;利用变密度模型对比了任意傅里叶变换算法与高斯快速傅里叶变换算法的计算效率与精度,在相近的数值精度下,本文算法波数选取少,效率高;测试CPU-GPU并行效果,结果表明相比CPU串行算法,CPU-GPU并行算法的计算效率大大提升,千万数量级节点数模型正演仅耗时数秒.最后利用实际地形数据进行三维重力异常场数值模拟,证明了新方法的高效性与实用性,对实现大规模复杂条件下重力异常精细化反演成像与综合解释有重要意义.