地震勘探中的三维随机介质模拟.docx

《地震勘探中的三维随机介质模拟.docx》由会员分享，可在线阅读，更多相关《地震勘探中的三维随机介质模拟.docx（9页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、地震勘探中的三维随机介质模拟摘要：常规的均匀层状介质理论已经不适用于高分辨率地震勘探，需要进一步发展随机介质理论，为研究复杂非均匀介质中的地震波场理论奠定基础。本文首先利用三维VonKarman型自相关函数模拟三维随机介质并给出了具体的建模流程，通过引入三维锥形函数，来压制建模时离散计算产生的误差，提高了模型的可信度。然后具体讨论了自相关长度、谱指数、标准差对建立的三维随机介质模型的影响。最后分析了模型内的速度变化和概率分布情况。模拟结果表明利用三维VonKarman型自相关函数能够更灵敏地模拟三维随机介质。关键词：三维随机介质；三维VonKarman型自相关函数；自相关长度；谱指数在地震勘探

2、中，长期以来都是以传统的均匀介质或层状均匀介质理论为基础，进而忽略了实际地震勘探中介质的局部非均匀性，但是随着高分辨率地震勘探的不断发展、勘探目的已由常规油气藏逐步转向复杂油气藏，介质中分布的大量微小异常会影响勘探结果1。用传统的介质模型理论很难准确地描绘这些异常，而以统计学理论为基础的随机介质模型能够较为灵敏、方便、有效地描绘这些介质的微小异常，进而便于研究局部非均匀介质的地震波场的变化特征2-4。国内外很多学者对随机介质理论已经进行了大量理论研究，并获得了一些成果。Ikelle等建立了二维随机介质模型并进行了正演模拟5-7。国内自2001年至今，姚姚等在这方面做出了很多成果，他们给出了随机

3、介质模型的建立方法、并对各种随机介质模型进行了波动方程正演模拟以及波场特征分析等方面进行了深化的研究8-19。本文利用三维VonKarman型自相关函数建立了三维随机介质模型，引入三维锥形函数来压制离散误差，并对自相关长度、谱指数、方差三个模型参数对三维随机介质模拟的结果和特点分别进行了具体讨论，并对模型内的速度和概率分布情况进行了分析，来验证模型的可行性。1三维随机介质的建模原理及步骤1.1随机介质的建模原理随机介质模型可由大尺度和小尺度的非均匀的两部分所组成，大尺度部分描绘介质的平均特性，即传统意义上的地质模型；小尺度部分是加在地质模型上的随机扰动，通常使用一个均值为零的、平稳的空间随机经

4、过来表示。以三维弹性随机介质为例。各向同性的弹性介质由其密度和拉梅系数、所确定。1.2随机介质的构建步骤构建平稳随机经过的步骤如下：1选择自相关函数。如今的自相关函数有高斯型、指数型、VonKarman型、混合型及一些改良型的自相关函数，它们分别具有不同的特点，适应不同的地质情况。本文采用VonKarman型相关函数，它描绘的随机介质具有多尺度、自类似的特点。2随机介质模型的误差压制由于在离散数据的计算经过中会产生误差，使得随机函数不知足假设条件，因而在计算随机功率谱时，引入三维锥形函数来消除这种误差。建立了大小为250250250个网格步长为1m的三维锥形函数见图1，锥形函数的半径为0.5m

5、-1，从图1中能够看出，锥形函数的作用是压制随机模型中的低波数分量。3三维随机介质的建立及特点分析为了便于比照分析，在引入三维锥形函数的前提下，假设各模型网格大小250250250步长为1m，背景速度v0=3500m/s，模型的其他特点通过改变标准差、谱指数、自相关长度来控制。3.1固定标准差、谱指数，研究自相关长度对模型的影响固定模型参量：标准差01，谱指数长度axayaz=2，改变自相关长度a、b、c和能够产生不同的三维随机介质模型，详情见图2。图2a和2b的自相关长度均相等，模拟的是各向同性时的三维随机介质模型。由图2a和2b能够看出，自相关长度越大，模拟的三维随机介质模型的尺度越大。图

6、2c和2d的自相关长度不同，增大某一方向的自相关长度，则该方向速度趋向于一致，随机模型的层状特点逐步明显，通过增加1个或2个方向的自相关长度到无穷大能够模拟各种各向均匀的薄互层状随机介质18。3.2固定标准差、自相关长度，研究谱指数对模型的影响固定模型参量：标准差01，自相关长度a=b=c=1m，改变谱指数ax，ay，az，能够产生不同的三维随机介质模型，详情见图3。图3a和3b同样模拟的是各向同性时的三维随机介质模型，由图3a和3b能够看出，谱指数的大小对应着随机介质尺度的大小，并且谱指数越大，所得到的随机介质模型越光滑。图3c和3d表明，通过增加1个或2个方向的谱指数到无穷大同样能够模拟各

7、种各向均匀的薄互层状随机介质。3.3固定谱指数、自相关长度，研究标准差对模型的影响固定模型参量：自相关长度a=b=c=1m，谱指数ax=ay=az=2，引入三维锥形函数，改变标准差，能够产生不同的三维随机介质模型，详情见图4。从图4中能够看出，其他参数不变的情况下，改变标准差，从所得的三维随机介质模型几乎看不出区别，但是从旁边的色标值能够看出，标准差越小，其速度扰动越小，即标准差不同，得到的随机扰动的范围就不同，进而能够反映出实际介质的非均质性。4随机介质模型中的速度变化及概率分布从上面所建立的模型图中能够非常明晰地看到随机介质模型的特点，改变自相关长度、谱指数和标准差，能够得到不同的介质模型

8、，但是无法了解模型内的速度变化和概率分布情况。从图2a中所建立的模型的x、y、z方向上随机抽取一维数据体，详情见图5。另外统计出各个速度概率分布，速度变化区间为500m/s，详情见图6。从图5和图6中能够得到，模型内的速度值在x、y、z方向上都是在背景速度v0=3500m/s上下浮动，与随机介质的基本假设比拟吻合，即随机扰动的均值为零，但扰动的范围与标准差有关。从图6中能够看出，模型内的速度在x、y、z方向上呈正态分布，与相差越小，速度分布样本数越多，概率越大。统计结果表明，随机介质模型内的速度大多分布在1倍标准差范围内，其概率大约为75%，2倍标准差范围内的速度概率大约为15%。因而可知，随

9、机介质模型内的速度符合统计规律，讲明所建立的模型具有可行性。5结论在前人的基础上，给出了三维VonKarman型自相关函数模拟随机介质的具体流程。在建模经过中，通过引入三维锥形函数，减小了随机建模时离散数据的计算误差。为了讨论自相关长度、谱指数等参数对模拟随机介质时的影响，通过改变这些参数来模拟三维随机介质模型，可知每种参数对模型结果的影响不同。从模拟的图能够看出，运用VonKarman型自相关函数模拟的三维随机介质模型与实际的地下岩芯图特别类似，通过改变各种参数(自相关长度、谱指数、标准差等)，能够灵敏地描绘多种实际介质在小尺度上的非均匀性的特点。随机介质模型理论为进一步研究实际三维地质体中

10、小尺度的非均质性对地震波传播特征的影响打下了一定的基础。参考文献：1李灿苹，王南萍，李志宏.VonKarman型自相关函数模拟随机介质J.物探与化探，2010，341：98-102.2李灿苹，刘学伟，李敏锋，等.非均匀地质体散射波特征初探J.石油物探，2006，451：134-140.3何又雄，姚姚.分形方法在模拟非均匀介质模型中的应用J.石油物探，2005，443：199-201.4何又雄，姚姚.各向异性分形介质模拟J.江汉石油学院学报，2004，264：67-69.5IkelleL.，YungS.，DaubeF.2-Drandommediawithellipsoidalautocorrel

11、ationfunctionsJ.Geophysics，1993，589：1359-1372.6KornM.，SeismicwaveinrandommediaM1993：247-269.7EngintavS.，CanitezN.，Modelingofmulti-scalemediaindiscreteformM1997:77-96.8奚先，姚姚.二维弹性随机介质中的波场特征J.石油地球物理勘探，2004，201：679-685.9奚先，姚姚.二维随机介质及波动方程正演模拟J.石油地球物理勘探，2001，365：546-552.10奚先，姚姚.非平稳随机介质模型J.石油地球物理勘探，2005，40

12、1：71-75.11姚姚，奚先.区域多尺度随机介质模型及其波场分析J.石油物探，2004，431：1-7.12奚先，姚姚.随机介质模型的模拟与混合型随机介质J.地球科学，2002，271：67-71.13吴何珍，符力耘，兰晓雯.基于随机介质模型的储层非均质性分析J.地球物理学进展，2008，233：793-799.14陈可洋.三维随机建模方法及其波场模拟分析J.勘探地球物理进展，2009，325：315-320.15王金山，陈可洋，吴清岭，等.随机介质模型的一种构造方法J.物探与化探，2010，342：191-194.16刘春园，朱生旺，魏修成，等.随机介质地震波正演模拟在碳酸盐岩储层预测中的

13、应用J.石油物探，2010，492：133-139.17殷学鑫，刘洋.二维随机介质模型正演模拟及其波场分析J.石油地球物理勘探，2011，466：862-872.18郭乃川，王尚旭，郭锐，等.地震勘探中三维小尺度非均匀性随机介质模型的建立及其特点分析J.石油天然气学报，2012，347：62-67.19徐振旺，安勇，支玲，等.运用分形方法模拟地震勘探中的三维随机介质J.地球物理学进展，2014，293：1319-1325.20BrichF.Thevelocityofcompressionalwavesinrocksto10kilobars-Part2J.JournalofGeophysicalResearch，1961，667：21992224.21MarpleS.L.，DigitalspectralanalysiswithapplicationM.NewJersey，1987