横观各向同性场地中埋置衬砌隧道对qP波的散射          下载免费PDF全文

巴振宁  安东辉  梁建文  高旭 《地震工程学报》2020,42(4):948-954

采用间接边界元法(IBEM)结合"分区契合"技术,研究qP波入射下横观各向同性(TI)场地中衬砌隧道的动力响应问题。方法充分利用半空间和全空间动力格林函数在分别构造含孔半无限空间域和闭合域内散射波场方面的优势,将含有衬砌隧道的层状弹性半空间分解为含孔半无限空间域和一个环形衬砌闭合域来分别进行波场构造,有效地降低了求解时间和存储量。文中验证方法的正确性,并以均匀TI半空间和基岩上单一TI土层为例,计算分析弹性半空间场地中隧道衬砌内表面动应力放大问题。结果表明qP波入射下,TI介质与各向同性介质场地中埋置衬砌隧道的动力响应差异显著,TI介质参数的变化导致场地动力特性的改变,进而改变场地与衬砌隧道的动力相互作用机制,显著影响着衬砌内表面动应力的大小及其空间分布。  相似文献   

水平层状非均质TI地层中仪器偏心情况下三维感应测井响应高效数值模拟与响应特征分析          下载免费PDF全文

林蔺  杨守文  白彦  陈涛  汪宏年 《地球物理学报》2017,60(5):2000-2010

本文针对中国石油集团测井有限公司(CPL)三维阵列感应测井仪器(Three Dimensional Induction Tool,TDIT)的真实结构与参数,利用2.5维数值模式匹配技术研究建立水平层状非均质TI地层中仪器偏心情况下三维感应测井响应的高效数值模拟算法.首先,利用2.5维数值模式匹配技术给出水平层状非均质TI模型中偏心磁流源并矢Green函数的半解析解.在此基础上,结合三维阵列感应测井仪器参数以及测量过程中出现的仪器偏心与旋转等实际情况,通过叠加原理推导出该仪器响应的有效计算公式.最后,通过数值模拟结果分析考察仪器偏心、仪器自旋角等参数对测井响应的影响以及变化特征,为多分量感应资料处理与解释提供理论基础.  相似文献   

基于Low-rank分解的复杂TI介质纯qP波正演模拟与逆时偏移   总被引：2，自引：2，他引：0       下载免费PDF全文

黄金强  李振春 《地球物理学报》2017,60(2):704-721

近年来,面向实际应用的TI介质准P波正演模拟与逆时偏移成像技术受到空前的关注.基于常规耦合型传播方程的正演模拟方法不仅存在伪横波及频散假象干扰,而且还遭受模型参数限制(η0)和不稳定影响;而纯qP波方程的推导繁琐,且由于方程中包含拟微分算子造成求解难度大且精度有限.为此,本文首先构建了一种适用于任意TI介质的纯qP波传播算子,然后借助Low-rank分解求取该算子中的空间-波数域矩阵,同时引入Cerjan衰减边界条件来压制边界反射干扰,最终实现了一种间接的纯qP波波场外推方案,并将其成功应用于复杂TI介质正演模拟与逆时偏移成像中.通过开展数值模拟,并与其他方法对比表明:①该方法既避免了纯qP波方程的繁琐推导,又克服了耦合型方程对模型参数的限制;②还彻底消除了残余伪横波噪音及数值频散;③且能适应较大时间或空间步长及高频震源,是一种相对准确且稳定的各向异性纵波正演与成像策略.  相似文献   

任意空间取向TI弹性张量解析表述   总被引：6，自引：4，他引：2       下载免费PDF全文

姚陈  蔡明刚 《地球物理学报》2009,52(9):2345-2348

本文理论给出任意空间取向TI(ATI)四阶弹性张量的解析表述,其以VTI弹性常数及其简单组合为系数,包括各向同性项、TI对称轴方向矢量分量的二次项和四次项,其中TI对称轴方向矢量可以在固定坐标系定义, 也可以相对三维倾斜界面甚至相对波传播方向.相比四阶张量变换法和Bond变换法,ATI弹性张量能简洁而透明地为本构关系和波动方程提供四阶张量的所有元素. ATI弹性张量能为诸多方面的理论研究提供支撑.  相似文献   

TTI介质的相速度研究     

轩义华  何樵登  林炎 《应用地球物理》2007,4(1):25-28

推导了二维TTI介质的相速度表达式,并且依据推导出来的相速度表达式,模拟并分析了二维TTI介质相速度的传播快照以及TI介质相速度的传播快照;对比并分析了TTI介质和TI介质模型的相速度理论计算值的X分量特征的差异。TTI介质的相速度研究具有较高的理论研究价值和实际应用价值．  相似文献   

3D P‐wave traveltime computation in transversely isotropic media using layer‐by‐layer wavefront marching          下载免费PDF全文

Jiangtao Hu  Junxing Cao  Huazhong Wang  Xingjian Wang  Renfei Tian 《Geophysical Prospecting》2018,66(7):1303-1314

Subsurface rocks (e.g. shale) may induce seismic anisotropy, such as transverse isotropy. Traveltime computation is an essential component of depth imaging and tomography in transversely isotropic media. It is natural to compute the traveltime using the wavefront marching method. However, tracking the 3D wavefront is expensive, especially in anisotropic media. Besides, the wavefront marching method usually computes the traveltime using the eikonal equation. However, the anisotropic eikonal equation is highly non‐linear and it is challenging to solve. To address these issues, we present a layer‐by‐layer wavefront marching method to compute the P‐wave traveltime in 3D transversely isotropic media. To simplify the wavefront tracking, it uses the traveltime of the previous depth as the boundary condition to compute that of the next depth based on the wavefront marching. A strategy of traveltime computation is designed to guarantee the causality of wave propagation. To avoid solving the non‐linear eikonal equation, it updates traveltime along the expanding wavefront by Fermat's principle. To compute the traveltime using Fermat's principle, an approximate group velocity with high accuracy in transversely isotropic media is adopted to describe the ray propagation. Numerical examples on 3D vertical transverse isotropy and tilted transverse isotropy models show that the proposed method computes the traveltime with high accuracy. It can find applications in modelling and depth migration.  相似文献   

任意空间取向TI介质中长排列非双曲动校正公式优化(英文)     

郝重涛  姚陈  张建中 《应用地球物理》2010,7(4):336-347

常规长排列非双曲动校正公式是在VTI介质中得到的,它不能满足任意空间取向TI(ATI)条件下的扩展.本文以VTI介质中非双曲动校正公式为基础,基于我们推导得出的ATI介质中精确四次时差系数解析解和NMO速度解析解,给出ATI介质中长排列优化的非双曲动校正公式.通过与各向异性射线追踪方法计算所得出的"精确走时"结果对比,研究表明优化后的非双曲动校正公式能精确地描述任意强弱、ATI介质中随测线方位变化的走时曲线,可以用来替代耗时、多偏移距、多方位的射线追踪方法正演拟合ATI介质中长偏移距反射走时,为利用非双曲时距的各向异性参数反演提供理论基础性认识。  相似文献   

Generalized moveout approximation for P–SV converted waves in vertically inhomogeneous transversely isotropic media with a vertical symmetry axis          下载免费PDF全文

Qi Hao  Alexey Stovas 《Geophysical Prospecting》2016,64(6):1469-1482

  相似文献   

一种TI介质纯qP波正演方法及其在逆时偏移中的应用          下载免费PDF全文

杨鹏  李振春  谷丙洛 《地球物理学报》2017,60(11):4447-4467

基于Tsvankin提出的精确频散关系,利用近似展开的方法,推导出解耦合的TTI介质纯qP波近似方程,并将方程中的偏微分算子分解成一个laplace算子和一个标量算子,用于代表qP波的精确传播方向,构建时间域二阶纯qP波方程.此推导过程无需设置横波速度为零,能够更加精确地描述qP波的运动学特征.这个方程相比于求解波数域二阶解耦qP波方程,计算效率高,存储需求小;相比于基于Alkhalifah频散关系推导的时间域二阶纯qP波方程,假象干扰压制好,数值误差小,更具一般性.但此方法求解波矢量时采用波场梯度一阶渐近近似,会造成垂直于对称轴方向的波场振幅不准确.为了较正振幅,将椭圆分解方法应用于此方程中,构建纯qP波椭圆分解方程,使得振幅更加均衡,并与Xu等提出的方程比较分析,应用本文构建的纯qP波椭圆分解方程得到的波场振幅值更加准确.本文首先选取了均匀TI介质模型进行了qP波正演模拟,并抽取波场单道波形进行振幅分析,验证了本文构建的纯qP波方程和纯qP波椭圆分解方程的正确性及有效性;然后选取BP TTI模型进行了qP波正演模拟,将其qP波正演结果和均匀TI介质模型振幅分析结果相结合,突出了本文构建的纯qP波椭圆分解方程的优势及适应性;最后选取逆冲模型和BPTTI模型,应用本文构建的纯qP波椭圆分解方程对其进行逆时偏移成像,验证了本文构建的纯qP波椭圆分解方程在逆时偏移中的可行性和适用性.  相似文献   

基于LS-RSGFD方法优化的横向各向同性(TI)介质一阶qP波高精度数值模拟   总被引：2，自引：2，他引：0       下载免费PDF全文

李振春  杨富森  王小丹 《地球物理学报》2016,59(4):1477-1490

为克服各向异性介质弹性波数值模拟中存在着计算量大和波场分离困难等局限,研究了声学近似的VTI介质和TTI介质一阶qP波数值模拟方法.首先对VTI介质弹性波方程进行声学近似,推导了VTI介质一阶qP波方程;然后基于精确的TTI介质频散关系,引入一个包含各向异性控制参数σ的新辅助波场,推导了稳定的TTI介质二阶耦合qP波波动方程,并通过引入波场的伪速度分量,推导了等价的一阶应力-速度形式.结合旋转交错网格有限差分(RSGFD)和基于最小二乘优化的有限差分(LS-FD)两种各具优势的方法,研究了最小二乘旋转交错网格有限差分(LS-RSGFD)方法,并用其数值求解VTI和TTI介质一阶qP波方程,然后通过构造其LS-RSGFD格式,实现了高精度的各向异性介质qP波波场数值模拟.数值模拟结果表明:TI介质一阶qP波方程能够准确地模拟各向异性介质中qP波的运动学特征,引入控制参数σ能够有效地减弱不稳定性问题,保证非均匀TTI介质中qP波场的稳定传播;利用优化的LS-RSGFD方法可以得到高精度的合成地震记录,同时还可以相对地提高计算效率.  相似文献