弹性波正演模拟中改进的非分裂式PML实现方法(英文)

在弹性波有限差分正演模拟中,吸收边界条件常用来吸收截断边界处引入的不期望边界反射,其中完全匹配层(PML)吸收边界条件被认为是目前最理想的吸收边界条件。但是PML吸收边界条件的传统实现却存在着很大不足:全局分裂式PML吸收边界条件实现简单但是需要占用太多内存;局部分裂式PML吸收边界条件需要考虑多个边界和角点区域,编程实现非常复杂;非分裂式PML吸收边界条件由于涉及卷积运算,计算量很大。本文基于非分裂式PML吸收边界条件,结合复频移伸展函数,提出了一种新的数值实现方法,其计算方程简单、占用内存小、编程实现容易,是对PML介质理论数值实现的改进和完善。     

弹性波正演模拟中PML吸收边界条件的改进

在弹性波有限差分正演模拟中, 完全匹配层(PML) 吸收边界条件是使用广泛、吸收效果最好的吸收边界条件.在目前的两种PML实现方法中, 分裂形式的完全匹配层(SPML) 方法计算存储量大、编程实现复杂; 非分裂形式的完全匹配层(NPML) 方法计算效率低、计算过程复杂.针对传统PML吸收边界条件在实现过程中存在的问题, 推导出了一种简洁有效的非卷积实现的NPML吸收边界条件, 既不需要对场分量进行分裂, 也不需要做复杂的卷积运算.分析结果表明, 本文实现的NPML吸收边界条件不仅具有良好的吸收衰减性能, 而且计算方程简单, 编程实现容易, 占有内存更小.      

地质雷达正演中的频散压制和吸收边界改进方法

从麦克斯韦方程组出发,建立了地质雷达的时域有限差分法（ＦＤＴＤ）数学模型,导出了理想 散关系和超级吸收边界条件。理想频散关系 考虑了ＦＤＴＤ法的收敛性和稳定性,也考虑了高频电磁波在Ｙｅｅ氏网格中的传播特点;超吸收边界条件则用磁场分量来提高电场分量精度。数值试验表明,理想频散关系能真实地反映雷达波在地下介质中的传播规律,超吸收边界条件能有效减小截断边界的伪反射,提高正演精度。将之应用于实际计算,取得     

弱各向异性介质弹性波的准各向同性近似正演模拟

准各向同性(QI)近似可用于弱各向异性介质的正演模拟.本文通过运用QI方法的零阶和一阶近似,计算了VTI介质模型的地震记录.得出的地震记录与标准各向同性射线理论(IRT)和基于伪谱法的三维地震正演模拟得出的地震记录作了比较,可以认为是精确的合成地震记录.     

探地雷达时域多分辨法(MRTD)三维正演模拟

应用小波伽略金方法,对Maxwell方程进行离散化,导出了DB2-MRTD算法的探地雷达3D差分公式、数值稳定性条件.在此基础上,开发了探地雷达MRTD(multi-resolution time domain)法正演模拟程序,该程序极大地提高了运算速度,改善了三维探地雷达正演方法,并利用该自制程序,对三角形金属体模型进行了正演模拟,得到了其相应的正演合成三维剖视图及切片图,通过对这些模拟结果进行分析,可以加深对三维雷达反射特征的认识,提高探地雷达探测的可靠性、准确度,同时也说明时域多分辨率法在探地雷达三维正演模拟中的有效性.     

地面瞬变电磁勘探中源外地形的影响

起伏地形为中国矿产资源的发现和勘查造成了极大困难,为了避免复杂的地形效应,在地面瞬变电磁勘探中尽量将发射源与接收点布置在地形平缓处,但源外起伏地形的影响仍然存在。目前,对地形影响的分析大多从测道曲线图或者视电阻率断面图进行定性判断,较少从整个观测时间段内量化评估地形影响的程度。基于有限元法实现了带地形条件下矩形回线源的瞬变电磁三维正演,在源外构造了山脊、山谷地形,将地形与发射源的距离、大地电导率、异常体电导率作为变化参数,给出了整个观测时间段内在不同接收点处地形造成的相对误差分布,从而分析了源外地形对电磁响应的影响规律。结果表明：源外地形对电磁响应的影响主要集中在某个时间段内,随着与发射源距离增大,地形作用减弱,受影响区域也推移到较晚期时间段上;随着大地电导率增大,地形主要影响区域会推移到较晚期时间上,但受影响程度基本一致;相同几何尺寸的山谷地形比山脊地形对电磁响应的影响更大;地形、异常体对电磁响应的影响可以分离,如果异常体产生的二次场比较强,地形效应的影响相对较小,可以忽略地形影响。研究结果为地面瞬变电磁勘探中地形效应的识别与后期处理解释提供了理论依据。     

二维有限元大地电磁正演模拟在Matlab上的实现

Matlab高级语言开发平台对矩阵的计算简便快速,并且绘图输出文件方便。以这些特点为契机,在Matlab平台上,讨论利用有限单元法实现二维大地电磁正演的问题,推导出相应的公式,计算地表或海底视电阻率及相位,并与解析解以及主流方法结果对比,证实了本文方法的正确性。     

复杂地形结构下瞬变电磁三维正演模拟-以大西洋洋中脊TAG热液区为例

如何确定海底多金属硫化物的空间分布并对其进行资源量评价至关重要.瞬变电磁法是陆地金属矿床勘探的重要手段,但海底的崎岖地形、硫化物的复杂结构以及不稳定的近底观测条件为瞬变电磁在海洋中的应用提出了挑战.为验证瞬变电磁方法对热液区硫化物的勘探应用潜力,本文采用了有限元方法,结合大西洋洋中脊TAG热液区实测地形数据和硫化物深部...     

Implications of the Precambrian Non-stromatolitic Carbonate Succession Making up the Third Member of Mesoproterozoic Gaoyuzhuang Formation in Yanshan Area of North China

A particular non-stromatolitic carbonate succession making up the third member of the Mesoproterozoic Gaoyuzhuang (高于庄) Formation might demonstrate that a stromatolite decline of the Mesoproterozoic occurring at ca. 1 450 Ma besides other three events of the Proterozoic,respectively,occurred at ca. 2 000 Ma,ca. 1 000 Ma,and ca. 675 Ma. The forming duration of this non-stromatolitic carbonate succession can be generally correlative to that of a similar depositional succession in North America,i.e. a non-stromatolitic carbonate succession made up by the Helena Formation of the Belt Supergroup,which suggests that the stromatolite decline occurring at ca. 1 450 Ma may be a global event. This information endows the non-stromatolitic carbonate succession making up the third member of the Gaoyuzhuang Formation in the Yanshan (燕山) area with important significance for the further understanding of Precambrian sedimentology. The Mesoproterozoic Gaoyuzhuang Formation in Yanshan area is a set of more than 1 000 m thick carbonate strata that can be divided into four members (or subformations). The first member (or the Guandi (官地) subformation) is marked by a set of stromatolitic dolomites overlying a set of transgressive sandstones; the second member (or the Sangshu'an (桑树鞍) subformation) is a set of manganese dolomites with a few stromatolites; the third member (or the Zhangjiayu (张家峪) subformation) is chiefly made up of leiolite and laminite limestones and is characterized by the development of molar-tooth structures in leiolite limestone; the fourth member (or the Huanxiusi (环秀寺) subformation) is composed of a set of dolomites of stromatolitic reefs or lithoherms. Sequence-stratigraphic divisions at two sections,i.e. the Jixian (蓟县) Section in Tianjin (天津) and the Qiangou (千沟) Section of Yanqing (延庆) County in Beijing (北京),demonstrate that a particularly non-stromatolitic succession making up the third member of the Mesoproterozoic Gaoyuzhuang Formation is developed in the Yanshan area of North China,in which lots of grotesque matground structures (wrinkle structures and palimpsest ripples) are developed in beds of leiolite limestone at the Qiangou Section and lots of molar-tooth structures are developed in beds of leiolite limestone at the Jixian Section. The time scale of the Gaoyuzhuang Formation is deduced as 200 Ma (from 1 600 Ma to 1 400 Ma). The duration of an obvious hiatus between the Gaoyuzhuang Formation and the underlying Dahongyu (大红峪) Formation is deduced as 50 Ma to 100 Ma,thus the forming duration of the GaoyuzhuangFormation is thought as 100 Ma (1 500 Ma to 1 400 Ma). Furthermore,the age of the subface of the third member of the Gaoyuzhuang Formation that is just in the mid position of the Gaoyuzhuang Formation can be deduced as about 1 450 Ma,which is the basis to infer a stromatolite decline of the Mesoproterozoic occurring at ca. 1 450 Ma. Importantly,several features of both the molar-tooth structure and the stromatolite,such as the particular forming environment,the important facies-indicative meaning,and the episodic distribution in the earth history,might express the evolutionary periodicity of the surface environment of the earth and can provide meaningful clues for the understanding of the Precambrian world,although their origin and forming mechanism is highly contentious. Therefore,like other three stromatolitic declines,respectively,occurring at ca. 675 Ma,ca. 1 000 Ma,and ca. 2 000 Ma,the identification of the stromatolite decline occurring at ca. 1 450 Ma during the Golden Age of stromatolites (2 800 Ma to 1 000 Ma) has important meaning for the further understanding of the evolving carbonate world of the Precambrian.     

高频大地电磁法有限元数值模拟

由于中深度工程地球物理勘探的迫切需要,以美国EH-4电导率成像系统为代表的高频率大地电磁法,在我国地球物理勘探行业的应用越来越广泛。其勘探深度在地下1000m以内,频率范围在10Hz～100kHz。这里研究了高频率大地电磁法有限元高精度,快速数值模拟。首先利用广义变分原理推导出电磁场的有限元方程,编制了双线性和双二次插值有限元程序,并进行了层状模型的验证性计算,然后对典型地堑模型的数值模拟结果进行了分析,得出了一些有意义的结论。