- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
针对裂缝储层叠前地震正演模拟方法研究
针对裂缝储层叠前地震正演模拟方法研究
针对裂缝储层的叠前地震正演模拟技术是裂缝储层地震预测的先导。有效地描述裂缝储层取决于准确地确定裂缝的方位分布和连通性以及与储层内非均匀体的空间关系。裂缝储层地震响应除了与裂缝密度、形状、饱和流体有关,还与岩石本身的岩性有关。裂缝会产生地震物性的各向异性,同时也会引起地震波传播的各向异性。在储层裂缝物理模型试验中,主要研究裂缝属性同纵横波属性(横波分裂、P波方位等)的关系。裂缝诱导的各向异性对介质中横波的传播产生了明显的影响,主要是引起横波分裂现象,并且横波各向异性的程度对着裂缝密度的增加而增大。在裂缝研究中,针对裂缝的叠前地震正演模拟技术是非常重要的,它可以帮助了解由于裂缝影响引起的地震波振幅、相位、频率随偏移距变化的规律。鉴于裂缝储层的复杂性,有必要设计一系列裂缝以及裂缝组合,分析并总结其地震响应特征,为叠前属性预测提供依据,并进一步指导裂缝预测。
1 基本原理
波动方程有限差分法正演模拟,对认识地震波传播规律、进行地震属性研究、地震资料地质解释、储层评价等,均具有重要的理论和实际意义。随着地震勘探的不断发展,油气勘探的重点正转向复杂地表条件和复杂地质条件的区域,希望有更准确的模拟方法来认识复杂介质的地震波传播规律,以便更好的地震勘探的处理和解释服务。针对这个问题我们拟提出了高精度差分数值模拟方法的研究与应用,以便更好的适应于当前地震勘探的需要。
1.1 高阶差分方法
用声波方程或弹性波方程进行有限差分法正演模拟时,如果使用二阶差分格式,网格间距必须取得很小,才能保证计算精度及稳定性。网格距的增大可以大大降低对计算机内存的要求、缩短计算时间。可实用化的三维地震波模拟方法为地震数据处理结果的地质解释的正确性提供了一种很好的检验手段。
1.2 交错网格技术
在交错网格技术中,变量的导数是在相应的变量网格点之间的半程上计算的。用下式计算一阶空间导数,可得到交错网格不同差分精度的差分权系数值。
1.3 完全匹配层吸收边界(PML)
??波方程完全匹配层吸收边界,其基本思想是在所研究区域的边界上引入吸收层,波由研究区域边界传到吸收层时不产生任何反射,在吸收层内按传播距离的指数规律衰减,不产生反射,从而达到吸收边界的效果。
1.4 频散校正―通量传输校正方法(FCT)
FCT方法是基于守恒型方程的差分格式。守恒型差分格式满足物理量守恒定律,这种格式可用控制体积法直接从物理守恒定律构造,也可以从散度型偏微分方程构造,在构造这种差分格式时,要求网格间交界面处守恒物理量的通量互相抵消,即交界面两侧一边流入和另一边流出的通量相等。这样内部单元边界上的截断误差被抵消,避免误差的积累。在交错网格有限差分法中,要求在网格边界上有物理量流入和流出,但要求其满足能量守恒关系,因此当两个网格的差分方程相加时,流过共同边界上的通量可相互消掉,所得的差分方程就和能量守恒关系直接应用于组合网格得到的差分方程相同。守恒格式保证差分方程与物理定律相互协调,可以很容易的处理交错网格的差分频散校正。下面给出FCT技术的具体实现过程。(1) 根据差分方程、初值和人为边界条件算出,其中,;(2)扩散计算。计算第n时间步的扩散通量,循环计算便可获得FCT后的差分结果。FCT方法的本质是对色散输运格式施加一个扩散项,使之能平滑通过冲击波区或流场陡变区且无数值振荡。反扩散计算则使得扩散误差得到补偿,然而在冲击波区或流场陡变区它却引起了振荡,这是因为反扩散之后出现了新的极值。采用限制反扩散通量条件后,消除了新极值,也消除了振荡,获得了稳定解。而且,修正扩散以守恒方式施加,即在某一网格点减去某个小量时,同样大小的小量会被加到另外一个网格点,这相当于一个小的修正值从一点挤向另外一点,从而不会导致整个系统的某个量的增加或减少。
2 正演模拟试算
为了对比裂缝参数不同情况下地震响应,我们建立了如下图所示的一系列速度场,取定上覆岩层纵波速度为3000m/s,密度为2.5下层纵波速度为4500m/s,密度为2.6裂缝内纵波速度为2300m/s,密度2.36,地层分界面位于第750个GDP点。根据模型,我们进行快了单炮正演。
2.1 不同裂缝组合。裂缝数量越多,其响应越强。
2.2 不同裂缝开度。用(A)(B)(C)三图分别代表在开度50cm、开度1m、开度2m情况下的地质模型。其延伸10m,密度20m,建立地质模型下的单炮地震记录,共1001道,采样率2ms,(A)(B)(C)分别代表开度50cm单炮记录、开度1m单炮记录、开度2m单炮记录,红色和黄色部分为裂缝响应,裂缝开度越大,其响应越强。
2.3 不同裂缝倾角角度。 (A)(B)(C)三图分别代表在倾角3
文档评论(0)