多源并发下Mur二阶吸收边界和非分裂递归卷积完全匹配层对比研究

多个激励源无延时发射(多源并发)相同中心频率脉冲会形成平面波束信号,增强数据记录质量。本文通过数值模拟对比分析在多源并发情况下,非分裂递归卷积完全匹配层作为吸收边界条件和Mur二阶吸收边界条件对电磁波的吸收效果。其研究结果表明,传统的Mur二阶吸收边界条件对多源并发、多角度掠射情况下电磁波的吸收效果不佳,在大偏移距下会造成波形畸变和形成虚假反射。而在多源并发情况下采用非分裂递归卷积完全匹配层作为吸收边界条件,将坐标伸缩因子引进时域有限差分算法中。通过傅里叶逆变换将频率域坐标伸缩变换PML方程转换到时域,对电场和磁场值在离散状态下进行递归卷积运算求解。从而避免了直接对卷积进行数值求解的复杂计算,在保证计算准确性的同时,节约了内存空间,提高了计算效率。在不分裂波场情况下,改善了网格截断位置对电磁波的吸收效果。     

A New Method for Matched Field Localization Based on Two-Hydrophone

The conventional matched field processing (MFP) uses large vertical arrays to locate an underwater acoustic target. However, the use of large vertical arrays increases equipment and computational cost, and causes some problems such as element failures, and array tilting to degrade the localization performance. In this paper, the matched field localization method using two-hydrophone is proposed for underwater acoustic pulse signals with an unknown emitted signal waveform. Using the received signal of hydrophones and the ocean channel pulse response which can be calculated from an acoustic propagation model, the spectral matrix of the emitted signal for different source locations can be estimated by employing the method of frequency domain least squares. The resulting spectral matrix of the emitted signal for every grid region is then multiplied by the ocean channel frequency response matrix to generate the spectral matrix of replica signal. Finally, the matched field localization using two-hydrophone for underwater acoustic pulse signals of an unknown emitted signal waveform can be estimated by comparing the difference between the spectral matrixes of the received signal and the replica signal. The simulated results from a shallow water environment for broadband signals demonstrate the significant localization performance of the proposed method. In addition, the localization accuracy in five different cases are analyzed by the simulation trial, and the results show that the proposed method has a sharp peak and low sidelobes, overcoming the problem of high sidelobes in the conventional MFP due to lack of the number of elements.     

商丘市灾害性天气配套预报方法

在上级指导产品和以数值预报格点资料制作预报系统的基础上,总结非格点气象资料的使用经验,形成本地灾害性天气的补充配套预报方法,从而提高预报准确率。     

二阶声波方程频域PML边界条件及频域变网格步长并行计算

研究了二阶声波方程频域PML边界条件和频域变网格并行计算技术。PML边界是一种较为理想的吸收边界方法,多用在求解时域应力速度方程中,但对于频域声波正演,二阶位移方程更常用。从一阶声波方程PML吸收边界条件导出频域二阶位移方程PML边界条件,模拟算例得到的频率切片、时间切片和地震记录对比都说明该边界条件吸收效果很好。频域单炮正演不同频率间是独立的,据此低频部分采用大网格计算,高频采用小网格,实现变网格步长计算技术,这是较时间域正演的一个优势,在保证模拟质量的同时,减少计算量和内存消耗。     

基于小波分析的磁测数据处理流程及解释方法

对某地区高精度磁测数据,为了有效分析地下场源的异常特征,作者采用小波分析方法对磁异常进行多尺度分析,获取不同深度层次的场源信息,结合平面、剖面资料的反演处理解释,利用2.5维模型反演建立了地下场源的地球物理模型.通过分析,认为小波多尺度分析可以更好的揭示深部场源的赋存信息,计算精度可靠,较传统的延拓、滤波等方法有更高的分辨率.     

二维高分辨自动地电阻率探测地下洞穴技术

高分辨地电阻率系统主要用于地下洞体探测。采用单极－偶极连续二维剖面观测，地采集的资料应用目标异常匹配滤波法或交汇法再现地下异常体形象。     

Comparison of perfectly matched layer and multi-transmitting formula artificial boundary condition based on hybrid finite element formulation

The theory of perfectly matched layer (PML) artificial boundary condition (ABC), which is characterized by absorption any wave motions with arbitrary frequency and arbitrarily incident angle, is introduced. The construc- tion process of PML boundary based on elastodynamic partial differential equation (PDE) system is developed. Combining with velocity-stress hybrid finite element formulation, the applicability of PML boundary is investi- gated and the numerical reflection of PML boundary is estimated. The reflectivity of PML and multi-transmitting formula (MTF) boundary is then compared based on body wave and surface wave simulations. The results show that although PML boundary yields some reflection, its absorption performance is superior to MTF boundary in the numerical simulations of near-fault wave propagation, especially in corner and large angle grazing incidence situations. The PML boundary does not arise any unstable phenomenon and the stability of PML boundary is better than MTF boundary in hybrid finite element method. For a specified problem and analysis tolerance, the computa- tional efficiency of PML boundary is only a little lower than MTF boundary.     

地震声波数值模拟中人工边界条件的差别与组合

地震数值模拟中选取的有限计算区域产生的边界反射会干扰正常波场模拟结果,因此引入人工边界条件来降低边界反射的影响。本文针对PML（perfectly matched layer）边界条件的解耦与非解耦差分形式,应用不同空间差分阶数进行地震波场数值模拟。空间差分阶数提高后,非解耦PML差分形式在计算效率和实现方式上均更具优势。针对CE（Clayton Engquist）边界条件受入射波入射角度限制、边界处精度低吸收效果不好等问题,将2阶CE边界条件和PML边界条件组合成一种新的组合边界条件,在保证吸收效果的同时减少衰减带厚度,从而达到提高计算效率的目的。数值模拟结果验证了算法的有效性。     

优化15点频率-空间域有限差分正演模拟

频率域正演是频率域波形反演的基础,有效快速的正演差分格式可以保证反演结果的精度和效率.本文以用较小的系数矩阵带宽来高效地压制频域正演频散为目标,综合利用加权平均算子、平均加速度项和优化系数三种方法,提出了优化15点差分格式;并且采用压缩存储方式来存放大型系数矩阵,极大地缩小了内存使用量;进而结合最佳匹配层(PML)边界条件,明显地压制了边界反射;最后,通过与前人方法的对比验证,证实了本方法可以在不明显增加计算量的情况下,较好地压制频散.     

Transient analysis of full-space unbounded domains by reciprocal absorbing boundaries

In this paper, we explore a combination of the reciprocal absorbing boundary condition with perfectly matched discrete layers for the time-domain analysis of wave motion in the full plane. Coupling of these absorbing boundaries leads to solutions for out-of-plane as well as in-plane excitations of the full plane. Together with numerical transient simulations, harmonic analysis is also carried out for the purpose of verifying the results by comparing with available exact frequency-domain solutions. The proposed combination is investigated through several numerical examples, which demonstrate its accuracy and efficiency.