Facies based Bayesian simultaneous inversion technology and its application: Case study of the north section of No.5 fault zone in Shunbei area, Tarim Basin

  • Fei CAO , 1, 2 ,
  • Zhiqiang LU 1, 2
Expand
  • 1. Research Institute of Exploration and Production,Northwest Branch of SINOPEC,Urumqi 830011,China
  • 2. Experimental Testing Technology Center,Northwest Branch of SINOPEC,Urumqi 830011,China

Received date: 2022-03-03

  Revised date: 2022-05-16

  Online published: 2022-09-28

Supported by

The National Natural Science Foundation of China: Enterprise Innovation and Development Joint Fund Project(U19B6003)

Abstract

Fault-karst carbonate reservoir in Shunbei area of Tarim Basin is controlled by deep strike-slip faults and forms fault-karst system. The reservoir space mainly includes hole and fracture, its strong anisotropism aggravates the complexity of the seismic response characteristics of the reservoir. The high-quality reservoirs in the fault-karst body in this area have different burial depth, which is not conducive to the establishment of low-frequency models in traditional inversion. The facies based Bayesian simultaneous inversion technology combines Bayesian classification with pre-stack simultaneous inversion, divides different facies based on multi-elastic parameters such as P-wave and S-wave velocity and density, and conducts in-depth trend analysis for each phase to establish initial model.Compared with the traditional inversion technology, this technology not only improves the inversion accuracy, but also increases the stability of the density inversion. Taking the carbonate fault-karst body in the northern section of the No. 5 fault zone in the Shunbei area as the research object, combined with the actual production situation, two facies, fractured-cavity limestone and tight limestone are divided by elastic parameters, and then are conducted with depth trend analysis and inversion. Through the single fracture-cavity equivalent model test and practical application analysis, the density data obtained by the Facies based Bayesian simultaneous inversion is highly consistent with the reservoirs drilled by Wells W3 and W3C in the northern section of the No. 5 fault zone, which verifies the applicability and reliability of the inversion technique in the study area and the reliability of the results.

Cite this article

Fei CAO , Zhiqiang LU . Facies based Bayesian simultaneous inversion technology and its application: Case study of the north section of No.5 fault zone in Shunbei area, Tarim Basin[J]. Natural Gas Geoscience, 2022 , 33(10) : 1702 -1711 . DOI: 10.11764/j.issn.1672-1926.2022.05.003

0 引言

近年来,在塔里木盆地顺北地区相继获重大油气发现,揭示了走滑断裂带对奥陶系碳酸盐岩断溶体具有明显“控储、控藏、控富集”作用。前人结合地质理论、地震技术方法等已开展了深入研究,鲍典等1应用相干、蚂蚁体、能量体、波阻抗反演等数据描述断裂的横向分隔性、断溶体的纵向结构变化与纵向分隔性及断溶体内部结构分隔性,阐述了纵、横向分隔是引起局部油水分异与生产动态差异的主要原因;伍齐乔等2、孙东等3认为断溶体油藏发育的必要条件是发育区域大型走滑—拉分断裂,而且断裂交会裂缝发育区对储层改造作用最为显著,是寻找高效井点的理想场所;刘军等4、马乃拜等5、刘宝增等6、李宗杰等7利用数值模拟、物理模拟技术,明确不同尺度断裂的地震响应特征,也从断溶体目标采集与成像出发,建立了地震识别模式,形成“张量属性定轮廓、振幅属性定边界、融合雕刻定体积”的圈闭落实与描述技术体系,在探区实现了油气重大突破。
5号断裂带是顺北一区内发现并落实的最长断裂带,累计长度达86 km,并且从南到北,由北东向逐渐过渡到北西向。目的层为奥陶系中下统,沉积相为开阔台地相,岩性以灰岩、泥质灰岩为主(图1),储层以裂缝型为主、次为溶蚀孔洞型。顺北5井是在该断裂带北段所钻的一口油藏探井,完钻层位鹰山组(鹰山组顶深7 434 m),该井钻进至7 514.30 m出现漏失,录井见油气显示,表明顺北5井断裂带成藏条件好,是油气富集区带8。顺北油气田1号、5号主干断裂带奥陶系油气藏的不同油气充注过程则说明顺北地区断溶体成藏的差异性9
图1 研究区构造位置与主干断裂分布图(a)及综合柱状图(b)

Fig.1 Structural location and distribution of main faults (a) and comprehensive histogram (b) in the study area

尽管针对5号断裂带断溶体的研究取得了很多成果,但在其北段勘探过程中仍存在一个严重问题,即对于1号断裂带断溶体非常有效的、储层预测吻合率近100%的内部结构三参数(即将剩余波阻抗、相干处理及断层自动拾取等数据融合后刻画断溶体结构与展布)表征方法不再适用,在地震反射特征以弱串珠和杂乱反射为主的5号断裂带北段储层预测吻合率仅为50%~60%,例如顺北5-3井酸压建产及后期累计产能与最初设计严重不符,说明利用叠后地震资料很难实现该区强非均质性的断溶体定量刻画。考虑到5号断裂带北段断溶体储层非均质性强,出现“串珠”、杂乱、弱、片状等多种反射特征,且钻井有直井、斜井、水平井等多种类型,目标埋深差异较大(图2),在这种地质背景和资料条件下,叠后地震技术的储层刻画精度相对低,且由于难以建立合理的低频模型,传统叠前反演方法也很难获得理想的断溶体预测结果。鉴于此,为提高5号断裂带北段断溶体的预测精度,本文对相约束贝叶斯叠前地震预测技术进行了研究,在贝叶斯框架下,引入“相”和深度趋势分析,针对不同的“相”建立纵波阻抗与横波阻抗、纵波阻抗与密度的关系,解决了传统反演中将所有岩相只建立一个关系的问题。在此基础上,对不同“相”进行相应弹性参数的深度趋势分析,并作为不同“相”的独立低频趋势模型,避免了传统低频趋势模型中不同岩相采用相同趋势的问题,应用期望最大化算法同时求解“相”和弹性参数。相较于传统方法,此技术不再依赖于初始井模型,突破了传统叠前反演低频模型的局限性,并将岩石物理分析贯穿整个反演过程,既增强了密度反演稳定性,同步获得纵横波速度、密度等弹性参数与“相”,也可以有效保证反演结果的一致性10-15
图2 研究区连井地震剖面

Fig.2 Seismic profile of wells in study area

1 相约束贝叶斯同时反演原理与方法

1.1 “相”及深度趋势的定义

相约束贝叶斯同时反演方法与传统叠前同时反演的主要不同是引入了“相”及深度趋势分析的概念。
在相约束贝叶斯同时反演中,“相”是与地质解释相关的岩性及所含流体,并且在弹性域中能够区分。基于相的反演,其目的是利用地震数据约束岩相,这就要求岩相在弹性域,尤其是在纵横波阻抗方面必须具有一定的区分度。这种定义与理想化的地质定义不完全匹配,主要是岩相划分的数量通常受到弹性域重叠程度的限制。因此,“相”应称为弹性相16。在贝叶斯分类框架下,引入“相”,能够将岩石物理分析贯穿整个反演过程,纵横波速度与密度及三者共同建立的“相”相互约束进行叠前同时反演,最终同步获得纵横波速度、密度等弹性参数数据体与“相”体。这不但保证了反演结果的一致性,更好地获得有利储层与各岩石物理参数的对应关系,也增强了密度反演结果的稳定性。与传统反演不同,所用井资料仅使用岩石物理特征进行分析,只是对反演结果进行质控而不直接参与反演计算。
传统地震反演过程中需要建立速度与密度的低频模型,而速度与密度等弹性参数都随着纵向压实和横向沉积的变化会产生时变与空变,很难建立准确的低频模型,从而导致强非均质储层预测的精度不高。因此,引入深度趋势分析可以避开低频模型建立的难题。深度趋势来源于碎屑岩,用于研究深度和压实作用对砂岩和泥岩孔隙度的影响,也能通过速度—深度趋势识别异常岩性或成岩区域(如碳酸盐岩的后期改造)17。在相约束贝叶斯同时反演中,深度趋势是建立“相”划分的基础,不同的“相”建立不同的深度趋势,具有2个作用。一是进一步调整相划分的标准。在实际反演过程中,不同“相”的各岩石物理参数曲线随深度的变化趋势拟合需要不断调整,随着深度趋势拟合的调整,各“相”在不同地球物理参数交会图上的分布区域也会发生相应的变化。调整的标准是使得各“相”在不同地球物理参数交会图上的分布区域能尽量地区分开来,调整的质控依据是使反演得到的井旁道岩相及各属性与实际井曲线越吻合越好。二是通过深度趋势拟合,建立不同“相”的纵横波速度、密度等不同弹性参数的深度变化趋势,作为相约束贝叶斯同时反演的初始模型,相当于传统叠前反演中的低频模型。这是深度趋势分析最重要的作用,也是与传统叠前反演差异最大的地方。针对不同“相”建立的多个独立弹性参数趋势模型,很好地避免了传统反演单一低频模型中不同岩相比例量化问题,解决了复杂地质背景下低频模型建立的难题。

1.2 方法的实现

在传统的叠前同时反演的算法中,依据FATTI等18方程推导,反射系数:
R P P θ = a R A I + b R S I + c R ρ
式中: R A I = V P 2 V P + ρ 2 ρ
           R S I = V S 2 V S + ρ 2 ρ
           R ρ = ρ 2 ρ
           a = 1 + t a n t 2 θ   b = - 8 V S V P   2 s i n 2 θ
通过与子波W(θ)褶积,将式(1)中的反射系数R pp(θ)转换为合成地震记录S pp(θ):
S P P θ = a ,   2 W θ L n A I   + b 2 W θ L n S I + c 2 W θ L n ρ   
简化的矩阵形式为:
S = J * L n z
其中: S是近、中、远等叠加地震道的块列向量。
考虑到地震信号的噪声分量及地震反演中岩石性质的不确定性等,同时反演可被视为一个统计问题,地震反演可以重新定义为贝叶斯问题。
贝叶斯反演实际上是一种使用概率约束的反演问题,贝叶斯理论是一种用于计算条件概率分布的统计学理论,表达式如下:
p m | d = p d | m p m p d
式中: m 为模型参数; d为观测数据;p m | d )是后验概率分布;p d | m )是观测数据的似然概率分布;p m )是模型参数的先验分布;p d )是边缘概率分布;一般为一个常数。所以可以将公式简化为:
p m | d p d | m p m
贝叶斯反演方法一般表示为后验分布等于似然函数与先验信息相乘,从式(5)可以看出,贝叶斯理论的关键在于先验信息和似然函数的构建,先验信息能否客观真实地反映真实地下介质的特性,似然函数的选取是否合适,二者决定了待反演参数的结果是否准确可靠。
先验信息是用来约束反演问题的,不同的先验概率分布对反演问题的约束效果也是不同的,因此需要选择合适的先验信息,这关系到反演结果的精度,甚至影响反演结果的正确与否。这也说明贝叶斯反演问题其实是依赖于先验模型的建立,一般来说反演的初始先验信息都是出现在真实值周围,先验分布多来源于测井信息,可能包含研究人员的经验知识,带有很强的主观人为性。
贝叶斯反演问题一般都是模拟出后验分布后,以最大后验概率解作为反演解,而以最大后验分布概率的解是把观测数据与模型参数都看作随机变量,求得的是概率分布,是一种不确定意义下的解。
因此,在相约束贝叶斯同时反演中,贝叶斯定理可以写成:
π Z | S r e a l L S r e a l | Z P ( z )
式中: π Z | S r e a l是后验分布; L S r e a l | Z P ( z )是似然函数;P(z)为先验分布。
先验分布Pz)可以从AI的线性深度趋势、AI与SI和AI与ρ之间的交会图中得到(通常来自井数据)。这些趋势拟合可以表示为一个多正态先验分布Pz),形式为:
P ( z ) ~ e x p - 1 2 Z - Z o T C p - 1 Z - Z o / C p - 1 / 2
式中: C p是描述阻抗之间的方差和相关性的协方差矩阵; Z o为初始模型。
再次使用贝叶斯定理,由于“相”的加入,后验分布的等式改写为:
π Z , F | S r e a l L S r e a l | Z P Z | F P F
P Z | F P F分别为“相”的先验分布、“相”的概率密度,其先验分布P Z |F)与式(7)非常相似,区别在于先验均值 Z o和协方差矩阵 C p现在取决于“相”:
P Z | F = e x p - 1 - 2 ( Z - Z o ( F ) ) T C p F - 1 Z - Z o F / C p F - 1 2
针对最终反演目标函数而言,相约束贝叶斯叠前同时反演与传统叠前同时反演不同。
π ρ s e i s m i c L s e i s m i c ρ p ρ f a c i e s p f a c i e s
q ρ s e i s m i c L s e i s m i c ρ p ρ
式(10)式(11)分别是以密度反演过程为例的相约束贝叶斯叠前反演、传统叠前反演的目标函数。
在相约束贝叶斯同时反演中,增加了“相”的先验信息,从而提高了反演结果的稳定性。
图3为相约束贝叶斯同时反演的实现过程,为了正确捕获地震反演问题的物理特性,对每个“相”的“相”和阻抗进行反演,并且每个“相”的岩石物理模型和低频背景模型都必须参与迭代优化,这种优化过程是期望最大化的一种形式。在相约束贝叶斯同时反演中,利用角度叠加地震数据,综合地质、测井等信息,通过岩石物理分析的各种不同岩相与各岩石物理参数之间的交会关系,采用贝叶斯分类技术对不同的岩相(如储层相、干层、非储层等)进行“相”划分,然后将不同“相”的岩石物理参数特征和概率统计关系代入到叠前同时反演过程中,在期望步骤中,最初给定阻抗,从地震中反演出“相”,然后在最大化步骤中,确定“相”,并从地震中反演出阻抗, 持续下去直到收敛,最终直接反演出各弹性参数和各种“相”的分布。
图3 相约束贝叶斯同时反演过程示意

Fig.3 The facies based Bayesian simultaneous inversion process diagram

1.3 模型分析

为了验证相约束贝叶斯同时反演算法的有效性,设计了一个如图4(a)所示等效洞缝储层模型,模型分别由纵波阻抗、纵横波速度比及密度表征,并划分了洞缝灰岩相与致密灰岩相。对模型分别进行传统叠前同时反演和相约束贝叶斯同时反演,在传统叠前同时反演中使用了一个低频背景模型,相约束贝叶斯同时反演中分别用致密灰岩和洞缝灰岩的深度趋势来代替低频模型,反演结果如图4(b)和图4(c)所示。从反演结果对比来看,相约束贝叶斯同时反演结果与模型更加吻合,边界更加清楚,尤其密度反演结果更精确,并同步获得可直接用于地质解释的储层的“洞缝灰岩相”体。
图4 基于等效缝洞模型的反演结果对比

(a)等效缝洞模型;(b)模型的传统叠前同时反演结果;(c)模型的相约束贝叶斯同时反演结果

Fig.4 Comparison of inversion results based on equivalent fracture-cavity model

2 实际应用及效果

将相约束贝叶斯反演方法应用于5号断裂带北段断溶体的预测,实现过程包括地震资料优化处理、相的划分与深度趋势的建立、井震标定和子波提取及反演运算、质控等关键技术环节,其中相的划分与深度趋势的建立不同于传统叠前反演,其他环节与传统叠前同时反演过程相似,本文不再赘述。研究区在目的层段(4.5 s附近)能量很强、信噪比较高,入射角最大可到40°,满足叠前反演的资料要求,选取5°~15°、15°~25°、25°~35°共3个限角度叠加地震数据用于相约束贝叶斯反演。

2.1 相的划分及标准的确立

通过薄片鉴定,研究区目的层岩性主要为泥(微)晶灰岩、亮晶—微晶颗粒灰岩、偶见亮晶藻黏结灰岩及少量云质灰岩等,原生孔隙不发育。为了在研究过程中保持岩性与弹性参数的一致性,保障后续建立“相”划分的同一标准及相约束贝叶斯同时反演的基础数据可靠,必须在岩石物理分析和参考地质录井资料的基础上,利用电性特征重新进行岩石类型的划分,最终确定目的层岩性主要为灰岩和白云质灰岩及少量灰质云岩。
在岩性确定的基础上,利用弹性特性(纵横波速度、密度、纵横波速度比等)、物性特征(孔隙度等)及参考录井、薄片、试验分析、试油等资料,对不同的岩石类型进行“相”的划分及其标准的建立。“相”划分是一个迭代过程,通过弹性参数的交会分析,不同的相在各个弹性参数空间上的差异来获得各个岩相的先验概率密度分布,并与实钻结果反复对比,确定有效值域。同时,在反演过程中,结合研究区生产实际情况,储层表现为低孔低渗特征,为了简化碳酸盐岩储层缝洞多重孔隙介质的复杂性,将碳酸盐岩地层划分为缝洞灰岩相与致密灰岩相。图5是W3井、W5井、W12井经过迭代确定不同“相”后的弹性特征参数的交会图,图5中紫色点为缝洞灰岩相、蓝色点为致密灰岩相,并以孔隙度1.5%为下限,明确储层与致密围岩的弹性参数变化范围,建立相划分的标准(表1)。
图5 不同相的弹性参数分布特征

Fig.5 Distribution characteristics of elastic parameters of different facies

表1 相划分标准

Table 1 The standard of facies division

序号 纵波速度/(km/s) 横波速度/(km/s) 密度/(g/cm3 有效孔隙度/%
1 碳酸盐岩 致密灰岩相 6.00~6.50 3.20~3.50 2.67~2.70 <1.5
2 缝洞灰岩相 5.10~6.40 3.00~3.47 2.50~2.70 >1.5

2.2 深度趋势的建立

图6为研究区目的层段碳酸盐岩缝洞储层及致密灰岩的纵、横波速度及密度等弹性参数的时间域深度趋势,从图6中可以更清楚地看到储层与致密层围岩的弹性参数差异随深度(时间)增大而减弱,4 400~4 500 ms的100 ms范围内,储层低纵横波速度、低密度特征明显,这可能与走滑断裂形成及演化的机制有较大关系19,走滑断裂由基底向盖层递进传播的发育过程中,纵向上呈现出逐渐向上散开,形成具有一定宽度的、由细小相互叠置不易区分的断裂构成的变形带。图7为目的层储层与致密层的纵横波速度、密度的深度趋势模型及迭代后的混合趋势模型,混合趋势模型为相约束贝叶斯同时反演所用模型。从图7看,由深度趋势分析得到的储层与致密层的纵横波速度与密度的初始模型,横向稳定、纵向上变化符合随深度加大,纵横波速度和密度增大的总体规律;而在经过迭代的混合趋势模型中,纵横波速度与密度的模型呈现纵横向较强的非均质性,储层为蓝色,表现为高纵横波速度、高密度背景下的低值特征,与地震剖面的波组特征有很好的对应关系,而且在井未钻遇的深层同样具有合理的变化趋势。
图6 不同相深度趋势

Fig.6 The depth trend of different facies

图7 不同相的深度趋势模型

Fig.7 Depth trend model of different facies

2.3 反演效果

图8分别为W11井—W3井—W3C井—W5井—W12井连井的纵波阻抗[图8(a)]、横波阻抗[图8(b)]、密度[图8(c)]及相[图8(d)]的反演剖面,暖色调代表灰岩缝洞储层,蓝色代表致密灰岩。储层的纵波阻抗范围为14 000~16 800 (m/s)·(g/cm3)、横波阻抗范围为8 000~9 100 (m/s)·(g/cm3)、密度范围为2.58~2.63 g/cm3。从图8中可以看出,反演结果对串珠地震响应的储层刻画更加清楚,储层的空间位置更加准确,浅层串珠小、成群分布,深层串珠大、呈片状。纵横波速度、密度及相反演结果的一致性非常好,对后续储层解释不会出现相互矛盾,相约束贝叶斯同时反演的优势得以体现,尤其“相”剖面,能够清楚展示储层的分布规律,其中暖色调代表储层,与井点吻合较好,可以直接作为储层的解释成果。
图8 W11井—W3井—W3C井—W5井—W12井连井相约束贝叶斯反演剖面

(a)纵波阻抗剖面;(b)横波阻抗剖面;(c)密度剖面;(d)相剖面

Fig.8 Interwell phase constrained Bayesian inversion profile from Wells W11-W3-W3C-W5-W12

图9是连井优化处理叠前时间偏移地震剖面[图9(a)]及相约束贝叶斯同时反演的密度剖面[图9(b)],密度异常对强“串珠”形态刻画效果很好,与对应的PSTM地震剖面“串珠”形态吻合,并且落实了储层发育的具体位置,可以有效刻画规模较大的断溶体。W3井对应地震剖面上的弱“串珠”反射在反演的密度剖面上几乎没有异常,验证其产能不理想的原因是储层不发育,而W3C井在地震剖面上钻遇的强“串珠”,发生放空漏失,放空漏失段及气测高异常段与密度剖面中明显的低值异常对应,与钻井钻遇储层及漏失情况吻合。钻井深度为7 613~7 619 m、7 651~7 659 m的2段,虽有荧光显示、全烃达10%左右,可能因为埋深大(超过7 000 m),并受地震分辨率的限制,在反演剖面上没有明显响应。
图9 W3井—W3C井连井精度分析

Fig.9 Precision analysis of connecting Wells W3-W3C

3 结论

(1)深度趋势分析和“相”的引入,是相约束贝叶斯同时反演有别于传统叠前反演的核心技术, 作为一种先验信息,不同的“相”将纵横波速度、密度等参数限制不同的范围进行迭代,而且对不同的“相”建立相应的纵波阻抗与横波阻抗、纵波阻抗与密度的关系,解决了传统反演中将所有岩相只建立一个关系的问题;深度趋势与常规反演中的低频模型作用相同,只是不再是单一的低频模型,而是不同的“相”具有各自的深度趋势,即各自独立的低频模型,避免了低频趋势模型中不同岩相比例量化问题。最终提高了反演密度的稳定性及反演结果与地质目标的一致性。
(2)相约束贝叶斯同时反演技术在塔里木盆地顺北地区5号断裂带北段的应用中,对强“串珠”和杂乱反射响应的断溶体有较好的预测效果,能够有效降低地震放大效应,尤其在井少区块的断溶体预测中可以广泛应用。
(3)相约束贝叶斯同时反演依赖于地震资料品质(分辨率和信噪比等),在反演过程中应结合地震相进行岩石物理划分“相”,明确弱反射和片状强反射与“相”的对应关系,可进一步提高反演的可预测性。
1
鲍典,张慧涛.塔河油田碳酸盐岩断溶体油藏分隔性描述方法研究[J].新疆石油天然气,2017,13(1):25-30.

BAO D,ZHANG H T.Description of the separation in fault-karst carbonate reservoirs in Tahe Oilfield[J].Xinjiang Oil & Gas,2017,13(1):25-30.

2
伍齐乔,李景瑞,曹飞,等.顺北1井区奥陶系断溶体油藏岩溶发育特征[J].中国岩溶,2019,38(3):444-449.

WU Q Q,LI J R,CAO F, et al.Characteristics of fault-karst carbonate reservoirs in the Shunbei No.1 well block, Tarim Basin[J]. Carsologica Sinica,2019,38(3):444-449.

3
孙东,杨丽莎,王宏斌,等.塔里木盆地哈拉哈塘地区走滑断裂体系对奥陶系海相碳酸盐岩储层的控制作用[J].天然气地球科学,2015,26(增1):80-87.

SUN D,YANG L S,WANG H B, et al.Strike-slip fault system in Halahatang area of Tarim Basin and its control on reservoirs of Ordovian marine carbonate rock[J].Natural Gas Geoscience,2015,26(S1):80-87.

4
刘军,任丽丹,李宗杰,等.塔里木盆地顺南地区深层碳酸盐岩断裂和裂缝地震识别与评价[J].石油与天然气地质,2017,38(4):703-710.

LIU J, REN L D, LI Z J, et al.Seismic identification and evaluation of deep carbonate faults and fractures in Shunnan area,Tarim Basin[J]. Oil & Gas Geology, 2017, 38(4): 703-710.

5
马乃拜,金圣林,杨瑞召,等.塔里木盆地顺北地区断溶体地震反射特征与识别[J].石油地球物理勘探,2019,54(2):398-403.

MA N B,JIN S L,YANG R Z, et al. Seismic response characteristics and identification of fault-karst reservoir in Shunbei area,Tarim Basin[J].Oil Geophysical Prospecting,2019,54(2):398-403.

6
刘宝增,漆立新,李宗杰,等.顺北地区超深层断溶体储层空间雕刻及量化描述技术[J].石油学报,2020,41(4):412-420.

LIU B Z,QI L X,LI Z J, et al. Spatial characterization and quantitative description technology for ultra-deep fault-karst reservoirs in the Shunbei area[J]. Acta Petrolei Sinica, 2020,41(4):412-420.

7
李宗杰,杨子川,李海英,等.顺北沙漠区超深断溶体油气藏三维地震勘探关键技术[J].石油物探,2020,59(2):283-294.

LI Z J,YANG Z C,LI H Y, et al, Three-dimensional seismic exploration method for ultra-deep fault-related dissolution reservoirs in the Shunbei desert area[J].Geophysical Prospecting for Petroleum, 2020,59(2):283-294.

8
朝鲁门,史继伟,王天柱,等.顺北5井产量递减原因分析及酸化压裂效果评价[J].化工管理,2018(18):143-144.

CHAO L M, SHI J W, WANG T Z, et al. Cause analysis of production decline and evaluation of acid fracturing effect in Shunbei 5 well[J].Chemical Management,2018(18):143-144.

9
曹自成,路清华,顾忆,等. 塔里木盆地顺北油气田1号和5号断裂带奥陶系油气藏特征[J]. 石油与天然气地质, 2020, 41(5): 975-984.

CAO Z C, LU Q H, GU Y, et al. Characteristics of Ordovician reservoirs in Shunbei 1 and 5 fault zones, Tarim Basin[J]. Oil & Gas Geology, 2020, 41(5): 975-984.

10
BULAND A,OMRE H. Bayesian linearized AVO inversion [J]. Geophysics,2003, 68(1):185-198.

11
SOMOZA A V,WATERS K,KEMPER M.Improved Seismic Inversion and Facies Using Regional Rock Physics Trends: Case Study from Central North Sea[C]. Third EAGE Workshop on Rock Physics15-18 November 2015, Istanbul, Turkey,2015.

12
WATERS K D, SOMOZA A V, BYERLEY G,et al. Detecting Production Effects and By-passed Pay from 3D Seismic Data Using a Facies Based Bayesian Seismic Inversion,Forties Field,UKCS[C].78th EAGE Confe-rence & Exhibition 2016,Vienna,Austria,30 May-2 June 2016.

13
KEMPER M,GUNNING J.Joint impedance and facies inversion-Seismic inversion redefined[J].First Break,2014,32(9):89-95.

14
SAMS M, SAUSSUS D. Rock Physics Constraints on Seismic Inversion[C]. EAGE Seismic Driven Reservoir Characterization and Production Management Symposium Kuala Lumpur, Malaysia, 2015:20-22.

15
MCLACHLAN G, PEEL D. Finite Mixture Models[M]. New York: John Wiley & Sons,2000:47-62.

16
SAMS M,WESTLAKE S, THORP J,ZADEH E.Willem 3D: Reprocessed, inverted, revitalized[J].The Leading Edge,2016,35(1):936-940.

17
AVSETH P, MUKERJI T, MAVKO G.定量地震解释[M].李来林,译.北京:石油工业出版社,2009:63-72

AVSETH P, MUKERJI T, MAVKO G.Quantitative Seismic Interpretation[M].LI L L,translated. Beijing: Petroleum Industry Press,2009:63-72.

18
FATTI J L, SMITH G C, VAIL P J, et al. Detection of gas in sandstone reservoirs using AVO analysis:A 3-D seismic case history using the Geostack technique[J].Geophysics,1994,59(9):1362-1376.

19
肖阳,邬光辉,雷永良.走滑断裂带贯穿过程与发育模式的物理模拟[J].石油勘探与开发,2017,44(3): 340-348.

XIAO Y,WU G H,LEI Y L.Analogue modeling of through-go-ing process and development pattern of strike-slip fault zone[J].Petroleum Exploration and Development,2017,44(3):340-348.

Outlines

/