薛东川
- 作品数:7 被引量:93H指数:4
- 供职机构:中国石油大学(北京)更多>>
- 发文基金:国家重点基础研究发展计划更多>>
- 相关领域:天文地球更多>>
- 起伏地表复杂介质波动方程有限元数值模拟方法被引量:44
- 2007年
- 波动方程数值模拟是深入研究地震波传播规律的有效方法.有限差分法因其方法简单、精度高而得到了广泛的应用.但其缺点是不能准确模拟具有复杂几何形态的物性界面.因而当遇到起伏地表或复杂构造时,求解精度低.为了准确模拟起伏地形、复杂构造和复杂介质条件下的地震波场,本文采用有限元法模拟二维声波方程.用三角形单元模拟地形和速度界面;把单元内的场和波速均看作单元上的线性函数,以适应复杂介质压制边角散射;采用吸收边界条件去除来自截断边界上的反射;采用集中质量矩阵和集中阻尼矩阵使得显式时间递推无需对矩阵求逆,提高了计算效率.对模型的计算表明该方法正确有效.
- 薛东川王尚旭焦淑静
- 关键词:起伏地表复杂介质有限元数值模拟吸收边界条件
- 利用有限元解波动方程求取初至时间被引量:2
- 2007年
- 叠前逆时偏移需要计算成像点的初至时间。射线追踪法因为计算效率高,常被用作旅行时的计算,但难以处理速度剧变和非规则网格情形。常用的有限元法在求解波动方程时,也会像有限差分法一样存在数据频散问题。本文采用有限元解波动方程获取地震波初至时间,与一般有限元法的主要区别在于利用集中质量矩阵压制时间频散使波前形态保真,并给出了根据震源子波第一极值出现时间求初至时间的方法。数值实验结果表明该方法能适用于非规则网格和速度剧变界面情形。
- 薛东川王尚旭
- 关键词:有限元法初至时间数据处理逆时偏移
- 利用组合质量矩阵压制数值频散被引量:8
- 2008年
- 波动方程数值模拟一直受到频散效应的影响。利用有限元模拟声波方程时,协调质量矩阵和集中质量矩阵对频散误差起着相反的作用:协调质量矩阵使频散误差领先于真实信号,而集中质量矩阵使频散误差落后于真实信号。因此利用两种质量矩阵的线性组合能有效压制数值频散。速度对组合系数的影响很小,在2500~4000m/s范围内,最优组合系数均约为0.5。这种压制数值频散的方法,具有不增加额外计算量且不导致算法复杂的特点,能够适应变速模型。
- 薛东川王尚旭
- 关键词:数值频散有限元
- 起伏地表复杂介质波动方程有限元叠前逆时偏移方法研究
- 随着地震勘探的不断深入和发展,油气勘探的重点正逐渐转向复杂地表条件和复杂地质条件的区域。而山地等复杂地表条件下的地震勘探是当今地球物理界所面临的世界性难题,其中最突出的两个问题就是信噪比低,静校正困难。为了在复杂地区获得...
- 薛东川
- 关键词:起伏地表叠前逆时偏移有限元地震成像地震勘探
- 波动方程有限元叠前逆时偏移被引量:39
- 2008年
- 为研究起伏地表、复杂构造以及复杂速度分布条件下地震数据的精确成像方法,本文采用有限元法完成观测波场的声波方程逆时外推,并提出了以节点波场最大振幅出现时间作为成像时间的最大振幅成像条件。该成像条件不同于激发时间成像条件,它基于波动理论而不是射线理论,可使能量更好地聚焦。有限元法采用六节点三角形单元,对单元内的场和速度均作二次插值以提高模拟精度;用对角化的集中质量矩阵代替协调质量矩阵,避免了矩阵求逆,提高了计算效率。通过对带地形的凹陷模型、Mormousi模型及带地形的Mormousi模型等三个理论模型的计算,结果表明该方法能够有效处理复杂地质条件下的地震成像问题,对于复杂地质体中的断层、断块、背斜构造、高速体和目标体成像都比较清晰,位置准确。尤其对于基于起伏地表直接进行偏移的模型数据也得到了较好的成像结果。
- 薛东川王尚旭
- 关键词:有限元叠前逆时偏移
- 电磁测深2.5维任意复杂介质数值模拟方法研究
- 数值模拟是勘探地球物理的一项重要工作,它既提供各种规律性的认识,又是反演解释的基础.2.5维电磁数值模拟研究三维源在二维构造中电磁场的分布.做法是将空间域的2.5维问题通过沿构造走向的傅立叶变换变换为波数域的二维问题求解...
- 薛东川
- 关键词:勘探地球物理电磁场
- 频率域2.5维电磁测深有限元模拟中的吸收边界条件被引量:7
- 2008年
- 针对频率域2.5维电磁测深问题,借鉴地震波模拟中吸收边界的处理方法,将波数域电磁场方程分解成两个传播方向相反的单程波方程,以沿边界外法向衰减的单程波方程作为该边界上的吸收边界条件。给出了全吸收边界条件的构造方法,并导出了15°吸收边界的具体形式。数值计算结果表明,该吸收边界条件使边界反射得到了有效压制,在相同的计算量下计算精度显著提高。
- 薛东川戴世坤
- 关键词:吸收边界条件电磁测深有限元方法