可控震源粘弹性波动方程有限差分模拟

英文篇名：Finite difference simulation of viscoelastic wave equation in vibroseis
中文刊名：煤田地质与勘探
英文刊名：Coal Geology & Exploration
作者：钟飞 ; 张伟 ; 焦标强 ; 钟约先
英文作者：ZHONG Fei1 ; ZHANG Wei1 ; JIAO Biaoqiang2 ; ZHONG Yuexian1(1.Key Laboratory for Advanced Materials Processing Technology ; Ministry of Education ; Department of Mechanical Engineering ; Tsinghua University ; Beijing 100084 ; China ; 2.Locomotive and Car Research Institute ; China Academy of Railway Sciences ; Beijing 100081 ; China)
中文关键词：有限差分 ; 可控震源 ; 粘弹性 ; 正演模拟
英文关键词：finite difference ; vibroseis ; viscoelastic ; forward modeling
出版日期：2011-04-25
机构：清华大学机械系教育部先进成形制造重点实验室;中国铁道科学研究院机车车辆研究所;
年：2011
期：02
出版单位：煤田地质与勘探

摘要

在一般的地震波计算中,地球介质可以当作各向同性的完全弹性体来对待。而对地表土壤和精细观测,粘弹性介质模型比完全弹性模型更符合实际。可控震源是近年来发展起来的精确可控、安全有效的新型震源,但国内对可控震源传播动力学的研究较少。对可控震源激发的地震波在粘弹性介质中的传播,用有限差分方法进行了数值模拟,给出了Kelvin一阶粘弹性波动方程组、吸收边界条件和震源处理。提出了检波器接收信号加速度和传播距离的三次多项式指数关系式,比线性指数关系更准确。使用该公式可以预测采集系统有效接收距离,指导测线布置。
Earth can be treated as full elastic isotropic medium in general calculation of seismic wave.Viscoelastic model is more realistic than elastic model for fine observation of surface soil.Vibroseis is a precise,controllable and safe new source developed in recent years.The study of vibroseis transmission dynamics is less mature than pulse source.This article uses finite difference to simulate propagation of seismic wave excited by vibroseis in viscoelastic medium.Kelvin viscoelastic wave equations with absorbing boundary condition and vibration source model are given.The acceleration of received signals and the propagation distance have the cubic polynomial ex-ponential relationship,which is more accurate than the linear exponential relationship,this formula can predict the effective received distance of certain acquisition system.

引文

[1]王忠仁,陈祖斌,张林行,等.可控震源非线性扫描地震响应的数值模拟[J].地球物理学进展,2006,21(3):756-761.
    [2]柴治嫒.可控震源伪随机扫描方法与地震响应的数值模拟[D].吉林:吉林大学,2007:2-3.
    [3]吴国忱.各向异性介质地震波传播与成像[M].山东东营:中国石油大学出版社,2006:5-98.
    [4]ROBERT L,HIGDON.Numerical absorbing boundary condi-tions for the wave equation[J].Mathematics of Computation.1987,49(179):65-90.
    [5]RENAUT R A,PETERSEN J.Stability of wide-angle absorbing boundary conditions for the wave equation[J].Geophysics,1989,54(9):1153-1163.
    [6]谭未一,赵兵,张中杰.粘弹性VTI介质地震波模拟特征分析[J].球物理学进展.2007,22(4):1305-1311.
    [7]杜启振,王延光,付水华.方位各向异性粘弹性介质波场数值模拟[J].球物理学进展,2006,21(2):502-504.
    [8]张毅.声波正演模拟中的人工边界问题[J].工程地球物理学报,2007,4(4):360-365.
    [9]孙卫涛.弹性波动方程的有限差分数值方法[M].北京:清华大学出版社,2009.
    [10]牛滨华,孙春岩.半空间介质与地震波传播[M].北京:石油工业出版社,2002:87-170.