基于adm谱反演的高分辨率裂缝预测技术研究及应用-张华.pdf

《基于adm谱反演的高分辨率裂缝预测技术研究及应用-张华.pdf》由会员分享，可在线阅读，更多相关《基于adm谱反演的高分辨率裂缝预测技术研究及应用-张华.pdf（10页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、第55卷第5期2016年9月石油物探GEOPHYSICAL PROSPECTING FOR PETROLEUMV0155，No5Sep，2016张华，贺振华，李亚林，等基于ADM谱反演的高分辨率裂缝预测技术研究及应用J石油物探，2016，55(5)：737745ZHANG Hua，HE Zhenhua，LI Yalin，et a1Research arld application of high resolution fracture prediction technology basedon趟)M spectral inversion口Gphysical Prospecting for Pe

2、tmleum，2016，55(5)：737745基于ADM谱反演的高分辨率裂缝预测技术研究及应用张 华123”，贺振华12，李亚林3”，张晓斌3”，何光明3”，李德珍3“(1“油气藏地质与开发工程”国家重点实验室，四川成都610059；2成都理工大学地球物理学院，四川成都610059；3中国石油天然气集团公司川庆钻探工程有限公司地球物理勘探公司，四川成都610213；4中国石油天然气集团公司山地地震技术试验基地，四川成都610213)摘要：从谱反演基本原理出发，推导出多层稀疏反射系数目标函数，利用基追踪交替方向法(Altemating Direction Method，ADM)谱反演算法求取高

3、精度反射系数，与谱模拟解析法提取的宽频子波褶积重构，得到保持低频信息的高分辨率剖面，利用改进型特征值相干体算法对此剖面进行裂缝预测。理论测试和实际资料应用结果表明，ADM谱反演高分辨率技术能得到极性、位置准确性较高的反射系数，而重构后的宽频剖面分辨率明显提高、同相轴连续性更好、断点更清晰、断层更明了，相干体切片精度更高、细节刻画更清晰，可作为裂缝预测的重要技术手段。关键词：分辨率；ADM；裂缝预测；相干体；反射系数中图分类号：P631 文献标识码：A文章编号：1000一1441(2016)05一0737一09 DOI：103969jissrL 1000一1441201605013Researc

4、h and application of high resolution fracture predictiontechnology based on A【M spectral inVersionZHANG Hual，23”，HE Zhenhual“，LI Yalin 3一，ZHANG Xiaobin 3一，HE Guangming 3”，LI Dezhen3t4(1Sn把K倒L乜6DmD珂o，az口耐s RPsP舢打G妞册dpzoi缸io竹，吼删“610059，吼i加；2CDzz昭e o，o一户矗ys洒，吼明“溉i醒叻o，nc行DZogy，吼已行鲥“610059，Ci船；3o户y5icnZp

5、Zorni彻C蝴加以y，吼抛挖qi咒gD，删聍g励gi九卯ri以g CDL酰，CNPC，饶翻”610213，吼i船；4A幻“刀缸锄opysicnZ n如noZo删n对(h抛r，CNPC，C明鲥M 610213，吼i加)Abshact：How to obtain the high resolution data vol啪e and high precision coherence slices have been the key research issues in the exploration of fractured reservoirsTherefore，the application

6、of high res01ution fracture prediction technologybased on舢temating Direction Method(ADM)spectr啪inversion is studied-F姗the basic principle of spectralinversion，the target function of the multi-layer sparse reflection coefficient is derived The high precision reflection coefficient is obtained谢th the删

7、spectrum inversion ofbasis pursuit algorith札The high resolution section谢th the low frequency information is reconstructed by the convolution of the high precision reflection coefficient with the broadband zero phase、，aveletextmcted fmm spectrum simuIation analytical1ethodThe fracture prediction for

8、this high resolution section using the improved eigenvalue algorithm is conductedThe theoretical experiments and actual data applicatian tests show出at fracture p廿diction tech的logy b芒ISt甜on ADM spect姗in、，erSion can get the renection coefficient with the acc啪te polaIity and position aS a11h唧吡mt th血que

9、 for the fractIlre prediction The reconStructed bmadband sections are chmct捌by higher resolution，bettercontinuous e咖ts and nDre clear bre如htS，and the 00h奠叫b()dy slice hasligher accuracy、itll fine reservoir characterization收稿日期：20150925；改回日期：2016一0412。作者简介：张华(1981一)，男，博士在读，工程师，主要从事地震资料处理、反演方法研究工作。基金项

10、目：中国石油天然气集团公司科学研究与技术开发项目(2013E_38一08)资助。Tks research is financially supported by the Science Research and Technology Development Pr两ect of China National Petr01eum Corporation(Grant No2013Ej38一08)万方数据738 石油物探 第55卷KeywoIds：resolution，ADM，fracture prediction，coherent body，refIection coefficient相干体裂缝特征

11、描述是裂缝性油气藏勘探的重要手段之一，而获取用于裂缝预测所需高分辨率数据体以及裂缝参数的求取一直是裂缝油气藏勘探中重点研究的问题。薄层反射系数谱反演13是近些年发展起来的一种能有效分辨薄层、提高分辨率的方法。该方法的主要特点是不需要任何先验信息和数学假设条件，利用反射系数奇偶分解理论和部分频率谱信息反演地层稀疏反射系数，进一步提高分辨低于调谐厚度的地震薄层能力。PURYEAR等4I、CHOPRA等5对谱反演算法理论进行了详细推导，模型测试和实际资料应用结果证明该方法能稳健地分辨薄层厚度。YUAN等6详细分析了谱反演的不适定性，并提出了快速的粒子群LevenbergMarquardt混合反演算法

12、，提高了反演的精度和效率。柴新涛等7利用最小二乘QR分解(LeastSquares QR，LSQR)算法进行谱反演，得到较好的模型测试效果。陈祖庆等8利用地震资料的谱信息和基追踪算法进行稀疏脉冲反射系数谱反演，模型试验结果表明该方法能分辨原始数据中无法识别的薄层。三维相干体技术是20世纪90年代后期兴起的一项十分有效的地震解释技术，该技术主要从相邻地震道相互之间的相干性出发，给出一种定量描述，能突出断层处地震波的变化情况，形成相干值异常区域，从而预测断裂展布。BAHORICH等9提出在振幅上进行互相关相干分析C。算法，但该算法计算出的断点非相干性点模糊，噪声大、能量弱。MARFURT等1 o提

13、出基于相似系数相干C2算法，具有稳定性好、断层刻画精度高的特点，但横向分辨率低的劣势无法改变。GERSZTENKORN等11提出基于本征结构分析的C。算法，利用Hil一bert变换计算相干体，提高其稳定性和抗噪性，但该算法用于高陡复杂地区时，如不考虑倾角则效果不如第二代算法的效果好，如考虑倾角则因计算效率低而无法应用于生产。张军华等1 2。、王玉学等13、刘振峰等1 4|、齐晴等151对C。算法及其优化算法在断层解释、裂缝预测方面进行了讨论和改进，取得不错的应用效果。本文在前人研究的基础上，利用基追踪算法是全局优化算法，能使地震信号的稀疏分解更加稳定的优势1 6|，进行基于基追踪算法理论的交替

14、方向法(ADM)谱反演算法研究，求取高精度反射系数，再利用该反射系数与含低频的宽频子波重构得到分辨率更高的三维地震数据体，在该数据体上利用改进的C。相干算法进行高精度断层识别和裂缝预测，避免了解释人员在断层解释和组合的随意性，使断层解释、裂缝预测的精度大大提高。1 ADM谱反演的高频拓展原理11 稀疏谱反演基本原理根据不含噪声的褶积模型，地震合成记录在频率域可表示为：S(厂)一V矿(厂)R(厂) (1)式中：S(厂)，w(厂)和R(厂)分别表示频率域地震记录、地震子波、反射系数。稀疏谱反演利用频率域地震记录与子波的丰富信息，采用合适的反演算法消除子波影响，从记录中得到所需的反射系数。利用反射系

15、数奇偶分解原理可得到频率域下谱反演目标函数表达式：o一肚Re器一Re附，)+口。Im器一-郴c纠(2)式中：n。和口。分别代表反射系数尺(厂)的偶、奇分量比例。为建立更符合实际地层情况的反射系数模型，必须推导多层反射系数模型下的目标函数，由反射系数奇偶分解原理得到多层反射系数偶、奇分量，表示为：g。(f)一90()一占卜孚卜外孚艿卜孚卜件孚出d(3)式中：T表示反射系数时间厚度re和ro分别代表反射系数，一()的偶、奇分量。对(3)式进行傅里叶变换得到其频率域实部与虚部表达式：ReR(，厂)一I r。()cos7c口()出 Jr (4)ImR(，厂)一I凡(f)sin兀，()出利用公式(2)和

16、公式(4)得到多层反射模型的目标函数表达式为：万方数据第5期 张 华等基于ADM谱反演的高分辨率裂缝预测技术研究及应用 739o饥，r。归怍Re潞一帕_h)cos丌口卜Im滁一_Jhsin兀口d)d厂 (5)式中：&表示时移量。为了便于求解，将公式(5)写成矩阵形式：。，roum乏二2东(6)其中，A。=cos(7cTl)cos(7rTl)cos(7【T1)A。一cos(7cT2，1)cOs(7c瓦，2)cos(7cT2厶)cos(7c丁2厂1)cos(丌瓦2)cOs(7c凡2厶)(7cL) sin(丌T2厂I) sin(7ch2)(7cTl Jf2) sin(7【T2厂2) sin(7cn2

17、厂z)； i i(7【T1) sin(丌T2) sin(7【n2)阢一re(f1) r。(f2) r。(“)T坑一r0(】) r。(f2) ro(N)TR。一R。一Re澈第删1Re器芬水1；Re器等e一qm器等刑1-m滟铬一叽1im器等e吨啪为了降低求解问题时的多解性和不确定性，需要结合实际地震反射系数为稀疏的这一假设条件。因此，加上一个稀疏化约束，将目标函数改为范数求解形式为：。c“，ro，一II：i乏二2姜乏j虻+A II r。+r0忆 (8)式中：A为稀疏因子；Il，表示L。范数；表示Euclidean范数的平方。根据基追踪去噪算法理论，可将目标函数(8)的求解问题转化为：min|舡一6

18、雌+l|x|I。 (9)J伊厶U式中：u为极小权重值；G为复数集合；行为向量x的维数。如果引入变量，9，那么公式(9)等价于公式(10)：min击|I r雌+11工眦：触+r一6 (10)jG“，(，I厶u J式中：m为向量r的维数。公式(10)对应的增广拉格朗日子问题可表示为：min|1 r|；+|x Il，一j，T(缸+r一6)+J伊，Gm I厶U譬ll触+r一6肥 (11)式中：yT表示乘子的共轭转置运算；p为罚参数。利用公式(11)进行基追踪ADM谱反演算法，实现步骤如下。1)令志一。对，矿赋初始值，并给定u，卢常数值，然后进行算法运算，如果满足终止准则，则完成运算，否则，进行步骤2)

19、；2)令x一，j，一矿，求解，的子问题得到12删谱反演算法求解目标函数(6)式)的算法有很多，如共轭梯度法、模拟退火法、匹配追踪法、基追踪法等，这些算法及混合算法各有优劣，其中，基追踪去噪算法由于其数值计算过程中表现出的稳健性、高效性和精确性越来越受到地球物理学者的青睐，本文采用以基追踪去噪理论为基础的ADM谱反演算法来求解目标函数最优解。由于该目标函数(6)式)有多个全局最优解，1一尚芳一c山)，3)令r一，+1，y一少则关于工的极小化问题公式(11)等价于公式(13)：圳圳，+钏触+，L 6一洲(13)那么可以通过公式(14)近似求解来完成对公式(13)的精确求解：f竺孥l|z|。+(矿)

20、T(工一)+去II工一rlI； J伊 L 厶c Jl矿垒AT(触6+，16一号)f1 4、万方数据740 石油物探 第55卷式中：r为大于零的邻近参数，矿为x一时二次项的梯度，则公式(14)通过公式(15)求解得到矿十1：一=Shrink(一小吉)垒max(1一rgl一吉，o)。sign(prg)(15)式中：“。”表示逐元素相乘；Shrink()表示一维收缩算子；sign表示符号函数。4)令x一+1，r一一+1，求解y的子问题得到y+1：j，抖1一矿一稻(触抖1+，16) (16)式中：y为大于零的常数。5)令惫一是+1，重复步骤1)至步骤4)进行迭代运算。13保持低频信息的宽频子波提取为得

21、到保持低频信息的高分辨率资料，需对求取的反射系数与保持低频信息的宽频子波重构。零相位子波是具有较高分辨率的子波，但采用不同主频的零相位子波，主频的变化会使其低频相应发生变化，重构后无法保持提频前的低频信息，所以必须获取能保持低频信息的宽频子波振幅谱，从而得到零相位宽频子波。我们利用谱模拟解析法17从提频前地震数据的振幅谱进行高频拓宽，把该宽频振幅谱作为零相位子波振幅谱，采用傅里叶反变换得到真正具有提频前低频信息的宽频子波。假设地震子波的振幅谱可以表示为下面光滑的解析函数：L硼(，)一尸exp(口。厂”) (17)n=O式中：厂表示频率；L为多项式阶数；n。为依赖于实际地震记录的待定常数。对于确

22、定的L，将计算出的口。(72一o，1，2，L)代人(17)式，可以得到一条光滑的拟合曲线l叫(厂)l，此曲线即是我们所认为的子波振幅谱。利用公式(17)求出提频前地震数据振幅谱作为子波振幅谱(图1a)，然后通过交互扩展该振幅谱高频部分，低频振幅保持不变，得到宽频子波振幅谱(图1b)，以该振幅谱与零相位谱重构经过傅里叶反变换后得到保持低频信息的宽频子波(图1c)。删 腓 k20 40 60 80 100 120厂Hzb图1保持低频信息的宽频子波提取a输人数据振幅谱及拟合曲线；b拓宽高频拟合曲线得到的振幅谱；c拓展变频拟合曲线得到的振幅谱与零相位谱重构得到的宽频子波2 改进型特征值法相干体裂缝预测

23、断裂、裂缝对碳酸盐岩储层有明显的改善作用，便于解释人员利用预测裂缝手段进行裂缝型储层预测。相干数据体技术主要利用相邻地震道的相似性原理，描述地层及岩性的横向不均匀性。图2中，地震信号上半部连续性好的同相轴对应于连续性较好的水平地层、倾斜地层等地质体，下半部连续性差的同相轴对应于连续性较差的断层、裂缝等地质体。采用相干数据体算法，对地震数据体的不连续性进行系统分析，使解释人员能有效识别构造和断层的分布情况，避免解释的随意性。万方数据第5期 张华等基于ADM谱反演的高分辨率裂缝预测技术研究及应用 741入图2相干体技术描述裂缝构造原理为了尽可能压制噪声干扰，提高相干体计算数据的分辨率，GERSZT

24、ENKORN等11提出了第三代特征值相干体计算法。首先利用MARFURT等10在第二代多道相似性相干体计算中定义的椭圆和矩形分析窗口(图3)中，以坐标为(z，y)的分析点为中心，时窗中心一扎位，则时间窗口内含2忌+1个元素的分析道与窗口内的，道数据排列成矩阵：【，= (18)式中：而。表示第歹道在时间f一优&一声她一qM沿着视倾角(户，q)插值后的振幅值。图3分析窗口示意a以分析点(红点)为中心的椭圆分析窗口，主轴为n，径向轴为6，主轴方位角为仇Ib以分析点(红点)为中心的矩形窗口，主轴为n，径轴为6，主轴方位角为他那么矩阵【，的第优列厦。一而，。，牙：。，研，。T是分析体内每个地震道在时间一

25、m&一乡纰一蚴(歹一1，2，)的振幅值，则沿着一对视倾角(户，q)，时窗中心在一优&的分析窗口内的2忌+1个采样点的数据协方差矩阵为：0(声，g)一UuT一，r女牙2，酒椭m=矿K芒 厶牙，属mm=rrK ，(19)设凡为分析窗口内地震数据协方差矩阵公式(19)的特征值，且有A。AzAJ，则在以分析点为中心的分析窗口内，沿视倾角(户，q)的，道C。特征值相干算法公式可表示为：1 C3(夕，g)一粤L (20)扎则在以分析点为中心的分析窗口内，沿视倾角(乡，go)的J道C。特征值相干算法公式表示为：C3一C3(户一O，gO) (21)则沿着所有倾角的C3特征值相干计算算法公式可表示为：C3一ma

26、xC3(户，q) (22)一彬一彬一彬MK=li一枞一万方数据742 石油物探 第55卷第三代特征值相干算法C3在子空间中进行计算，比前两代相干算法能更好地压制噪声，但计算量大、效率慢，影响了实际应用时的经济效益。为此，我们对该算法进行改进，利用基于相似性算法的高效率性估算各点的视倾角和方位角，在大于10倍的分析窗口内计算视倾角和方位角的平均值(声，g)；再在以分析点为中心的分析窗口内，沿着区域平滑后的视倾角(石，石)计算地震数据协方差矩阵的特征值及相干值：己一c3(石，石) (23)数值实验及实际资料应用31 断层模型谱反演试验利用一个断层模型(图4a)进行ADM谱反演实验，图4b是对图4a

27、断层模型正演得到30 Hz主频的合成记录，图4c为本文方法反演的反射系数，图4d为主频为50 Hz重构记录。分析图4可知，ADM谱反演算法能求出极性和位置准确性较高的反射系数，而且通过重构能得到断层分辨更清晰的合成记录。图4断层模型反演实验a断层模型；b 30 Hz合成记录；c谱反演反射系数；d 50 Hz重构合成记录32 Ma咖。吣i模型谱反演试验利用SEG提供的Marmousi模型进行谱反演反射系数实验，以验证本文方法对纵、横向均存在变化的地震数据反演的适应性。图5a为MarmOusi模型反射系数，图5b是35 Hz主频的合成记录，图5c为本文方法反演的反射系数。分析图5可知，ADM谱反演

28、算法能精确求出含断裂带、深部盐丘两侧高速体等复杂数据体的反射系数，证明了该方法的正确性和对复杂模型较好的适应性。33实际资料应用将基于ADM谱反演重构后的高分辨率实际地震数据用于裂缝预测中，以验证该方法在储层预测中的有效性。对图6a叠前时间域偏移剖面目的层进行ADM谱反演实验，得到反射系数剖面(图6b)，该反射系数与测井系数(蓝色曲线)匹配较好，各层位置和极性求取较佳。用提取的宽频子波进行高分辨率重构，得到高分辨率剖面(图6c)，其分辨率明显比图6a的高，且浅、中、深层的同相轴(如白框处)连续性更好、断点更清晰、断层更明了，更有利于后续解释人员进行裂缝预测。分析高分辨率处理前、后的振幅谱(图7

29、)可知，提频后频带得到明显拓宽，而且由于使用保持低频的宽频子波重构，提频后振幅谱的低频信息得到完整保留，为后续反演等解释提供了丰富的频带信息。为了分析高分辨率处理前、后的效果，分别提取了20 ms的C。相干切片(图8)，分析图8可知，利用ADM谱反演重构的高分辨率剖面比处理前剖面能获取分辨率更高的相干切片效果，且切片上裂缝走势更清晰、裂缝刻画细节更精细。利用过井成像测井资料(图9)进一步验证高分辨率处理前、后裂缝发育所在断层的真实性。从图9中可以看出，在2 2782 283 m段有1个裂缝发育段，在2 279 m处裂缝发育，通过合成记录精细标定裂缝发育带相当于须三底向下20 ms，对最小负曲率

30、属性数据体进行须三底向下20 ms切片(图loa)，与过井相干剖面(图10b)对比后可知，过井点处发育黑色异常，也即存在裂缝发育带，这也印证了提频后剖面(图6c)箭头所指位置存在明显断层，证明本文方法得到的高分辨率剖面为后续裂缝预测提供了精度更高的依据。万方数据第5期 张华等基于ADM谱反演的高分辨率裂缝预测技术研究及应用 743图5 Mamousi模型反演实验a Mamousi模型反射系数；b 35 Hz合成记录；c谱反演反射系数691 741 791 84l 891 94l道号a691 741691 74l 84l 89l 941 99I道号C841 89l 941 99l道号b图6提频前

31、实际地震数据(a)、谱反演反射系数剖面(b)和宽频子波重构后高分辨率剖面(c)万方数据744 石油物探 第55卷a图7提频前(a)、后(b)数据振幅谱a图8提频前(a)、后(b)数据须三底向下20 ms相干切片“然伽一5 密度(gcTTr：J 浅删同jU 4H半 删外扎隙 删什渗透率 试AP【 (n-ml h 岩性特征 度， 10Lun!) 油储 有测试1 1jO25中子一02 2 000 地 0 20 D00l 100000含水饱和 效Jr径cm 深删向电阻半m 质 嬲剁强露 岩心孔隙 井心渗透牢度 厚 射 24 嘉波7(岫鼍3fnm1 分 度， (10urll2) 层 度 孔 情况2 2

32、000 展 1 20 D 0 0l 100 000 100 0 段户隧黼?。j曩基鬃巍翻r豳 2+乏 坼 淄 2 220 ) 2 279+3)须三一*g缚 7 2 230j 至2 240- 鋈。j群ii!i!i鞫 t) 三， ；誊巍攀2 280颁 -尹霹譬懋毖鹭牙 ， 事_ 乏+I，j鎏三；醯；14xr8懑群蟊獭溺 。i。j2 250 卜2 1602 2332 287m羹一 2 28l 2 5m3d匿 器：芦_ 产气1 270-105m3dI6 。l 2 280一 熏B粼鋈l隧剿2 2822 290 醮黼糍罄誊鬻鲺ij_警2 300：二 | 4S-C 叠墓鲼端珏2 283)C 锌一 *冒MK

33、InP 露黼 夕万方数据第5期 张华等基于ADM谱反演的高分辨率裂缝预测技术研究及应用 7454 结论a图10须三底向下20 ms最小负曲率属性(a)和过井数据相干剖面(b)本文研究了ADM谱反演高分辨率技术在裂缝预测中应用情况。模型试验和实际资料应用结果表明，利用ADM谱反演算法求取精度高、极性位置正确、能量强弱分明的地层反射系数；与谱模拟解析法提取的宽频子波褶积重构，得到分辨率高、同相轴连续性好、断层清晰、低频信息丰富的剖面；在此剖面上能得到精度更高、细节刻画更清晰的C。相干切片，便于断层和裂缝的有效识别，避免了解释人员在断层解释和组合时的随意性，大大提高了断层解释、裂缝预测的精度。参考文

34、献1 PORTNIAGUINE 0，CAsTAGNA J PInversespectral decompositionJExpanded Abstracts of74由A1111ual hternat SEG M唔，2004：1786一17892 PORTNIAGUINE 0，CASTAGNA J PSpectral inversion：lessons from modeling and王bonesville casestudyJExpanded Abstracts of 75th Annual Intermt SEG Mtg，2005：163816413 CHOPRA S，CA盯AGNA J

35、 P，PORTNIAGUINE0Thin?bed reflectivity inversionJExpanded Abstracts of 76小Annual htemat SEG Mtg，2006：205720614 PURYEAR C I，CASTAGNA J PLayer_thicknessdetennimtion and stratigraphic interpretation usingspectral inversion：theory and applicationJ一physics，2008，73(2)：R37一R485 CHOPRA S，CASTAGNA J，XU YRelat

36、ive acoustic impedance application for thinbed reflectivity inversionJExpanded Abstracts of 79【h，岫ual Intemat SEG Mtg，2009：355435586YUAN S Y，WANG S X，TIAN NSwaml intelli-gence optimization and its applicatiomn in geophysical data inversionJApplied Geophysics，2009，6(2)：1661747柴新涛，李振春，韩文功，等基于LsQR算法的谱反

37、演方法研究口石油物探，2012，51(1)：1118CHAI X T，LI Z C，HAN W G，et aL Spectralinversion method analysis based on the 15QR algoritJGeophysical Prospecting for Petr01e啪，2012，51(1)：11一188陈祖庆，王静波基于压缩感知的稀疏脉冲反射系数谱反演方法研究J石油物探，2015，54(4)：459466CHEN Z Q，WANG J BA spectral inversion method of sparsespike reflection coeffic

38、ients based oncompressed sensingJGeophysical Prospecting forPetroleum，2015，54(4)：4594669BAHORICH M S，FARMER S L 3一D seismic coherency for faults and stratigraphic features：the coherence cubeJThe Leading Edge，1995，14(10)：1053105810 MARFURT K J，KIRLIN R L，FARMER s L，et a13一D seismic attributes using a

39、 semblance-based coherency algorithmJGeophysics，1998，63(4)：1150一116511 GERS2rTENKORN A，必LI心URT K JEigenstructurebased coherence computations as an aid to 3一Dstructural and stratigraphic mappingJGeophysics，1999，64(5)：1468147912张军华，王月英，赵勇C3相干体在断层和裂缝识别中的应用J地震学报，2004，26(5)：560一564ZHANG J H，WANG Y Y，ZHA0

40、 YApplication ofthe thhd generation of coherent cube in recognizingfaults and fracturesJActa seismologica Sinica，2004，26(5)：56056413王玉学，韩大匡，刘文岭，等相干体技术在火山岩预测中的应用J石油物探，2006，45(2)：192196WANG Y X，HAN D K，LIU W L，et a1The appli-(下转第763页)万方数据第5期 侯伯刚等变差函数的参数和井数对随机反演精度影响的分析 7631213141516刘勇，雷茂盛地震资料随机反演方法J石油物

41、探，1998，37(增刊A)：3l一35LIU Y，LEI M SA stochastic inversin method ofseismic dataJGeophysical Prospecting for Petroleum，1998，37(S1)：3l一35印兴耀，孙瑞莹，张广智基于分形高频初始模型和低频先验信息的物性参数随机反演J石油物探，2014，53(5)：537544Y1N X Y，SUN R Y，ZHANG G ZStochastic inversion of reservoir physical property parameters basedon highfrequenc

42、y initial nlodel from fractal and low_frequency prior infolmationJGeophysical Prospecting for Petmleurn，2014，53(5)：537544刘文岭油藏地球物理学基础与关键解释技术M北京：石油工业出版社，2014：158175LIU W LReservoir geophysics foundation and keyinterpretation technologyMBe巧ing：Petroleum Industry Press，2014：158175侯伯刚，刘文岭，罗娜地震反演中测井数据的预处

43、理J物探与化探，2009，33(3)：33l一336HOU B G，LIU W L，LU0 NThe preprocessing ofwell log data for seismic inversionJGeophysicaland(kochemical Exploration，2009，33(3)：331336王玉学，侯伯刚，牛滨华地震储层特征反演技术在冀中南马庄地区的应用J现代地质，2003，17(增刊A)：219223WANG Y X，HOU B G，NIU B HApplication ofseismic reservoir inversion in the area of NMZ

44、ofHebei ProvinceJGeoscience，2003，17(S1)：21922317杨丽兵，张虹，谯述蓉随机地震反演在川西坳陷DY构造储层预测中的应用J石油物探，2010，49(4)：373379YANG L B，ZHANG H，QIA0 S R Application ofstochastjc seismic inversion on reservoir predictionof DY structure in western Sichuan depressionJGeophysical Prospecting for Petr01e眦，2010，49(4)：37337918侯伯

45、刚，刘文岭，王经荣港东一区井震联合剩余油富集规律研究J石油物探，2009，48(2)：175180HOU B G，LIU W L，WANG J R Study on the enrichrnent regularity of residual oilin well and seismicjoint in GD yiqu口GeophySical Prospecting for Petroleum，2009，48(2)：17518019高云，朱应科，赵华叠前同时反演技术在砂砾岩体有效储层预测中的应用J石油物探，2013，52(2)：224228GA0 Y，ZHU Y K，ZHAO HApplica

46、tion of prestack simultaneous inversion technique in effectivereservoir prediction of sandgravel bodyJGeophysical Pmspecting for Petroleum，2013，52(2)：224228(编辑：陈 杰)-+一一-+_一+_一-+_一-卜一卜一卜一-卜一-卜一卜一一一-+_一卜-卜-卜-卜一卜-卜*卜-+一-卜一卜-卜一+-+-+-+-+一+-+-+一+一+-卜一+-_卜-卜一卜一卜-卜一卜-卜w卜-卜一卜(上接第745页)cation of coherence tech

47、nology in the volcanic reservoir predictionJGeophysical Prospecting for Petr01eum，2006，45(2)：19219614刘振峰，曲寿利，孙建国，等地震裂缝预测技术研究进展J石油物探，2012，51(2)：191198LIU Z F，QU S L，SUN J G，et a1Progress of seismic fracture characterization technologyJ(；eophysicaI Prospecting for Petroleum，2012，51(2)：19119815齐晴，曲寿利，李振春，等也门S2区块叠前、叠后裂缝综合预测方法应用研究J石油物探，2014，53(5)：603608QI Q，QU S L，LI Z C，et a1Fracture prediction byusing prestack and post_stack seismic technique