共查询到19条相似文献，搜索用时 171 毫秒
1.
本文将离散的局部透射边界与集中质量有限元方法相结合,用于模拟粘弹性成层介质中的稳态波动。并以一个受集中线力作用的粘弹性土层的出平面稳态波动为例,将有限元模拟与解析方法的计算结果进行比较,验证此模拟方法的有效性,着重显示在频域中实现高阶透射的可能性。 相似文献
2.
据逆时偏移的基本思想,首次引用傅氏变换和差分法相结合的方法(FDM)数值求解全声波方程,实现对垂直地震剖面(VSP)记录数据的波场外推。同时,利用射线方程实现对震源的向下延拓。从而解决了非零井源距VSP数据的偏移。 所用方法对地下界面倾角较大、速度横向变化没有限制,因而适应于地下构造比较复杂的地质情况。此外,还可用于地面地震记录的迭前深度偏移。 相似文献
3.
井外爆炸点源激发井孔声波场 总被引：1，自引：0，他引：1
利用井外爆炸点源激发的声波场与井孔的相互作用模拟井间地震的物理过程,推导了全频域井孔内外声波场的理论计算公式,并针对不同震源主频、不同源距(震源到井轴的径向距离)和软硬地层进行了井内声场的全波模拟与频率-相速度分析。结果表明:井外爆炸点源激发井孔内外的声场等效于无穷多级声场迭加的结果;震源靠近井壁时的频率-相速度图与震源在井内时相似;震源的主频越高、源距和接收器离井轴的径向距离越大、地层越硬,则顾及的多级声场的项数就要越多;震源主频越低、源距越小,井中激发的斯通利波则越强。 相似文献
4.
5.
AVO简化方程的物理意义及其在油气识别中的应用 总被引：3，自引：2，他引：3
地震反射振幅随炮检距变化研究(AVO)的基本原理是反射系数随入射角的变化而变化。这里从精确的Zoeppritz方程出发，根据Shuey及Hilterman提出的P—P波反射系数计算公式，设计多层地质模型进行AVO正演模拟并合成地震记录，利用地震记录对比讨论了Shuey及Hilterman的简化方程各项的物理意义及其在不同入射角范围内对地震记录的贡献。 相似文献
6.
井间地震弹性波传播特征数值模拟 总被引：1，自引：0，他引：1
采用弹性波方程数值模拟和Zoeppritz方程解相结合的方式,研究了井间地震固-固界面上以及低速和高速夹层中弹性波场特征及传播规律.给出了井间地震弹性波的主要类型,其中包含一种特殊面波,该面波由激发井和接收井内的地层界面处产生的二次源形成,并说明了反射波以广角反射占优势;利用正演模拟记录对广角反射现象进行了分析,得出了在临界角处广角反射波的振幅、相位突变,极性反转,这与反射系数变为复数相对应;提出利用临界反射波特征和临界反射点位置来确定地层速度.震源在高速夹层中激发时,其中传播的弹性波场简单;而在低速夹层中,则弹性波场非常复杂.低速夹层中的导波存在明显的频散现象,其振幅和频率与低速层的速度和厚度密切相关.提出利用导波可以确定低速层的横向变化和连通性. 相似文献
7.
8.
多个激励源无延时发射(多源并发)相同中心频率脉冲会形成平面波束信号,增强数据记录质量。本文通过数值模拟对比分析在多源并发情况下,非分裂递归卷积完全匹配层作为吸收边界条件和Mur二阶吸收边界条件对电磁波的吸收效果。其研究结果表明,传统的Mur二阶吸收边界条件对多源并发、多角度掠射情况下电磁波的吸收效果不佳,在大偏移距下会造成波形畸变和形成虚假反射。而在多源并发情况下采用非分裂递归卷积完全匹配层作为吸收边界条件,将坐标伸缩因子引进时域有限差分算法中。通过傅里叶逆变换将频率域坐标伸缩变换PML方程转换到时域,对电场和磁场值在离散状态下进行递归卷积运算求解。从而避免了直接对卷积进行数值求解的复杂计算,在保证计算准确性的同时,节约了内存空间,提高了计算效率。在不分裂波场情况下,改善了网格截断位置对电磁波的吸收效果。 相似文献
9.
对遭遇浅层气灾害地质的工程,一般防治措施多为建设前期布设超前排气井进行控制性排气。深入了解排气过程中气藏内部的气水运移规律,对于排气井井距的设计和高效率排气有重要意义。基于非饱和土两相渗流原理和杭州湾富浅层气区的典型概化地质模型,利用数值模拟方法研究了有控排气条件下气藏内部气水界面运移过程,探讨了影响超前排气井布设的因素,提出一种井距优化方法,并在实际工程中应用。结果表明：不同井口排气速率下,气藏内部气水分布形态与气体排放效果显著不同;工程超前排气井井距的设置可以淹井时间最长和剩余气藏厚度最小进行优化;原始气藏压力是影响井距优化的最显著因素,而储气层下部渗水系数、气藏厚度等并不影响最优井距的确定;调节井距对改善竖向排气井的排气效果随原始气藏压力系数减小而减弱;原始气藏压力系数大,井距适当减小,反之,可适度增加;对于气藏压力系数小的工程,建议发展水平排气井布设技术,以改善排气效果。此外,工程超前排气井的井距设置还应考虑与建(构)筑物的合理安全距离,达到高效、安全和避免投资浪费的目的。 相似文献
10.
11.
正交各向异性介质中地震记录的有限差分数值模拟 总被引：9，自引：3，他引：6
本文首先介绍了正交各向异性介质弹性波动方程的建立过程,然后讨论了高阶有限差分算法的实现原理、算法的稳定性和收敛性以及吸收边界条件的选取,并对正交各向异性介质作了正演模拟。通过波场快照、共炮点地震记录和ＶＳＰ模拟记录等正演分析,结果表明该算法在模拟各向异性介质中弹性波的传播是有效的 相似文献
12.
分析了复杂地表地区复杂地震波场的形成原因,探讨了基于高角度时空域单程波动方程适用于起伏地表之下构造成像的逆时叠前深度偏移方法。从声波方程有限差分法合成起伏地表上的炮集记录知,当激发点位于低速层中,共炮点道集中低速层表面和高速层表面接收到的信号强度相差较大:低速层表面的信号很强,而高速层表面的信号很弱,低速层中产生了较强的槽波。将该逆时叠前深度偏移方法应用于起伏地表地震波场的偏移处理。虽然在偏移成像前既没有压制与低速层有关的槽波,又没有压制随机噪声,但偏移剖面上界面清晰、位置正确,断层面也得到了很好的聚焦。 相似文献
13.
14.
转换横波共炮记录单程声波方程正演模拟 总被引：7，自引：1，他引：6
本文分别把分离的纵波和横波方程分裂为各自的两个单向旁轴方程（Ｊ．ＦＣｌｅａｒ－ｂｏｕｔ）,然后用单程纵波和横波方程,结合相位移加有限差分混合波场延拓法（Ｄ．Ｒｉｓｔｏｗ）,在共炮记录叠前深度偏移原理的基础上,借鉴几何地震学的射线法两点问题原理,提出合成转换Ｐ—ＳＶ波共炮记录的时间一致性原理和逆过程方法,从而推导出精确的转换Ｐ—ＳＶ波共炮记录正演模拟算法,由于单程延拓波场,可完全解决多次反射波问题。 相似文献
15.
有限元法是复杂介质地震模拟的有力工具,它能比较客观地反映地震波传播过程,比较细致地再现地震图像。但由于有限元法地震模型所需网格节点数多、数据规模大,使它的计算受到限制。采用对角化集中质量矩阵、刚度矩阵变二维为一维存贮,并通过使用有限差分吸收边界条件来处理地质模型的人为边界反射问题,从而使计算所需内存降低,计算时间减少,并能比较好地实现地震声波用有限元法计算,理论模型的计算结果验证了算法的正确性和可行性。 相似文献
16.
17.
复合单元法动力分析初步研究 总被引：1，自引：0，他引：1
复合单元法是一种适用于不连续介质结构分析的数值计算方法,优点是建立网格时不用考虑结构中不连续面,计算时能离散模拟不连续面。在复合单元法静力计算模型的基础上,提出了复合单元的动力分析方法。给出了复合单元动力分析的质量矩阵和阻尼矩阵的具体表达式,通过质量矩阵、刚度矩阵和阻尼矩阵将分布于复合单元内部的惯性力和阻尼力分别转移到复合单元的广义节点上。建立了惯性力、阻尼力、动力荷载及弹性力作用下的复合单元动力系统的控制方程组。以一个基岩内含有两条断层的混凝土重力坝模型地震为例,采用复合单元法动力分析程序进行计算,其结果与有限单元法解吻合得较好,说明所提出方法的正确性和有效性。 相似文献
18.
在勘探地球物理领域中,能快速、有效地数值模拟2D和3D的地震全波场,包括同时模拟P波、S波、面波多分量波场特征的软件并不是很多;但在结构力学领域中,已有多种成熟的多波多分量有限元数值模拟大型商用软件用于求解弹性本构方程,如ANSYS,PLAXIS等。从方程结构来看,弹性本构方程与完全弹性波动方程相比,只多一项对时间的一阶微分项--阻尼项。因此,调整合理的介质参数--Rayleigh系数,令阻尼项近似为0,弹性本构方程就蜕化为完全弹性波动方程。隧道地震波场是在岩体中传播,具有弹性参数较好、无常规地震勘探中覆盖层等低速带干扰等的优点,采用ANSYS软件对隧道地震波场进行数值模拟,研究了频散、阻尼与吸收边界等数值模拟的相关技术。针对介质吸收,比较了数值模拟中有无阻尼的时间记录和波场快照,对隧道地震预报来说,将阻尼系数设定为0是一种合理的假设;比较了实例中不同网格长度与频散的关系,当网格长度小于波长的1/π倍时,才能消除频散;通过由波方程的推导,引入了黏弹性边界条件,通过实例计算证明它可以有效地吸收边界反射;最后对隧道地震预报进行了实例计算,算例表明采用ANSYS软件可以有效地模拟隧道复杂地质条件下全波场的激发传播过程。 相似文献
19.
根据白山抽水蓄能泵站地下厂房开挖过程中的变形观测数据,提出了一种基于响应面法的岩体力学参数反演方法。该方法利用响应面函数建立了岩体力学参数与围岩变形之间的非线性关系。通过有限元数值模拟确立了响应面函数中的系数。定义参数反演的目标函数,将参数反演问题转化为优化问题。分别采用拟牛顿优化算法和遗传算法求解参数反演的目标函数,得到了地下厂房的岩体力学参数。根据反演确定的岩体力学参数,对地下厂房围岩的开挖变形进行了数值模拟,研究表明,有限元模拟的地下厂房与现场观测值基本一致,验证了反演方法的有效性。 相似文献