共查询到20条相似文献，搜索用时 15 毫秒
1.
《物探化探计算技术》2017,(1)
利用交错网格有限差分技术实现了频率域海洋可控源电磁三维正演算法,并探讨了海底地形对海洋可控源电磁响应的影响和校正方法及效果。通过与海底二维山峰地形的二维自适应有限元解对比,验证了文中算法的正确性和模拟海底起伏地形的可行性,并分析了海底二维山峰地形对三维油气储层电场响应的影响。关于海底起伏地形对海洋可控源电场响应的影响,采用比较法对其进行了地形校正,结果表明,在深水区比较法可以很好地消除地形影响,恢复海底高阻油气储层的电场响应。 相似文献
2.
时域瞬变场电磁场有限差分法 总被引：1，自引：2，他引：1
模拟电磁场在地下介质中的传播规律是理解地下介质响应的重要手段，直接应用有限差分法在时间域对二维、三维瞬变场进行数值分析是一种有效方法。应用反映电磁场基本规律的麦克斯韦方程组，导出时域场量的齐次扩散方程，对所研究的空间区域作差分离散，源作为初始条件加入，采用合适的边界条件，利用差分方程进行计算，从而得到场量的数值结果，展现瞬变场在地下随时间扩散的全过程，对电磁场的瞬态响应有更加直观的了解。 相似文献
3.
本文研究了利用阻抗信息进行可控源电磁勘探有限内存拟牛顿法三维反演的技术。首先用理论模型来验证有限内存拟牛顿法反演的准确性和可行性。观测参数为复阻抗Zxy分量,采用交错网格有限差分方法计算模型响应,反演采用有限内存拟牛顿法。数值模拟结果表明：有限内存拟牛顿法反演迭代速度较快,每4 min迭代一次,拟合差由146.00下降到1.78,收敛稳定;异常体的位置与理论模型吻合较好,有效地验证了有限内存拟牛顿法可控源三维反演的正确性。为了进一步验证该方法的实用性,将其应用到隐伏钼矿可控源电磁勘探工作中。工区的反演结果显示：在工区北西段深部存在高阻异常,其上为低阻异常。截取过钻孔的3号测线发现,-500~-100 m的位置表现为低阻,东侧存在向上涌起的高阻。推测此低阻为矿化蚀变带,延伸较深。该异常与钻探资料揭示的钼矿脉一致,证明了反演结果的准确性。因此,利用可控源观测得到的阻抗信息进行有限内存拟牛顿法三维反演,可以获得可靠的三维电阻率分布。 相似文献
4.
波动方程数值模拟是研究地震波传播机理的重要工具,有限差分求解波动方程是当前地震波数值模拟的主要方法之一。当地下介质中的地震波速度较低或地震波高频成分丰富时,常规有限差分技术常常产生严重的数值频散误差,这种误差会降低数值模拟的精度,影响对地震波传播机理的分析。为压制地震波数值模拟时产生的数值频散误差,提高波场模拟精度,提出了基于NAD算子的时间四阶精度波动方程差分格式。根据对应的差分格式,分析了该差分格式的数值频散关系。与常规四阶精度差分算法的频散曲线相比,基于NAD时间四阶精度差分方法不但能够实现时间频散的有效压制,同时其基于更多网格点的位移分量和位移梯度分量空间微分求解方法还能够实现空间频散的有效压制。另外在相同模型条件下,基于NAD算法的声波方程时间四阶差分解法可采用大网格对模拟空间进行差分离散,减少网格数,提高计算效率。 相似文献
5.
起伏地形对可控源音频大地电磁(CSAMT)响应具有强烈的影响,因此在CSAMT数据处理解释时需要考虑地形。同时,实际的地下地质情况和地表的地形情况通常比较复杂,地质结构和地形大部分情况下都是三维的。在水平地表三维有限差分CSAMT数值模拟算法的基础上,推导了利用地下交错网格采样点处的总磁场计算起伏地形下空气-地下介质分界面处的总电场和总磁场的表达式,从而实现了起伏地形下三维CSAMT数值模拟算法。在算法实现过程中,采用伪δ函数代替麦克斯韦方程中的场源项和直接计算总场的策略,避免了原有的将总场分离成背景场和二次场的策略在复杂地质条件下难以选择合适背景电阻率的问题。为了直接模拟总场,起伏地形下三维CSAMT数值模拟算法给出了新的三维正演方程的边界条件。将模拟水平地表三维异常体和三维山峰地形两个理论模型得到的响应结果与前人算法的计算结果进行对比,验证了所实现算法的正确性和有效性。 相似文献
6.
在井间地震有限差分数值模拟中,用离散化的高阶差分方程近似连续导数的波动方程时,不可避免地会产生数值频散,而数值频散程度则直接影响到地震波数值模拟精度,因此为了得到清晰准确的地震波场记录,必须尽可能地压制数值频散。这里在一阶速度应力弹性波方程的基础上,利用两个约束条件构造拉格朗日函数获取优化差分系数,与泰勒展开差分系数下的交错网格高阶差分模拟结果比较,发现改进的优化交错网格差分算子的高阶差分数值模拟能更有效地压制数值频散,进一步提高交错网格高阶差分数值模拟的精度,为高精度井间地震数据的波场成像、纵横波联合解释等提供可靠依据。 相似文献
7.
殷文 《吉林大学学报(地球科学版)》2008,38(1):144-151
在弹性波频率空间域有限差分数值模拟方面,差分网格及边界条件是影响弹性波模拟成功与否的关键,为了压制数值模拟中的网格频散,采用25点有限差分算子,建立了有限差分矩阵方程,且借鉴匹配层衰减边界条件思想,设计了弹性波频率空间域有限差分数值模拟算法。由于采用高阶有限差分法来提高差分格式的精度,将会导致计算量显著增加,为此,对频率空间域有限差分弹性波数值模拟方法,采用流水线技术与分治策略进行了并行算法研究,提高了计算效率,使得在合理的计算时间内更精确地模拟弹性波在弹性介质中的传播过程。 相似文献
8.
依据电磁场理论，构建地电模型，利用有限差分对矿井三维地电异常体进行数值模拟，得到异常体在电阻率不同、至巷道距离不同、及测点不同情况下的层测深曲线，并对其分布规律进行分析研究。通过模拟矿井三维地质异常体地电异常响应特征，并与井下实测层测深曲线对比和分析，论证了矿井三维地电体数值模拟的正确性。 相似文献
9.
长接地导线源的一维解析解可以通过对电偶极源解沿导线方向积分获得,而三维模型中则难以直接加载。这里将长接地导线源看作多个电偶极源的组合,利用交错采样有限差分法,进行了直接计算电场总场的有限长接地导线源电磁场三维正演研究,并验证了算法的正确性。正演模拟时为消除场源奇异性,用电流密度呈伪delta函数分布的等效源来代替点源。 相似文献
10.
在基于波动方程的有限差分数值模拟中,会不可避免地出现数值频散（也称网格发散）问题。数值频散问题通常会给数值模拟的结果造成严重影响,因此在数值模拟中,应尽量设法消除这种现象。这里在前人的基础上,在基于各向同性介质的弹性波方程的数值模拟过程中,通过引入通量校正方法来解决数值频散问题。由数值模拟结果表明,该方法可以有效地消除数值频散现象,大大改善数值模拟的结果。 相似文献
11.
基于球坐标系下有限差分的地磁测深三维正演 总被引：2，自引：0，他引：2
为了计算全球尺度电磁感应的响应,本文介绍地磁测深频率域三维正演。正演算法采用球坐标系下的交错网格有限差分方法,从Maxwell方程的积分形式出发,采用PARDISO对离散后的方程组求解,避免了迭代求解的散度校正。为了验证本文结果的正确性和精度,与前人的有限元和有限差分方法进行了对比,一维层状模型的三维交错网格有限差分数值结果和解析解相对误差小于5%,双半球模型的计算结果与前人的计算结果完全吻合。三维"棋盘模型"计算表明磁场分量对异常体的大小和位置具有很好的分辨能力。 相似文献
12.
《物探化探计算技术》2017,(5)
为了消除源点的奇异性问题,在源点利用delta函数离散。设计一维的层状模型,验证海洋可控源2.5D有限元数值模拟算法的正确性和精度;考虑到在海中洋流、海底地形起伏变化大以及海底介质电性结构复的影响,海洋可控源电磁的发射源会发生倾斜、旋转、位置偏离测线,因此定义了任意取向的发射源。实际应用证明,发射源姿态变化以及不同发射频率和油气层参数,对海洋可控源电磁场响应都有很大地影响。 相似文献
13.
巷道影响下三维全空间瞬变电磁法响应特征 总被引：3，自引：0，他引：3
使用磁偶极源模拟三维均匀介质内部低阻体的全空间响应特征,采用非均匀网格和时间步长推导了瞬变电磁场的三维时域有限差分格式,并根据多个测点计算的数据,绘制了磁场和磁场对时间导数的等值线图,分析了瞬变电磁场的传播规律,以及巷道和交界面对场的影响特征。结果表明:在磁场穿过低阻体与大地的交界面时,磁场的扩散速度降低,衰减变慢;随着迭代次数的增加,磁场以低阻体为中心形成等值线环。巷道在早期对场影响较大,晚期则不明显,为解释全空间TEM异常提供参考。 相似文献
14.
姚伟华王鹏郭恒王信文安又新 《物探化探计算技术》2018,(3):353-357
从回线源瞬变电磁法三维有限差分正演算法出发,与均匀半空间解析解作比较,验证了数值模拟算法的正确性。利用该正演算法分别研究了巷道、背景地层电阻率、掘进机对矿井瞬变电磁超前探测响应的影响。结果表明,背景地层电阻率越低,响应幅值越强,巷道影响越显著;掘进机距探测线框超过10m时,影响幅度明显变小;矿井瞬变电磁法的方向性具体表现为在不同的探测角度上测量的异常幅度不同。通过研究这些影响因素,对分析矿井瞬变电磁超前探测资料具有重要的借鉴作用。 相似文献
15.
从回线源瞬变电磁法三维有限差分正演算法出发,与均匀半空间解析解作比较,验证了数值模拟算法的正确性。利用该正演算法分别研究了巷道、背景地层电阻率、掘进机对矿井瞬变电磁超前探测响应的影响。结果表明,背景地层电阻率越低,响应幅值越强,巷道影响越显著;掘进机距探测线框超过10m时,影响幅度明显变小;矿井瞬变电磁法的方向性具体表现为在不同的探测角度上测量的异常幅度不同。通过研究这些影响因素,对分析矿井瞬变电磁超前探测资料具有重要的借鉴作用。 相似文献
16.
17.
18.
《物探与化探》2015,(5)
采用三角网格剖分的有限元法,研究了2.5维航空瞬变电磁法正演模拟问题。利用时频变换数值方法将时间域电磁场转换到拉氏域,再利用傅里叶变换将三维问题降维变为2.5维问题,然后由有限元法求解得到拉氏域二维电磁场,逆拉氏变换后得到时间域航空瞬变响应。为了回避正演模拟中总感应磁场在场源处的奇异性问题,采用异常场算法,场源响应通过在微分方程中施加背景电磁场实现。由于瞬变电磁信号具有较大的动态范围,而且需要经过两次正、逆拉氏变换和傅里叶变换,每个环节的计算精度和速度要严格控制在较高的水平上,否则积累误差会非常大。模型计算表明均匀大地和层状大地模型解析解与数值解吻合很好。这证明该算法是正确可行的,可作为研究二维复杂地质体的方法手段。 相似文献
19.
在地震多发的中国内陆环境下,场地稳定性评价越来越引起重视。为了研究北京延庆新城规划区大型场地动力响应特征,在分析延庆新城规划区第四纪覆盖层和断裂的空间分布规律基础上,基于数值模拟方法,运用有限差分数值模拟软件FLAC3D建立了规划区三维仿真模型,并对新城规划区场地动力响应进行了仿真,最后通过定义场地速度响应系数来刻画场地动力响应特征。结果表明：对于“土岩组合”的典型场地,在低频剪切波的激励下,第四纪覆盖层的厚度将引起不同程度的场地速度响应,第四纪覆盖层越厚,场地速度响应系数越大,动力放大程度越高;对于下伏断裂的复杂场地,剪切波传播路径变得极其复杂,断裂对地表速度响应系数有放大作用,放大程度与断裂分布形式和发展规模有关。 相似文献
20.
本文以实例结合数值模拟说明CSAMT中非平面波效应对资料解释的重要影响。实际数据测量于河南卢氏县,在同一条测线测量,发射源长度1.3 km,收发距离分别为4.5 km和7 km。通过对数据进行必要的校正后进行了光滑模型反演,为了对比反演结果,整个处理过程的参数都保持一致。反演的电阻率断面显示,两种收发距离反演的电阻率断面其异常形态基本一致;但在电阻率的数值上,4.5 km的收发距离反演的电阻率较7 km的大。分析发现上述现象主要是源的非平面波效应造成的,并通过三维正演数值模拟得到了验证。由于在有限收发距离时,计算卡尼亚视电阻率靠近源受非平面波效应影响偏大,远离源后,该视电阻率趋于平面波视电阻率。因此,为了采用目前的可控源音频大地电磁测深解释理论,要求收发距离尽可能大。否则,必须采用带源效应的真正可控源电磁反演技术完成数据分析。 相似文献