地震波各向异性走时层析成像有限元模型仿真
|
摘要
为解决当前地震波各向异性走时层析成像方法中存在的成像准确性差、耗时长的问题,提出基于有限元的地震波各向异性走时层析成像方法。利用地震波应力张量分量、地震波应变张量分量、矢量松弛子函数值分量等参数计算线性粘弹性介质的应力应变关系,并根据所得结果获取地震波各向异性介质本构方程的应力应变之间的关系,完成地震波本构关系的分析。依据分析结果,得到地震波三维各向异性粘弹介质运动方程。基于所得本构关系与运动方程,对地震波应变位移关系进行分析,并得到质量矩阵、阻尼矩阵、刚度矩阵和力矩阵等参量。利用所得参量设置切向边界和法向边界,将地震波各向异性走时有关参数输入至ANSYS软件中的Plane182单元,完成地震波各向异性走时层析成像。实验结果表明,该方法模拟地震波各向异性走时层析成像结果与实际结果拟合程度高,成像耗时短。
In order to solve the low accuracy of image formation in the current seismic anisotropy travel-time tomography method, this paper puts forward a method of anisotropic travel-time tomography method based on finite element. Firstly, the stress-strain relationship of a linear viscoelastic medium was calculated through some parameters such as the stress tensor component of seismic wave, the strain tensor component of seismic wave and the component of phasor relaxation subfunction, and then based on the obtained results, the relationship between stress and strain of constitutive equation of seismic wave anisotropic medium was obtained, so as to complete the analysis of seismic wave constitutive relationship. Based on the analysis result, the motion equation of three-dimensional anisotropic viscoelastic medium of seismic wave was obtained. On the basis of constitutive relation and the motion equation, the strain-displacement relation was analyzed and some parameters such as mass matrix, damping matrix, stiffness matrix and force matrix were obtained. Finally, the obtained parameters were used to set the tangential boundary and the normal boundary. Then, relevant parameters of anisotropy travel time of seismic wave were inputted to Plane182 unit in ANSYS software. Thus, the seismic anisotropy travel-time tomography was completed. Simulation results show that the simulation result of proposed method has high fitting degree with the actual results. Meanwhile, the image formation takes a short time.
In order to solve the low accuracy of image formation in the current seismic anisotropy travel-time tomography method, this paper puts forward a method of anisotropic travel-time tomography method based on finite element. Firstly, the stress-strain relationship of a linear viscoelastic medium was calculated through some parameters such as the stress tensor component of seismic wave, the strain tensor component of seismic wave and the component of phasor relaxation subfunction, and then based on the obtained results, the relationship between stress and strain of constitutive equation of seismic wave anisotropic medium was obtained, so as to complete the analysis of seismic wave constitutive relationship. Based on the analysis result, the motion equation of three-dimensional anisotropic viscoelastic medium of seismic wave was obtained. On the basis of constitutive relation and the motion equation, the strain-displacement relation was analyzed and some parameters such as mass matrix, damping matrix, stiffness matrix and force matrix were obtained. Finally, the obtained parameters were used to set the tangential boundary and the normal boundary. Then, relevant parameters of anisotropy travel time of seismic wave were inputted to Plane182 unit in ANSYS software. Thus, the seismic anisotropy travel-time tomography was completed. Simulation results show that the simulation result of proposed method has high fitting degree with the actual results. Meanwhile, the image formation takes a short time.
引文
[1] 孔寅,等. 地震波层析成像技术在冲击地压防治中的应用[J]. 煤炭工程, 2016,48(6): 83-86.
[2] 刘成禹,刘汗青,黄振. 基于实测地震波波速的围岩各向异性分析——以赣龙铁路汀州隧道为例[J]. 岩土力学, 2016,37(5): 1451-1457.
[3] 田原,等. 中国东北地区高精度相速度层析成像及其构造意义[J]. 地球物理学报, 2017,60(5):1659-1675.
[4] 黄国娇,白超英,钱卫. 一般TI介质中多震相旅行时层析成像——以井间成像为例[J]. 石油地球物理勘探, 2016,51(1):115-126.
[5] 王华忠,等. 特征波反演成像理论框架[J]. 石油物探, 2017,56(1):38-49.
[6] 杨修超,等. 利用全波形反演提高逆时偏移成像的精度[J]. 地球物理学进展, 2016,31(3):1257-1261.
[7] 蔡方伟,卞春江,冯水春,等.基于CUDA的实时红外辐射传输及成像仿真设计[J].电子设计工程,2018, 26(1):60-65.
[8] 翟桐立,等. 全方位高密度单点接收地震采集技术[J]. 石油学报, 2016,37(b12):56-63.
[9] 王书文,冯美华. 观测系统对回采工作面地震CT成像质量影响规律研究[J]. 中国煤炭, 2017,43(10):33-38.
[10] 郭强,等. 非理想情况下星载微波辐射计成像建模仿真[J]. 计算机仿真, 2017,34(5): 44-48.
[2] 刘成禹,刘汗青,黄振. 基于实测地震波波速的围岩各向异性分析——以赣龙铁路汀州隧道为例[J]. 岩土力学, 2016,37(5): 1451-1457.
[3] 田原,等. 中国东北地区高精度相速度层析成像及其构造意义[J]. 地球物理学报, 2017,60(5):1659-1675.
[4] 黄国娇,白超英,钱卫. 一般TI介质中多震相旅行时层析成像——以井间成像为例[J]. 石油地球物理勘探, 2016,51(1):115-126.
[5] 王华忠,等. 特征波反演成像理论框架[J]. 石油物探, 2017,56(1):38-49.
[6] 杨修超,等. 利用全波形反演提高逆时偏移成像的精度[J]. 地球物理学进展, 2016,31(3):1257-1261.
[7] 蔡方伟,卞春江,冯水春,等.基于CUDA的实时红外辐射传输及成像仿真设计[J].电子设计工程,2018, 26(1):60-65.
[8] 翟桐立,等. 全方位高密度单点接收地震采集技术[J]. 石油学报, 2016,37(b12):56-63.
[9] 王书文,冯美华. 观测系统对回采工作面地震CT成像质量影响规律研究[J]. 中国煤炭, 2017,43(10):33-38.
[10] 郭强,等. 非理想情况下星载微波辐射计成像建模仿真[J]. 计算机仿真, 2017,34(5): 44-48.