基于地震波形指示的碳酸盐岩储集体反演方法——以塔河油田中-下奥陶统为例-韩长城.pdf

《基于地震波形指示的碳酸盐岩储集体反演方法——以塔河油田中-下奥陶统为例-韩长城.pdf》由会员分享，可在线阅读，更多相关《基于地震波形指示的碳酸盐岩储集体反演方法——以塔河油田中-下奥陶统为例-韩长城.pdf（9页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、磊沾-5天然气必彦第38卷第4期 OIL&GAS GEOLOGY 2017年8月文章编号：02539985(2017)04082209 doi：101 1743o9920170419基于地震波形指示的碳酸盐岩储集体反演方法以塔河油田中一下奥陶统为例韩长城1，林承焰12，任丽华12，董春梅12，魏婷41中国石油大学(华东)地球科学与技术学院，山东青岛266580；2山东省油藏地质重点实验室，山东青岛266580；3中国石油油气储层重点实验室中国石油大学(华东)研究室，山东青岛266580：4中国石油新疆油田分公司石西采油厂，新疆克拉玛依834000摘要：塔河油田奥陶系碳酸盐岩岩溶储集体埋深大，分

2、布不均，类型复杂，具有较强的非均质性，储集体类型识别一直是该区油气勘探和开发的难题。为此，基于贝叶斯判别理论和马尔科夫链蒙特卡洛抽样算法，实现了一种新型波阻抗反演方法。在储集体类型识别、地震反射特征和岩石物理分析基础上，利用地震波形相似性优选样本井，参照样本空间分布距离和曲线分布特征建立初始阻抗模型。通过建立建议分布与先验信息之间的关系，实现先验解空间的有效采样，同时采用MetropolisHastings抽样算法对贝叶斯后验概率分布抽样，从而获得最大后验概率解。将该反演方法应用到塔河油田奥陶系岩溶储集体类型识别中，在反演精度得到明显提高的同时，充分利用了地震波形横向变化特征反映岩溶储集体空间

3、分布位置和横向展布规律，有效识别了岩溶储集体类型，为非均质极强的岩溶缝洞储集体的精细表征提供有力支撑。关键词：波形指示反演；贝叶斯理论；岩溶储集体；奥陶系；塔河油田；塔里木盆地中图分类号：TEl222 文献标识码：AWaveformindicationbased seismic inversion of carbonate reservoirs：A case study of the Lower-Middle Ordovician in Tahe oilfield，Tarim BasinHan Changchen91，Lin Chengyanl，2一，Ren Lihual2”，Dong Chu

4、nmeil，2，Wei Tin941School of Geosciences，China University ofPetroleum，Qingdao，Shandong 266580，China；2Reservoir Geofogy Key Laboratory ofShandongProvince，China University ofPetroleum(East China)，Qingdao，Shandong 266580，China；3PetroChina Key Laboratory ofOilGasReservoir，China University of Petroleum(Ea

5、st China)，Qingdao，Shandong 266580，China；4Shixi Oil Production Plant，Xinjiang Oilfield Company，PetroChina，Karamay，Xinjiang 834000，ChinaAbstract：Carbonate karst reservoirs of the LowerMiddle Ordovician in Tahe oilfield feature in large burial depth，irregular distribution，complex reservoir types and st

6、rong heterogeneity，thus reservoir type identification is always a challengingduring oil and gas exploration and development in the areaThis paper documented a new impedance inversion method，iewaveformindicationbased inversionbased on Bayesian discriminant theory and Markov Chain Monte Carlo sampling

7、algorithmThe method，based on reservoir type identification，seismic reflection characteristics and petrophysical analyses，esn be used to optimize sample wells by observing seismic waveform similarities，and establish an initial impedancemodel by referring to sampling spacing and curve characteristicsT

8、he relationship between proposed distribution and priorinformation is firstly established to achieve efficient sampling of a priori solution spaceA MetropolisHastings sampling algorithm is then used to sample the posterior probability distribution，SO as to obtain the maximum posterior probability so

9、lutionThe application of the method in identifying the types of the Ordovician karst reservoirs in Tahe Oilfield shows signifi-cant improvement of inversion accuracy and a better utilization of horizontal variation of seismic waveforms to reveal effectivelv the types and spatial and lateral distribu

10、tion of the reservoirs，and provides fine characterization of the reservoirsKey words：waveformindicationbased inversion，Bayesian theory，carbonate reservoir，Ordovician，Tahe oilfield，Tarim Basin收稿日期：2017030l；修订日期：20170630。第一作者简介：韩长城(1984)，男，博士研究生，碳酸盐岩油藏精细描述。E-mail：517892849qq通讯作者简介：任丽华(1979一)，女，副教授，碳酸盐

11、岩油藏精细描述。E-mail：rlh-79163con。基金项目：国家科技重大专项(2017ZX05009001)。万方数据第4期 韩长城，等基于地震波形指示的碳酸盐岩储集体反演方法塔里木盆地奥陶系碳酸盐岩岩溶储集体分布不均，储集体类型复杂，具有较强的非均质性Ll引。碳酸盐岩基质具有低孔、低渗的特点，储集空间主要为后期溶蚀作用形成的洞穴、孔洞和裂缝，根据储集体的形态和组合方式不同，将其分为溶蚀洞穴型、溶蚀孑L缝型和裂缝型3类H J。缝洞储集体纵、横向变化快，规模大小不一，分布规律复杂，主要依靠地球物理资料进行缝洞储集体识别预测。但由于奥陶系埋深大，一般大于5500 in，碳酸盐岩致密，地震波传

12、播速度快。5j，反射信号弱，地球物理特性非常复杂，给地震勘探开发带来了很多困难。近年来，随着地球物理技术的发展，很多学者尝试利用多种地球物理方法地对碳酸盐岩缝洞储集体进行了预测。诸如，正演模拟缝洞储集体反射特征，建立缝同储集体形态及岩石物理参数与地震响应之间的关系7 o；利用地震属性预测缝洞体引，如振幅变化率属性、相干属性集多属性优化技术等均取得了较好效果；通过地震波形特征建立储集体和地震之间的响应关系。10。111；利用波阻抗反演。1纠和地质统计学反演方法对缝洞储集体进行预测孓14 J。这些研究成果均较好指导碳酸盐岩缝洞储集体勘探，但仍然存在较多问题，如地震属性多解性较强；波阻抗反演分辨率太

13、低，无法满足小型缝洞体预测需求；测井约束反演高频分量主要来源于初始阻抗模型，多解性强；地质统计学反演随机性较强，对井位分布要求较高。地震波形指示反演是贝叶斯判别理论n5。171和马尔科夫链蒙特卡洛8191抽样算法实现的一种反演方法。在贝叶斯理论框架下，通过给定地质认识、测井解释等建立先验解空间，建立建议分布与先验解空间之间的关系，实现先验解空间的有效采样，采用MetropolisHastings抽样算法对后验概率分布抽样，选取目标函数最大值时的解作为可行随机实现，求取多次可行实现均值作为期望输出。该反演结果在纵向上分辨率得到提高，同时在横向上对储集体边界精细刻画。1反演方法原理11 目标函数建

14、立根据贝叶斯条件概率定义，事件的后验概率公式为：如2帮2嚣嚣恭oc P(m)p(n m) (1)式中：m模型参数，对应的先验分布为P(m)，先验分布P(m)与观察数据z无关，主要由测井资料得到；n与待求模型参数m有关的观测样本；P(n，m)m和n的联合分布；P(n)n发生的概率；P(nIm)似然函数，表示在模型参数m下的条件概率；P(m n)在观察数据n下模型的后验概率。一般由声波曲线所计算的波阻抗曲线是由两部分组成：Y=z+ (2)式中：)，由声波计算得到的波阻抗值，gcm一ms；z待求解的地下各地层的实际波阻抗值；N背景噪音，一般情况下，假设背景噪音N服从高斯分布心0I，其均值为0。由于(

15、zz。)2与一样均符合高斯分布，则目标函数可以表示为：Jl(z)：业i出 (3)式中：z，波阻抗z的期望值，邻主要通过选取的样本井的波阻抗曲线求得；盯均方差。波阻抗反演过程就是通过公式(3)研究z为何值时从而使，。达到最小。波阻抗反演存在不稳定性和多解性，使用一种方法很难获得较理想解；同时由于地震资料带限宽度的限制，只能提供波阻抗反演的中频成分，而高频和低频成分则主要通过测井数据、地质信息获得心1I，因此可以用贝叶斯推理把反演问题与观测数据信息、模型信息以及先验信息联系起来，引入先验信息对其进行约束，从而获得较稳定的反演解心2。在贝叶斯理论下，具有先验信息条件时，经过贝叶斯推理，得到待反演模型

16、参数的目标函数为：J(z)=J。(z)+A 2厶(z) (4)式中：z波阻抗，gcm一ms。厶(z)势函数的和，是与先验信息如测井资料、地质信息等有关的函数。J。(三)与某些后验信息有关的函数。常数A为平滑参数，用于协调，。(彳)和，：(z)之间相互影响。公式(4)进一步改写为：J(z)：业暑止+粕乜)万方数据824 石油与天然气地质 第38卷+A2咖(D：(z)届) (5)C式中：西势函数；A正规参数，一种平滑参数；6正规参数，一种刻度参数，调节不连续处调节梯度值。先验项西。(z)中的k可取1，2，3。当k=1时表示先验项求和是在目标点最邻近的水平和垂直方向所有区域进行。k=2，3时表示在平

17、面或二次曲面上对邻近点之间的梯度矢量和Hsssian矩阵求和。2川。公式(4)由两部分组成，第一部分是关于测量模型和数据之间的某种一致性，当其值达到最小时，一致性会变得很好，但解稳定性变差；第二部分是关于解“光滑性”的测量标准，当其值取最小时，获得的近似解较稳定。12 MCMC方法和Metropolis-Hastings抽样准则MCMC方法是在贝叶斯框架下，用已有资料进行约束，既可使最优解满足参数的统计特性，同时融入了先验信息，解的精度得到提高旧3I。由于反演问题的非线性性，后验分布常常没有显式的表达式，因此，引入MCMC方法对后验概率分布进行采样，从而获得参数的估计值。其基本思想是：通过构建

18、一个平稳分布为7r(x)的马尔科夫链来获得7r(x)的样本，被估计的模型参数值是链的状态空间，被估计的后验分布为链的贝叶斯后验分布。利用马尔科夫链蒙特卡罗(MCMC)方法对贝叶斯理论下的后验概率分布进行随机抽样，反演问题希望得到的是所估参数的后验概率分布，因此各个参数的马尔科夫链应收敛于所估计参数的后验概率分布旧4I。设M为某一空间，n为产生的总样本数，m为马尔科夫链达到平稳时的样本数，则MCMC的思路为2 5I：(1)构建一条Markov链，使其收敛至平稳分布仃(石)，常用的构造转移核的是Gibbs抽样和MetropolisHastings抽样；(2)产生样本，由M中的某一点戈出发，用(1)

19、中的Markov链进行抽样模拟，产生点序列咒l，一，戈。；(3)蒙特卡罗积分。任一函数八石)的期望估计为1“科八戈)=二：八石)(6)|Lnt；i+I式中：n总样本数；m马尔科夫链达到平稳时的样本数。MetropolisHastings抽样的基本思路。1引是：选择一转移函数q(x；xi。1)和初始值算叭，若第i次迭代开始时的参数值为戈“。1，则第i次迭代过程为：(1)从q(x；戈“一)中抽取一个备选值戈；(2)计算接受概率刊n君篆b)(3)以概率Ot(戈。1，z 7)，置戈)-戈7，以概率1一d(戈H，戈)，置戈i=戈H；(4)重复(1)一(3)步t次，则可以得到后验样本戈，石扪，戈扪，石。从

20、而可计算后验分布的各阶矩以进行统计推理。三维地震资料是分布密集的空间结构化数据，地震波形的横向变化反映不同沉积环境和岩性组合在空间上的变化规律。因此地震波形指示反演是在沉积学基本原理的基础上，利用地震波形相似性优选相关井样本，参照样本空间分布距离和曲线分布特征建立初始阻抗模型。同时通过对同一相带内声波阻抗曲线的分析，提出了“相截频率”，进一步拓宽频带分布，统计样本井纵波阻抗建立先验概率函数。将初始模型与地震波阻抗体进行匹配滤波，求得似然函数，对样本进行多尺度分解，逐步滤除高频成分；最后基于贝叶斯理论，联合似然分布与先验分布得到后验概率分布，并将其作为目标函数；基于波形指示优选的样本在空间上具有

21、较好的相关性，采用Metropolis-Hastings抽样算法对后验概率分布抽样，选取目标函数最大值时的解作为可行随机实现，求取多次可行实现的均值作为期望输出。因此该反演方法在高频成分确定过程中充分利用了地震波形的横向变化特征，在提高垂向分辨率的同时，随机性减小。2应用实例21地质概况塔河油田位于塔里木盆地北部沙雅隆起中段的阿克库勒凸起西南端。2 6I，先后经历了加里东、海西、印支一燕山和喜马拉雅期构造升降运动，遭受多期剥蚀及岩溶作用。塔里木盆地奥陶系主要为一套台地相碳酸盐岩沉积，发育有下奥陶统蓬莱坝组、中一下奥陶统鹰山组、中奥陶统一间房组，上奥陶统切尔巴克组、良里塔格组、桑塔木组，主要目的

22、层为中、下奥陶统的鹰山组和一间房组。鹰山组主要岩性为褐灰色泥微晶灰岩、白云质灰岩和灰质白云岩；一间房组主要岩性为灰色砂屑灰岩、生物灰岩和泥微晶灰岩等。加里东期和万方数据第4期 韩长城，等基于地震波形指示的碳酸盐岩储集体反演方法海西期发生两期重要的构造运动，加里东中期第1幕构造运动形成了一间房组顶面平行不整合面，主要发生大气淡水溶蚀作用；加里东中期幕构造运动，在良里塔格组与上覆桑塔木组形成角度不整合面，遭受强烈剥蚀；晚奥陶世末期，发生加里东中期幕，构造抬升，中、上奥陶统遭受剥蚀；海西早期构造运动幅度较大，抬升剥蚀较剧烈，形成大型不整合面，岩溶作用较强烈(图1)。22储集体类型及特征通过对研究区1

23、3口井的岩心、薄片与录井等资料的观察分析，塔河油田奥陶系岩溶储层储集空间类型多样，主要包括孔隙、裂缝和大型溶洞。根据它们的组合关系不同，储集体类型主要有溶蚀洞穴型、裂缝一孑L洞型和裂缝型3类4。溶蚀洞穴型储集体的储集空间多为未充填或半充填的大型洞穴系统。岩溶垮塌充填的沉积物中可见分选较差，棱角一次棱角状的灰岩角砾(图2a)。钻井过程中发生井漏和放空漏失，声波时差呈现明显“U”型高值，密度显著降低，中子孑L隙度增大，电阻率曲线明显降低，深浅电阻率出现正幅值差，在泥质充填处自然伽马值升高，在地震剖面上表现为强的“串珠”状反射。裂缝一孑L洞型储集体储集空间多为裂缝和溶蚀孔洞，孔洞多为沿先期裂缝、微裂

24、隙、缝合线等局部溶蚀扩大形成，裂缝和溶孔对油气的储集、渗流都有重要的作用，后期多被方解石或有机质全充填一半充填(图2b，c)。常规测井曲线特征为：密度稍微降低，声波时差发生轻微跳变，电阻率值降低，在地震剖面上多表现为杂乱反射特征。裂缝型储集体的储集空间是构造缝、溶蚀缝和成岩缝等。裂缝的成因主要受构造作用和成岩作用控制，经历多期次构造挤压作用，可见多组不同角度裂缝相互交叉形成裂缝网络。具有多期次发育特点，后期被方解石、有机质等充填一半充填，部分直立缝未充填，裂缝面存在油质(图2d，e)。23地震反射特征分析在正演模拟和前人研究的基础上。27之引，结合塔河油田碳酸盐岩储层实际地质情况及地震反射特征

25、，将碳酸盐岩岩溶储集体的地震反射结构分为“串珠”状反射、杂乱反射和弱反射3类。“串珠”状反射是塔里木盆地碳酸盐岩勘探和开发的主要研究对象，且已证实具有较好产能(图3a)。当碳酸盐岩缝洞与不整合面间距大于30 m时，溶洞与致密基岩之间会形成多次绕射，绕射波经过叠前偏移成像可在垂向上形成多个短轴强反射悼4。通常是垂直洞穴、大型椭圆型储集体及多层储集体的响应特征。当风化壳表层储集体不发育，内幕储集体发育时，表现为表层强反射、内幕“串珠”状反射(图3b)；当风化壳表层和内幕区储集体均发育时，且之间有较厚致密层隔挡时，表现为表层弱反射、内幕“串珠”状反射(图3c)。杂乱反射常由一些连续性较差、方向不定、

26、形状各异、振幅有强有弱的短同向轴组成。当纵向厚度较小的孔洞和裂缝组成的横向随机分布孔洞储集体时，地震反射结构特征表现为强杂乱反射(图3d，e)。圃地名 口断层曰盆地边界习二级构造单元边界三三级构造单元边界研究区匡虽泥岩臣丑灰岩巨丑颗粒灰岩匡虿自云岩图1 塔河油田构造位置和地层综合柱状冈Fig1 Stiu(ture units and stlatigraphic column in Tahe oilfie万方数据826尘0娑1 50AC(gsm瑟旧 l1。1) |岬l蜃岩性削面RO(nIn)20 00(RS(f2m120 00岩心及微观照片r1 r一1_i：磕| 基釜|l匿虿_ l洞穴型 裂缝一

27、孔洞型 裂缝型 放空漏失 灰色泥品灰岩 褐色泥隔灰岩图2塔河油田Tk210井地层综合柱状图Fig2 Composite columnar section of Well Tk210 in Tahe oilfield弱反射在地震剖面上表现为振幅值小或空白的反射值。内幕区发育规模较小的孑L洞型储集体、裂缝型储集体或缝洞储集体距风化壳距离较近时表现为弱反射特征(图3f)。24储层岩石物理分析首先对研究区测井资料进行标准化预处理，主要针对的是因测井仪器误差等原因导致测井数据出现的偏差。然后对测井资料进行统计分析，碳酸盐岩基岩致密，地震波传播速度和密度均较高，地震波的平均速度达到6 000 ms。地层中

28、发育裂缝、溶孑L和洞穴时导致声波时差明显增大、密度降低，即波阻抗值显著降低，同时深侧向电阻率也明显降低。因此波阻抗与储集体之间存在明显的对应关系。图4给出了不同储集体类型纵波阻抗与深侧向电阻率交汇分析结果。溶蚀洞穴型储集体波阻抗值最低，从12 000第38卷15 500 gemms，裂缝一孔洞型储集体波阻抗值域为14 10016 400 gemms，裂缝型储集体波阻抗值较高，从15 70017 500 gemms，而碳酸盐岩基岩波阻抗值多为大于17 000 gemms。因此波阻抗值能较好区分各类型储集体。25 子波提取及频率参数设置子波的提取主要采用循环迭代法。其提取流程为：利用地震道的自相关

29、性计算初始子波的振幅谱和相位谱；利用振幅谱和相位谱的信息合成一个理论雷克子波，用理论子波合成的记录进行初步标定，提取井旁地震道实际地震子波，用该子波制作合成记录并重新调整时深关系。如此反复，直至得到相位振幅谱变化稳定的精细子波。地震数据的频带范围影响反演结果的分辨率。地震资料主频为30 Hz，有效频带范围为1070 Hz，主要集中在中频范围；而反演中的低频信似锄眦似觚叭一蓦墨一万方数据韩长城，等基于地震波形指示的碳酸盐岩储集体反演方法羹I二一II 喈：喜，。 丽， 喜，。量 lI 吾、 T、 l3 700一固 l 冈”00斤L1F 眇 IT1 bk260一羹西暑念rllr 1师一I t移 一一

30、一一、斤1一一一_1 500-匣_l黧一I裂缝j图4塔河油田奥陶系储集体纵波阻抗与深侧向电阻率测井值交汇图Fig4 Crossplot of Pwave impedance VSRD of the0rdovician reservoirs in Tahe oilfield息和高频信息主要由测井曲线提供。因此在地震波形指示反演中，低频成分主要从测井资料中获取，中频成分通过地震资料求取，高频成分是在地震波形指示下进行的随机模拟结果，整个过程是由确定到随机的逐步过渡，随机成分较少。本次反演的低频参数和高频参数优选如下：低频参数中高通频率取6 Hz，高截频率取10 Hz；高频参数中低截频率取60 Hz

31、，低通频率取80 Hz，高通频率取200 Hz，高截频率取250 Hz。26反演效果分析通过对塔河油田的地质、地震以及测井等资料分析，确定了碳酸盐岩储集体类型及对应的地震反射波特征，得到了岩石物理参数模型。在贝叶斯理论基础上，以地质认识、测井数据和地震解释层位作为先验信息来约束反演结果，综合其概率密度函数得到岩溶储集体发育情况的后验概率分布函数，即缝洞储集体的空间发育规律。图5a是原始地震剖面、测井约束反演和波形指示反演结果。在原始剖面上，缝洞储集体为大片的“串珠”状强震幅反射以及弱振幅反射，能大概确定储集体分布的位置，但对纵向上储集体的叠置关系、横向上储集体的边界范围以及储集体的具体形状无法

32、确定。测井约束地震反演是基于模型的反演，从图5b中可以看出，在井点处反演结果垂向分辨率较高，但由于其初始波阻抗模型是由井资料内插和外推实现的，其横向连续性较好。但塔河油田奥陶系碳酸盐岩缝洞储OOOOOOOOO0OOOOO8406284，27一llg五已、Q世万方数据石油与天然气地质 第38卷集体埋藏深度大，储集体分布受断裂、岩溶作用等多因素控制，储集体形态复杂，纵横向分布变化快，非均质性极强，因此对于碳酸盐岩缝洞型储集体这类横向变化较快的储层，测井约束地震反演井间预测较好较差。地震波形的横向变化可以较好的反映缝洞储集体横向变化规律，不同地震反射特征所反映的储集体的类型不同，如大型溶洞储集体多为

33、“串珠”状反射，而杂乱反射多为裂缝一孔洞型储集体，弱反射多为裂缝发育带。因此在反演过程中基于地震波形特征优选样本，利用样本井的原始数据和空间结构特点，对未知样点进行线性无偏、最优估计。在地震波形约束条件下，有效的提高高频信息，反演结果能较好地区分各类储集体的空间分布规律(图5c)。波阻抗值低的部分(即波阻抗值小于15 500 gcmms。1部分)为规模较大、物性较好的溶洞型储集体。对于非溶洞型储集体，波阻抗变化主要反映在物性差异上，波阻抗值较低处为裂缝一孔洞型储集体(即波阻抗值域为14 10016 400gcmm-s。1)，如Tk211井测井解释在表层发育一段大约30 m的裂缝一孔洞型储集体，

34、在地震剖面上无法看出，而反演剖面上有较好的反应。裂缝型储集体主Tk211 Tk445ib形； ：譬k一 一：。 二弋J- ，广懈蕊、f、：静J孓l厂一 一、夕 一 r-沙、图5 塔河油田过Tk21 1一Tk4J45井地震剖面波和反演剖面Fig5 Seismic and inversion sections crossing Tk21 1 and Tk445 wells in Tahe oilfieha原始地震剖面：b测井约束波阻抗反演剖面；(波形指示反演剖面日OOO剖似似乇一一一O2468彬m如如如如如j，44667788波：I万方数据第4期 帏成，帑肇J-J断史形指示的碳酸盐岩储集体反演方法

35、 829l抗n1S000500图6塔河油田波形指示反演结果平面属性Fig6 Map showing wavef01rnindication-basedinversion results in Tahe oilflled要是波阻抗值较高部分，主要发育在风化壳表层和构造变形部位。波阻抗值最高处为碳酸盐岩基岩部分，在反演剖面上主要为蓝色一浅蓝色部分。在反演平面图上也可以看出缝洞储集体多为管道状暗河和孤立型溶洞，在溶洞周围有孑L洞一裂缝型储集体发育，储集体主要分布在研究区的北部(图6)。缝洞储集体横向分布较自然，与地质认识相符。3 结论1)波形指示反演是在贝叶斯框架下，通过建立建议分布与先验解空间之间

36、的关系，实现先验解空间的有效采样，然后采用MetropolisHastings抽样算法对后验概率分布抽样，从而实现的新型波阻抗反演方法。该过程中地质认识及测井资料起着软约束作用，在尊重地震本身趋势的情况下，根据波形特征储集体纵向分辨能力得到提高的波阻抗数据体，更加突出了地震资料在井间储集体预测作用。2)缝洞碳酸盐岩储集体非均质性强，分布规律复杂，储集体预测识别一直是一个难题。溶蚀洞穴型储集体规模较大，物性好，利用该方法识别效果较好，裂缝一TLN型储集体多是成带分布，也能较好的识别；裂缝型储集体与基岩的波阻抗值差异较小，预测结果准确性仍需要进一步确定。同时由于引人的高频信息是在波形指示下进行的随

37、机模拟，因此，需要注意样本井的选取和高频信息的选取。参考文献1 韩长城，林承焰，鲁新便，等塔河油田奥陶系碳酸盐岩岩溶斜坡断控岩溶储层特征及形成机制J石油与天然气地质，2016，37(5)：644652Han Changcheng，Lin Chengyan，Lu Xinbian，et a1Characterizationand genesis of faultco,m-oiled karst reservoirs in Ordovieian carbonate kalst shoe of TaheOilfieldTarim BasinJOilGas Geology2016，37(5)：64465

38、22 刘立峰，孙赞东，杨海军，等缝洞型碳酸盐岩储层地震属性优化方法及应用J石油地球物理勘探，2009，44(6)：747754，Liu Lifeng，Sun Zandong，Yang Haijun，et a1Seismic attribute opti-mization method and its application for fracturedvug对carbonatereservoirJOil Geohysical Prospecting，2009，44(6)：7477543 张学丰，李明，陈志勇，等塔北哈托哈塘奥陶系碳酸盐岩岩溶储层发育特征及主要岩溶期次J岩石学报，2012，28(3

39、)：815826Zhang Xuefeng，Li Ming，Chen Zhiyong，et a1Characteristics andkarstifieation of the 0rdovician carbonate reservoirHalahatang area，north TarimBasinJActa Petrologica Sini(-a，2012，28(3)：8】58264 韩长城，林承焰，任丽华，等塔里木盆地塔河lO区奥陶系断裂特征及对岩溶储层的控制作用J天然气地球科学，2016，27(5)：790798Han Changeheng，1in Chengyan，Ren Lihua

40、，et a1Characteristics ofOrdovician fault in the Block 10 of Tahe Oilfield，Tarim Basin and itscontrolling effect on karst reservoirsJNatural Gas Geoscience，201627(5)：7907985赵裕辉，胡建中，鲁新便等碳酸盐岩缝洞型储集体识别与体积估算J石油地球物理勘探，2010，45(5)：720724Zhao Yuhui，Hu Jianzhong，Iu Xinbianet a1Identification of carbonate fract

41、uredvuggy reservoir and volumeestimationJOil Geohysical Prospecting，2010，45(5)：7207246 马灵伟，顾汉明李宗杰，等正演模拟碳酸盐岩缝洞型储层反射特征J石油地球物理勘探，2015，50(2)：290297Ma Lingwei，GuHanming，Li Zongjie，et a1Simulation of carbonatefracturecaverll reservoir retie(-tion characteristics with forward modeli“gJOil Geohysical Prospe

42、cting，2015，50(2)：2902977 李倩，狄帮让，魏建新碳酸盐岩储层孔洞体积的地震物理模拟估算J石油地球物理勘探，2014，49(6)：11471156Li QianDi Bangrangwei JianxinCarbonate reservoir cave volumeestimation with seismic physical simulationJOil Geohysical Prospeeting，2014，49(6)：1147一11568 韩革华，漆立新，李宗杰，等塔|giflt田奥陶系碳酸盐岩缝洞型储层预测技术J石油与天然气地质，2006，27(6)：860870

43、Han Gehua，Qi lixin，Li Zongjie，et a1Prediction of the Ordovicianfractured-vuggy carbonate reservoirs in TaheoilfieldJOil&GasGeoloD,，2006，27(6)：8608709王光付碳酸盐岩溶洞型储层综合识别及预测方法J石油学报，2008，29(1)：4751Wang GuangfuIntegrative identification and prediction methods forcarbonate rock cave reservoirJActa Petrolei

44、Sinica，2008，29(1)：475110温志新，王红漫，漆立新，等塔河油田奥陶系缝洞型碳酸盐岩储层预测研究J地学前缘，2008，15(1)：94100Wen Zhixin，Igang HongmanQi Lixin，et a1Research on the prediction of 0rdovician carbonatereservoir of fissurecave type in the TaheOilfieldJjEarth Science Frontiers，2008，15(1)：9410011罗亚能，黄捍东，王玉梅，等地震波形反演与测井联合的i维建模方法!J石油地球物理勘

45、探2016，5I(5)：947954iiJ5nn71S8rI万方数据830 石油与天然气地质 第38卷一-一LuoYaneng，Huang Handong，Wang Yumei，et a1Three-dimensional and Geochemical Exploration，2015，39(6)：12111216modeling based on integration of full wave inversion and well logging 20 Rimstad K，More HImpact of rockphysics depth trends and MarkovdataJ0i

46、l Geophysical Prospecting，2016，51(5)：947954 random fields on hierarchical Bayesian lithologyfluid predictionJ12孙素琴，杨玉杰塔里木盆地顺西区块地震反演及储层预测J Geophysics，2010，75(4)：R93一R108石油物探，2012，51(3)：304311 21 张宏兵，杨长春正则参数控制下的波阻抗约束反演J地球物Sun Suqin，Yang YujieSeismic inversion and reservoir prediction of 理学报，2003，46(6)

47、：827834Shunxi block in TarimBasinJGeophysical Prospecting for Petroleum， Zhang Hongbing，YANG ChangehunA constrained impedance inver-2012，51(3)：304311 sion method controlled by regularizedJ1Chinese Journal of Geo13孙炜，王彦春，李玉凤，等基于地质统计学反演的碳酸盐岩孔缝physics，2003。46(6)：827834洞预测研究J现代地质，2012，26(6)：12581264 22宋维琪，朱海伟，姜宇东，等地面微地震资料震源定位的贝叶斯Sun Wei，Wang Yanchun，Li Yufeng，et a1Fractures and solu