计算水平层状VTI介质中P波旅行时的近似公式.pdf
《计算水平层状VTI介质中P波旅行时的近似公式.pdf》由会员分享,可在线阅读,更多相关《计算水平层状VTI介质中P波旅行时的近似公式.pdf(6页珍藏版)》请在咨信网上搜索。
1、2023 年 8 月第 58 卷 第 4 期计算水平层状VTI介质中P波旅行时的近似公式魏建,孙祥娥*(长江大学电子信息学院,湖北荆州 434023)摘要:基于泰勒级数的 P 波旅行时方程的计算精度受到具有垂直对称轴的横向各向同性(Transverse Isotropy Medium with Vertical Symmetry Axis,VTI)介质的影响,为了改善这种情况,提出一种基于平方处理与系数匹配的计算方法。首先,利用平方处理将非平方形式的基于泰勒级数的旅行时公式转换为包含高阶项的平方形式;再使用系数匹配法处理旅行时平方公式的高阶项,保持方程炮检距的最高阶为常用的 4阶;然后,在简化
2、参数形式与优化旅行时平方公式计算的基础上,得到与各向异性参数相关的新系数 y;最后,形成了含有 y的基于平方处理与系数匹配的 4阶 P 波旅行时计算方法。基于水平层状 VTI介质模型的实验结果表明,与有理近似算法、三射线广义时差近似计算法和扩展广义时差计算法相比,该方法计算误差更小,同时远炮检距处的计算能力得到一定提高,从而给以旅行时方程为基础的应用提供了更多选择。关键词:系数匹配,平方处理,旅行时,泰勒级数,非均质性参数中图分类号:P631 文献标志码:A doi:10.13810/ki.issn.1000-7210.2023.04.006Approximate equation for c
3、alculating Pwave traveltimein horizontally layered VTI mediaWEI Jian,SUN Xiang e(School of Electronics and Information,Yangtze University,Jingzhou,Hubei 434023,China)Abstract:The calculation accuracy of Pwave traveltime equations based on the Taylor series is affected by the transverse isotropy medi
4、um with vertical symmetry axis(VTI).In order to improve this situation,a calculation method based on square processing and coefficient matching is proposed.In this method,the square processing is used to transform the traveltime equation of the nonsquare form based on the Taylor series into the squa
5、re form containing highorder terms,and then the coefficient matching method is used to deal with the highorder terms of the traveltime square equation so that the highest order of the equation offset is the commonly used fourth order.Then,on the basis of simplifying the form of parameters and optimi
6、zing the traveltime square equation,a new coefficient y,which is related to the anisotropic parameter,is obtained.Finally,a fourthorder Pwave traveltime calculation method with y based on square processing and coefficient matching is formed.The experimental results based on the horizontally layered
7、VTI medium model show that compared with the rational approximation method,the generalized moveout approximation method based on three rays,and the extended generalized moveout approximation method,this method has smaller errors,and the calculation ability at the far offset is improved to some exten
8、t,which provides more options for the applications based on traveltime equations.Keywords:coefficient matching,squaring processing,traveltime,Taylor series,heterogeneity parameter魏建,孙祥娥.计算水平层状 VTI介质中 P波旅行时的近似公式 J.石油地球物理勘探,2023,58(4):812817.WEI Jian,SUN Xiang e.Approximate equation for calculating Pw
9、ave traveltime in horizontally layered VTI media J.Oil Geophysical Prospecting,2023,58(4):812817.0引言旅行时近似是时差校正、速度分析、反演等地学领域的重要研究课题之一13。炮检距函数应用于具有垂直对称轴的横向各向同性(Transverse Isotropy Medium with Vertical Symmetry Axis,VTI)介质时,常使用非双曲线形式。在早期研究中,学者们提出了多种适用于VTI介质的 P波旅行时近似方程。Taner等4 建议在旅行时公式中加入泰勒级数项得到非双曲线型公式;
10、May等5 提出一种基于正交多项式的非双曲线近似方法计 处理技术 文章编号:1000-7210(2023)04-0812-06*湖北省荆州市南环路 1号长江大学电子信息学院,434023。Email: 本文于 2022年 7月 17日收到,最终修改稿于 2023年 5月 10日收到。本项研究受大庆油田有限责任公司勘探开发研究院项目“井中 VSP 地震资料保真处理技术流程及质控方法研究”(DQYT12010022018 JS647)资助。第 58 卷 第 4 期魏建,等:计算水平层状 VTI介质中 P波旅行时的近似公式算四阶时差谱。由于VTI介质的影响,这些方法计算的旅行时会随炮检距的变大而出现
11、偏差。随后,学者们提出了其他基于泰勒级数的方法提高远炮检距旅行时的计算精度。Alkhalifah等6 利用声学近似理论提出了基于泰勒级数的 4 阶旅行时近似方程;胡中平7 提出一种优化炮检距 6阶项的计算方法;Stovas等8 在泰勒展开式的基础上提出了三种旅行时近似公式;之后,Ursin等9 利用泰勒级数高阶项进行系数匹配得到了4阶旅行时计算方程。以上几种方法使用的均是高阶泰勒级数的截断形式,远炮检距处旅行时的计算精度并不理想10。通过进一步分析旅行时公式在各向异性介质中的表现,学者们提出了基于泰勒级数与其他理论相结合的 P波旅行时计算方程。邓怀群等11 研究了适用于弱各向异性 VTI介质中
12、的旅行时计算方法;之后,Song等12 利用基于高阶泰勒级数系数的 Pad近似方法估计远炮检距时差;Sripanich等13 通过特定参数选择提出了简化的旅行时近似方程;Xu等14 基于泰勒级数的炮检距扩展和弹性正交介质中的声学近似优化了 P 波旅行时计算公式;另外,Abedi等15 与 Blias16 提出了基于广义时差近似的旅行时计算方法;随后,Abedi 等17 利用微扰理论推导出新的旅行时近似方程;2020年 Abedi18 提出了 P波运动学有理近似方法;之后,Farra等19 推导了基于弱各向异性参数的非双曲线型 P波时差方程;2021年Abedi等20 又提出了一种基于声学假设的
13、六参数旅行时计算方程;之后,一种利用降阶补偿的处理方法开始用于计算介质中的旅行时21。随着旅行时公式的进一步完善,能够提高旅行时计算精度的方法的种类也随之增加。针对常用的基于泰勒级数的 P波旅行时近似方程在水平层状VTI介质中的计算有待提高的问题,本文先利用平方处理将其转换为旅行时平方近似方程,再结合系数匹配方法依次处理炮检距的高阶项,由此得到一种含有与各向异性相关的新系数y的P波旅行时计算方法。经数据测试与方法对比表明,该方法可以改善远炮检距处旅行时的计算且误差较小。1理论方法1.1基于泰勒级数的炮检距 4阶旅行时公式Ursin 等9 在常规基于泰勒级数的旅行时公式基础上,提出一种基于系数匹
14、配方法的P波旅行时计算方程t(x)=t0+x22v2NMOt0+dx4v4NMOt30()1+Dx2v2NMOt20(1)其中d=-S28(2)D=2S22-S32S2(3)式中:t0为垂直双程旅行时;x为炮检距;vNMO为正常时差速度;S2、S3为非均质性参数,其中 S2与各向异性参数的关系为9,2223 S2=1+8(4)通过推导出的 S2与 S3的近似关系可进一号简化式(3)。首先写出基于泰勒级数定义的垂直慢度的系数aj(j=0,1,2,)9,22 a0=2a1=2(-)1-20(1-20+2)a2=4(-)(1-20)2(1-20+2)-(1+20)(5)式中:0=00表示S波与P波的
15、速度之比22;、0和0为 Thomsen参数6,9,其中0和0分别表示 P波与S波的垂直速度,与分别表示P波各向异性差异和NMO速度与反射振幅的变化6,9,23。S2、S3与系数aj的关系为S2=1+4a1(1+a0)2S3=1+4a1(1+a0)2+8a2(1+a0)3(6)结合声学近似理论22,即0=0,此时非均质性参数S3可用S2表示S3=1+8+322=1+S222(7)再将式(7)代入式(3),即可得到参数 D 的近似形式D=3S22-14S2=3(1+8)2-14(1+8)(8)813石 油 地 球 物 理 勘 探2023 年1.2基于平方处理与系数匹配的旅行时平方公式与旅行时平方
16、公式相比,式(1)是一种不常作为研究对象的旅行时方程。首先,通过将式(1)两边平方将其转换为平方形式。然后,按炮检距阶数从小到大排列t2=t20+x2v2NMO+Ax4+Bx6+Cx8(9)其中A=1+8d4v4NMOt20()1+Dx2v2NMOt20B=D+4d4v6NMOt40()1+Dx2v2NMOt20C=d2v8NMOt60()1+Dx2v2NMOt202(10)由式(9)和式(10)可见,三个高阶项系数 A、B和C均包含相对复杂的结构。为了保持相对简洁的结构以及体现高阶泰勒级数的意义,根据式(1)的组成形式,可使式(9)中炮检距的最高阶数仍为常规形式的四阶。在推导基于泰勒级数的旅
17、行时公式过程中,常对高阶项采用近似处理来简化方程形式。常用的处理方式有:省略高阶项的截断近似;利用高阶项进行系数匹配10 的近似。分析式(9)可知,截断近似并不能使其满足应用要求,而系数匹配近似是更适合的处理方法。首先,将系数匹配方法用于处理式(9)的 8阶项,即利用B和C可以得到t2=t20+x2v2NMO+Ax4+(D+4d)x64v6NMOt40()1+Dx2v2NMOt20(11)式中=CBv2NMOt20()1+Dx2v2NMOt20=4d2D+4d=S324(S22-1)(12)在式(11)的基础上,继续利用系数匹配方法处理6次方项。为了简化旅行时公式的参数结构,可先得到中间参数,
18、即=-BAv2NMOt20=S22-14S2(S2-1)=1+S24S2(13)结合式(13)即可得到最终形式的计算方程,即基于新系数y的炮检距四阶旅行时公式t2=t20+x2v2NMO-8x4v2NMO()4v2NMOt20+yDx2(14)其中y=-1=(1+S2)(S22-1)S42(15)基于平方处理与系数匹配方法得到的旅行时平方方程(式(14),可定义为 SCM(Squaring and Coefficient Matching)法。2实验结果从实验角度分析SCM与有理近似法(TR)18、三射线广义时差近似法(TG)24 和扩展广义时差近似法(TE)1,20 等 4种方法与精确时间的
19、计算误差,实验使用水平层状VTI介质模型,参数如表1所示。介质中P波的准确时间可由 x(p)=iux,iwit(p)=iut,iwi(16)计算1,20。式中ux,i=pv2NMO,iut,i=(1-2ip2v2NMO,i)2+2ip4v4NMO,iwi=t0(i)(1-2ip2v2NMO,i)3(1-p2v2h,i)(17)x(p)、t(p)分别表示通过射线参数 p计算出的炮检距与旅行时;t0,i、vNMO,i、vh,i、i分别表示第 i层的双程旅行时、正常时差速度、水平速度、各向异性参数。定义实验误差为TMO=|Ta-TxTa|(18)式中:Ta表示准确时间;Tx表示通过 TR、TG、TE
20、 和SCM法所得计算值。2.1模型一模型一为椭圆各向异性介质模型。图 1是本文方法(SCM)、TR、TG、TE法的旅行时计算误差曲线。由图可见,当炮检距与深度的比值为02.5时,SCM814第 58 卷 第 4 期魏建,等:计算水平层状 VTI介质中 P波旅行时的近似公式法的计算精度高于其他 3 种方法。在图 2 中,TR、TG、TE与SCM法的TMO之差均为正值,可知SCM法的计算误差最小。TE法作为公认的一种较准确的计算方法,其误差仅小于SCM法。虽然TG与TR法的 TMO 差值非常小,实际上 TR法仍具有更小的误差。从数值上分析表明,当炮检距与深度的比值等于2.5时,即远炮检距处,SCM
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 计算 水平 层状 VTI 介质 旅行 近似 公式
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。