排序方式： 共有70条查询结果，搜索用时 15 毫秒
31.
本文研究声电效应测井波场的有限差分模拟算法.忽略井外地层中诱导电磁场对孔隙弹性波的影响,将求解动电耦合波方程组的问题解耦,先计算孔隙弹性波,再计算其诱导电磁场.基于轴对称柱坐标系下的速度-应力交错网格,采用时域有限差分计算井孔流体声波和井外地层孔隙弹性波.将电磁场近似看作似稳场,基于轴对称柱坐标系下的5点式有限差分网格,求解不同时刻的电位Poisson方程,计算诱导电场.结果表明:本文算法可准确模拟频率6.0 kHz的声电效应测井全波;在声波测井频率范围内,电导率、动电耦合系数和动态渗透率的低频近似对伴随电磁场的计算影响不大;地层水平界面导致伴随反射斯通利波的电场和显著的界面电磁波,后者对于探测地层界面具有潜在的应用价值. 相似文献
32.
本文针对声-弹耦合介质,为尽可能的减少频率域正演模拟的计算内存,提高计算效率,在一阶非均质位移-应力波动方程的基础上,借助等效交错网格思想并充分考虑密度参数空间变化对地震波传播的影响,推导了声-弹耦合地震波波动方程.在流相介质和固相介质中分别采用非均质情况频率域二阶声压标量波、二阶纯位移控制方程,为保证流、固相介质间地震波能量的稳定传输和有效交换,提出了声-弹耦合界面转换过渡层方法,并详细阐述了过渡层与上下介质空间差分具体耦合方法.在与非均质纯位移波动方程正演结果对比分析的基础上,首先采用各向同性单层流相介质模型进行正演模拟验证了声-弹耦合方程数值模拟中过渡层策略的有效性和准确性,随后又数值模拟了地震波在声-弹耦合介质简单模型和复杂Marmousi2模型中的传播,验证了本文方法稳定性和准确性,同时该方法可以简单的推广到三维情况. 相似文献
33.
基于应力、速度混合变量弹性波方程及任意四边形网格差分算子,给出了交错计算应力及速度的非规则网格弹性波应力一速度差分法该方法融合了有限元法能适应复杂形状边界及差分法无需计算刚度阵的特点,具有较高的计算精度,所需计算机存储空间较少,计算效率也很高．基于积分平衡方程引入了任意形状自由表面的边界条件,且通过局部滤波改善了自由表面边界条件的稳定性,使得该方法可应用于考虑地表形状影响的地震波数值模拟 相似文献
34.
35.
非均匀介质一般为多尺度介质,对非均匀介质进行正演模拟需要多尺度的网格剖分.碳酸盐岩缝洞介质的尺度一般为厘米级甚至毫米级.使用有限差分方法对其进行精细模拟需要差分步长达到缝洞介质的尺度.为了提高有限差分数值模拟方法的精度和效率,使之可以应用于非均匀介质的正演模拟,本文推导了基于PML边界的空间和时间步长同时变化的高倍数可变网格差分格式,步长变化倍数可以达到百倍以上.并且在一般意义的变网格算法的基础上,改进了变网格算法的网格剖分方式,进一步减小了精细尺度模型数值模拟的内存消耗.数值试验表明,该方法可以精细描述毫米尺度的地质体,提高有限差分方法模拟精度,同时也节约了内存,提高了模拟效率. 相似文献
36.
2013年6月2日台湾省南投县发生一次里氏6.5级地震,这次地震是最近几年内发生的最大的一次地震,台湾大部分地区以及中国大陆东南沿海城市均有明显震感.为了更深入的理解此次地震的地震波能量传播过程以及强地面运动过程,本文采用三维有限差分方法对此次地震的强地面运动过程进行了模拟.结果显示,台湾中央山脉起伏的地形对此次地震的强地面运动分布特征具有较大的影响,出现了比较明显的地形放大效应;此外,台湾的平原和盆地中的沉积层也对地震波有较强的放大效应. 相似文献
37.
有限差分法广泛应用于地震波数值模拟、成像和波形反演中,差分数值解的精度直接影响着地震成像和反演的效果.因为有限差分算子可以通过截断伪谱法的空间褶积序列得到,而截断窗函数的属性影响有限差分算子逼近微分算子的精度.具体地讲,窗函数的幅值响应的主瓣和旁瓣决定了有限差分算子逼近的精度,主瓣越窄,旁瓣衰减越大,则有限差分算子逼近微分算子的精度越高,更好地压制数值频散.基于此认识,本文提出了一种基于Chebyshev自褶积组合窗截断逼近的有限差分算子优化方法.Chebyshev自褶积组合窗的主瓣较窄,且旁瓣衰减大,其可通过只调节三个参数,更直观和可视化地控制主瓣和旁瓣的形状,改变有限差分算子逼近微分算子的精度;该窗函数截断逼近的有限差分算子不仅有较大的谱覆盖范围,而且精度误差波动较小,这表明低阶的差分算子可以达到高阶算子的精度,且逼近误差更稳定;从经济上来讲,将有效地减少模拟计算花费,提高计算效率. 相似文献
38.
本研究实现了一套基于有限差分（FD）方法的大地电磁测深数据带地形三维反演算法及代码.其中,在大地电磁场正演数值模拟方面,开发了起伏地形条件下基于交错网格剖分、有限差分方法的大地电磁测深三维正演代码;在满足平面波场假设的前提下,使用长方体网格剖分模拟三维起伏地形,实现了带地形三维正演计算;并设计理论模型进行试算,经试算结果与前人的有限元法计算结果对比,验证了所研发的带地形三维正演计算的正确性与可靠性.在反演方面,本研究基于非线性共轭梯度方法编写了大地电磁测深带地形三维反演代码,试验了不同的共轭梯度搜索因子β,避免了目标函数对海森矩阵（参数二次导数矩阵）的显式计算和存储,初步实现了大地电磁资料的带地形三维反演.最后,对一系列理论模型进行正演计算,利用其生成的合成数据模拟实测数据进行反演,并与现有的不带地形大地电磁测深三维反演结果比较,检验了所研发的带地形三维反演计算的可靠性与稳定性. 相似文献
39.
气候和天气数值预测本质上都需建立在离散空间的基础上，并要求有高精度的空间离散化计算方案。然而，解析性微分算子和离散性差分算子的差异却是造成数值计算的原始误差源。作根据约束性代价函数原理提出了能够使得截断误差最小化的中心差分优化算子。从最优化角度看，该算子本质上是以一种简单但却相当创新的方式尝试应用变分连续同化（VCA）技术来最小化大气模式中的空间截断误差。此外，该优化差分算子的设计不依赖于网格结构，即可应用在绝大多数Arakawa网格上，诸如非跳点网格（Arakawa-A)以及常规跳点网格（Arakawa-B,Arakawa-C,Arakawa-D)。但为了对优化差分算子精确性进行严格测定，该算子的基准试验结果是通过在截断误差最易被激发出来的非跳点网格上进行严格理论和实际数值模拟测试得到的。两的结果都表明该优化差分算子截断误差造成的数字噪音能被减小到最小量。 相似文献
40.
A numerical experiment carried out to investigate the structural model of the Domenico lignite site is discussed. The model is a 2D structure containing several lignite layers at different depths, and a low-velocity layer at the top of the model. The experiment consists in simulating a measured CDP section by two independent techniques, based on completely different concepts: the finite-difference method and the ray method. Due to the incompleteness of the ray synthetic wave field, as well as to numerical problems of the finite differences at higher frequencies, the agreement between the synthetic seismogram sections for the individual shot points is poor. However, the CDP stacked sections modelled by the ray and finite-difference methods agree rather well. This is because the main differences between the wave fields computed by the two methods are due to the presence of the low-velocity layer (ground roll, head waves, etc.), and just these parts of the wave field can be suppressed by routine data processing such as f–k filtration. Synthetic ray and finite-difference CDP stacks agree relatively well with the observed data. They confirm three lignite seams and a fault in the shallower one. The synthetic data also indicate that many apparent horizons of the measured section may be due to the multiple reflections within the subsurface low-velocity layer. 相似文献