Application of Double Width Seismic Data to Channel Sand Body Prediction of X Gas Field in the East China Sea
-
摘要: 东海X气田主力储层埋藏深、横向变化大,常规地震资料品质较差、分辨率低,满足不了勘探开发中日益精细的地质需求。采用三船四源斜缆的采集方式获得宽频宽方位的地震资料,其具有高分辨率、高信噪比、高保真度的特征。应用宽频宽方位地震资料高分辨率、各向异性的优势信息,结合叠前同时反演,求取分方位的河道砂体敏感弹性参数反演体,并且把垂直于河道走向的多个方位反演体进行叠加,从而进行河道砂体的精细预测。相比于常规地震资料,基于宽频宽方位地震资料的储层反演提高了河道砂体的预测精度,为东海X气田的滚动勘探与开发开采奠定了基础。Abstract: Due to that the main reservoirs of the X gas field in the East China Sea were deeply buried and have large lateral changes, the conventional seismic data resulted in poor quality and low resolution, which couldn’t meet the increasingly refined geological requirements in exploration and development. Seismic data with wideband and wide azimuth was obtained by using the acquisition method of three ships and four sources with oblique cables, which held the characteristics of high resolution, high signal-to-noise ratio and high fidelity. By taking advantage of the superior information of wideband and wide azimuth seismic data which was high-resolutional and anisotropic, combined with the simultaneous prestack inversion, the inversion body of sensitive elastic parameters of channel sand bodies in different azimuths can be obtained, and superimpose multiple azimuth inversion bodies perpendicular to the direction of the river channel to carry out fine predictions of channel sand bodies. Compared the conventional seismic data, reservoir inversion based on wideband and wide azimuth seismic data improved the prediction accuracy of channel sand bodies, laying a foundation for the progressive exploration and development of the X gas field in the East China Sea.
-
河流相储层是一类重要的含油气储集层,其精细描述的核心是有效识别出河道。但是此类储层纵向上常表现为砂体多期叠置,可识别厚度小于地震分辨率;横向上物性变化快,非均质性特征强,河道边缘特征不清晰。为了提高此类储层的识别精度,众多专家学者进行了深入研究,相干体[1-2]、地震波形结构属性[3]、甜点属性[4]、谱分解和分频属性融合[5]、概率神经网络[6]、相层双控智能识别描述[7]等方法技术得到广泛应用。
上述研究主要集中在河道构型、地震相整体特征、含油气分布预测等相对宏观的信息,对河道边缘特征的精细刻画还较少,在河道信号微弱、强屏蔽干扰或者横向不连续等情况下难以精准刻画储层。而图像边缘检测可通过捕捉图像局部特征的细微突变实现边缘的定位[8],其本质是通过突出图像的有效边缘信息、弱化非边缘信息来提取边缘细节。在地震解释领域,已有学者将此项技术应用于断层与裂缝预测中[9-11]以及砂坝、砂体厚度的描述中[12-14]。然而现有相关文献对河道边缘特征的数学意义以及算子处理的本质含义分析不够深入,还有待精细解剖。
本文从河道边缘几何特征与数学含义出发,从图像处理角度分析相干体边缘特征识别优势及存在的不足,以Sobel算子为例探究算子处理边缘识别的本质。对于图像灰度不均导致小河道识别困难的问题,提出基于直方图均衡化和模糊集理论的图像增强方法。上述技术方法均用三维河道理论模型和实际资料加以检验,应用效果显著,相关认识与结论对开展河流相储层精细描述有一定的借鉴作用。
1. 河道边缘特征的几何表征与三维建模
1.1 河道边缘几何特征与数学含义
河流相储层中的河道,其垂向厚度一般小于调谐厚度,其正演模型的左右两边缘振幅为零,中心位置振幅最大(图1(a))。若提取振幅切片,河道内有幅值,呈“实心”特征。对图1(a)截面计算一阶导数,导数从左到河道中心为正,边缘特征为“右负左正”或者“左正右负”。如果对一阶导数求模,左右两边界均为正值,边缘呈“左右双正”特征(图1(c))。二阶导数是一阶导数的导数,所以它的值与原始振幅相比,符号相反,但河道空间形态会更加锐化(图1(d))。
1.2 河流相储层三维建模
为便于对本文方法进行验证,构建200×300×100(线×道×时间)三维模型。模型有5个水平层,设计了4条河道,速度结构如图2(a)所示,河道1到河道4的砂岩速度分别为3600、3200、3150和3300 m/s,围岩速度为3100 m/s。从速度结构上看,各支河道反射系数由强到弱分别为河道1、河道4、河道2、河道3。河道1和河道2的埋深设计在47~57 ms之间,河道3和河道4在55~60 ms之间,河道1和河道4宽度为10道,河道3宽度为15道,河道2宽度为5道。正演模型采样间隔为1 ms,子波采用主频为30 Hz的雷克子波,加入最大振幅5%的随机噪声(图2(b))。
2. 基于相干属性的经典边缘识别技术
2.1 方法原理概述
对于地质异常体的边缘特征识别,常用的地震属性有均方根振幅、相干、曲率、甜点等,尤以相干属性为最常用[15]。对于以特征值为基础的第3代相干,其表征公式常有以下几种(表1)[2],物理意义最明确、应用最广的是最大特征值相干。通过相干属性计算,河道特征比原始切片会更加清晰。
表 1 特征值相干不同表征公式Table 1. Different characterization formulas of eigenvalue coherence时间 作者 表征公式 物理含义 1999 Gersztenkorn
和Marfurt${C_{31}} = \displaystyle \dfrac{{{\lambda _1}}}{{\sum\limits_{j = 1}^J {{\lambda _j}} }}$ 用最大特征值在所有特征值中的占比来表示相干 2000 Randen等 ${C_{32} } = \displaystyle \dfrac{ {2{\lambda _2} } }{ { {\lambda _1} + {\lambda _3} } } - 1 = \dfrac{ { {\lambda _2} - {\lambda _3} - ( { {\lambda _1} - {\lambda _2} } )} }{ { {\lambda _1} + {\lambda _3} } }$ 被称为“chaos”的相干属性 2002 Bakker ${C_{33} } = \displaystyle \dfrac{ {2{\lambda _2}( { {\lambda _2} - {\lambda _3} })} }{ {( { {\lambda _1} + {\lambda _2} } )( { {\lambda _2} + {\lambda _3} } )} }$ 重点考虑第2特征值和第3特征值的差异 2007 Donias等 ${C_{34} } = 1 - \displaystyle \dfrac{3}{2}\dfrac{ { {\lambda _2} + {\lambda _3} } }{ { {\lambda _1} + {\lambda _2} + {\lambda _3} } } = \dfrac{ { {\lambda _1} - {\lambda _2} + {\lambda _1} - {\lambda _3} } }{ {2( { {\lambda _1} + {\lambda _2} + {\lambda _3} } )} }$ “disorder”相干属性 2017 Wu ${C_{35}} = \displaystyle \dfrac{{{\lambda _1}{{ - }}{\lambda _2}}}{{{\lambda _1}}}$ 利用第1特征值和第2特征值的差异 2.2 模型测试与效果分析
对图2(b)所示模型,计算C31相干,提取时间为55 ms处的相干切片与原始切片进行比较(图3)。可以看到:①原始切片(图3(a))展示的是振幅信息,所以河道内是包含振幅信息的,河道显示是“实心”的;②相干(图3(b))突出的是异常信息,所以在相干体上展示的是河道边缘信息,河道是“空心”的;③无论是原始切片还是相干切片,显示的幅值范围都需要人工设置比较合理的范围,才能做到基本清晰(如图3(c)调整了色标幅值范围,河道边缘对比度有了一定的提高,但受背景噪声影响,即使精细调整色标范围仍难以精准刻画能量较弱的河道边缘)。
2.3 实际资料应用与存在问题分析
选取胜利油田罗家高密度资料作为实例数据,数据范围400线和600道,时间950~1050 ms。研究区曲流河、辫状河等河流相储层比较发育,河道相互叠置,横向变化大,砂体厚度小,河道砂体的精细描述是勘探开发的重点。
通过观察原始振幅切片(图4(a),T0=984 ms)可以看到:①研究区中间偏右侧南北向发育有一条宽度较大的河道;②北部1条河道虽然不是很宽,但振幅较强。研究区大量发育的是河道较小的曲流河,它们的相干值差异不大,大量值分布在白色的小值上(图4(b),显示范围0.45~1.0),此时需要人为设置显示范围才能较好地展示全区河道分布(图4(c),显示范围0.8~1.0)。且即使色标做了调整,对河道的的刻画效果仍然欠佳,个别振幅较小的河道的边缘特征不明显,进行基于图像处理的目标处理具有客观必要性。
3. 基于算子处理的边缘增强技术
3.1 方法基本原理
计算一个点的导数,容易出现异常值。利用多个点特殊的组合,进行数值或导数计算,可以得到具有某些特殊特征的图像,这就是算子处理的功效。在图像处理领域,常用算子包括Canny、Prewitt、Sobel、Roberts、Laplace、LoG算子等,其中Sobel算子在梯度运算中应用了局部平均的思想,使其具备较强的抗噪能力,检测结果不易受噪声干扰,在地震边缘检测中已有应用[12-14]。
下面以Sobel算子为例,阐述其物理含义与计算公式。如图5,对于中心点(x, y),将第3列的和(局部平均)减去第1列的和,构建的是该点x方向的梯度Gx。求和运算可以均衡幅值,起到平滑的作用,计算时距离中心点最近的两个点(x+1, y)和(x − 1, y)给予2倍的权重;相减(差分)操作可以突出异常的边缘信息,同时起到去除相似背景的作用。类似地,可以得到y方向的梯度Gy。上述操作既可以用计算公式(1)实现,也可以核函数卷积实现,核函数设置见图5。
$$ \left\{ {\begin{array}{*{20}{l}} {{\boldsymbol{G}}_x} = & \{u (x + 1,y - 1) + 2u (x + 1,y) + u (x + 1,y + 1)\} -\\ & \{u (x - 1,y - 1) + 2u (x - 1,y) + u (x - 1,y + 1)\} \\ {{\boldsymbol{G}}_y} = &\{u (x - 1,y + 1) + 2u (x,y + 1) + u (x + 1,y + 1)\} - \\ &\{u (x - 1,y - 1) + 2u (x,y - 1) + u (x + 1,y - 1)\} \end{array}} \right. 。 $$ (1) 实际应用时既可以用x、y方向的梯度单独表示,使其具有方向性,也可以用1范数、2范数或
$ \infty $ 范数(极大模)表示。$$ {{\boldsymbol{G}}_1} = \Big| {{{\boldsymbol{G}}_x}} \Big| + \Big| {{{\boldsymbol{G}}_y}} \Big| \text{,} $$ (2) $$ {{\boldsymbol{G}}_2} = \sqrt {{\boldsymbol{G}}_x^2 + {\boldsymbol{G}}_y^2} \text{,} $$ (3) $$ {\boldsymbol{G}}_{\infty} = \max \Big( {\big| {{{\boldsymbol{G}}_x}} \big|,\big| {{{\boldsymbol{G}}_y}} \big|} \Big) 。 $$ (4) 3.2 模型测试与效果分析
为检验本文方法有效性,对图6(a)所示切片,用Sobel算子进行处理,结果如图6(b)~图6(d)所示。梯度处理后河道两侧出现了不同极性的边缘,以能量较强的河道1为例,边缘右负左正、下负上正,这与前文中图1的理论分析一致。经Gx方向处理后,线方向的河道边缘特征得到较好展示;Gy方向处理后,道方向的河道边缘特征有所加强,但南北向的河道空间特征变得模糊。图6(d)是计算Gx和Gy模极大值结果,可以看到由于已计算了梯度的绝对值,河道出现类似于相干检测的边缘特征。总的来看,原始切片在采用Sobel算子做相关处理后相比于传统相干属性其边缘检测效果获得一定的增强,但对于能量较弱的河道边缘检测效果仍有待提高。
3.3 实际资料应用与尚存在的问题
对图7(a)实际资料进行处理,结果如图7(b)~图7(d)所示。对比分析可知,不同方向Sobel算子对南北向和东西向河道识别能力不一样,这与模型讨论认识一致;相对方向梯度,使用模值表征更利于刻画细小河道。
需要指出的是,以上只是对切片进行Sobel卷积操作,它没有构造导向和挖掘数据子体内部特征的功能,效果还比较受限,分辨率还有待提高。
4. 基于图像灰度处理的边缘增强技术
经相干体处理后的切片其数值分布范围很不均匀,且该经典方法与本文后续提出的基于算子处理的边缘增强技术对能量较弱地质异常体边缘显示效果仍然不佳。针对这样的问题,我们采用基于图像灰度处理的边缘增强技术来加以改善。
图像边缘的抽象定义为图像灰度发生空间突变的像素集合,以灰度值表征的图像在均匀区域之间的幅值突变被称为边缘。若将地震数据视为图像,地震振幅视为灰度,地震边缘对应的幅值突变区域往往是波阻抗界面产生的位置,它是地下构造或岩性信息的一种反映。
4.1 方法基本原理
4.1.1 灰度均衡图像处理技术
为解决相干处理后切片数值分布范围不均匀的问题,文章引入图像处理中的灰度直方图均衡化技术。首先统计灰度图像中元素具有的相同灰度值i出现的次数ni,计算灰度值i出现的概率:
$$ P(i) = \frac{{{n_i}}}{n},\begin{array}{*{20}{c}} {}&{} \end{array}i \in 0,1,\cdots,{L_0} - 1 , $$ (5) 式中,n为图像像素总数,L0为图像灰度矩阵中非重复灰度级的个数。再计算灰度值i的累计概率:
$$ c(i) = \sum\limits_{k = 0}^i {P(k)} 。 $$ (6) 按以下映射关系对原始像素进行处理,其中
$ {C_{\text{p}}}(i) $ 为当前元素灰度值的均衡化结果,Le为图像处理预期的总灰度级(取值256):$$ {C_{\text{p}}}(i) = {\rm{round}}\Big( {c(\,i\,) \times ({L_e} - 1)} \Big) \text{,} $$ (7) round为四舍五入取整函数,用
${C_{\text{p}}}(i)$ 替代原始灰度值为i的像素,得到最终均衡结果。4.1.2 基于模糊集理论的边缘增强技术
为了进一步分离出相干等切片上的有效边缘信息,对均衡化后的图像做进一步模糊增强处理。方法核心思想是利用隶属度函数将空间域图像变换为模糊域内的模糊特征平面,在模糊特征平面上对模糊特征进行非线性变换,最后将其变换回空间域,得到增强处理的图像。
在模糊集理论中,将一幅灰度级为L的M×N图像u视为一个模糊集,通过隶属度函数将其由空间域转换到模糊特征域,并由模糊矩阵P进行表征:
$$ {\boldsymbol{P}} = \bigcup\limits_{i = 1}^M {\bigcup\limits_{j = 1}^N {{\mu _i}_j} } /{u_i}_j;\begin{array}{*{20}{c}} {}&{} \end{array}i = 1,2, \cdots ,M;\begin{array}{*{20}{c}} {}&{} \end{array}j = 1,2, \cdots ,N \text{,} $$ (8) 式中,
${\mu _i}_j\left( {0 \leq {\mu _i}_j \leq 1} \right)$ 表示图像中第(i,j)个像素点的灰阶,$ {u_i}_j $ 为特定灰度级的隶属度,一般这个特定灰度级取图像最大灰度$ {u_{\max }} $ ,$ {\mu _i}_j $ 组成模糊特征平面$ \left\{ {{\mu _i}_j} \right\} $ 。采用Pal-King算法中的隶属函数[16]:$$ {\mu _i}_j = F({u_{ij}}) = {\Big( {1 + \big( {( {{u_{\max }} - {u_{ij}}} )/{F_d}} \big)} \Big)^{ - {F_e}}} \text{,} $$ (9) 式中,Fd为倒数模糊因子,Fe为指数模糊因子,这里取Fd=128,Fe=2。对
$ {\mu _i}_j $ 进行模糊增强,非线性增强变换函数为:$$ {\mu _i}_j' = {E_r}({\mu _i}_j) = {E_1}\left( {{E_{r - 1}}({\mu _i}_j)} \right);\begin{array}{*{20}{c}} {}&{} \end{array}r = 1,2, \cdots ,$$ (10) $$ {E_1}({\mu _i}_j) = \left\{ {\begin{aligned} &2{\mu _i}{{_j}^2},&0 \leq {\mu _i}_j \leq {\mu _c} \\ &1 - 2{{\left( {1 - {\mu _i}_j} \right)}^2},&{\mu _c} \leq {\mu _i}_j \leq 1 \end{aligned}} \right. \text{,} $$ (11) 式中,
$ {E_r} $ 表示对函数E的r次迭代运算,$ {\mu _c} $ 为渡越点,一般取0.5。利用$ {F^{ - 1}} $ 对$ {\mu _i}_j' $ 进行逆变换即可得到增强后的图像$ {l_{ij}} $ :$$ {l_{ij}} = {F^{ - 1}}({\mu _i}_j') = {u_{\max }} + {F_d}\left( {1 - {{({\mu _i}_j')}^{ - \frac{1}{{{F_e}}}}}} \right) 。 $$ (12) 4.2 模型测试与效果分析
4.2.1 灰度均衡直方图对比及分析
从前面分析已知,相干切片上能量较弱的河道2与河道3的边缘特征依然不够清晰(图8(a)),视觉分辨率不高。分析其原因主要是灰度值过于集中(图8(c),灰度值大多聚焦在白色的大值上),所以与背景值接近的河道不易识别。图8(b)为灰度均衡化后结果,此时直方图主灰度级已分散到6个以上(图8(d)),由此河道弱边缘特征得到较大增强,不过背景噪声也被同步放大,需要后续进一步处理。
4.2.2 模糊边缘增强及迭代次数评价
针对均衡后噪声同步放大的问题(图9(a)),用模糊增强方法做进一步处理,图9(b)~图9(e)为1~4次模糊增强结果。可以看到:①1次模糊增强处理后,背景噪声得到较大程度的压制,弱边缘得到进一步增强;②模糊次数增加到 4次,噪声确实已得到较好去除,但边缘信息损失也很大。总的来看,1~2次模糊增强效果较好。
4.3 实际资料应用与效果评价
利用以上方法对实际资料相干切片,做直方图均衡化和模糊增强处理,结果如图10所示。可以看到:①原始相干切片的灰度范围较为集中,对比度不高,直方图均衡化后河道边缘特征变得比较清楚,由此也说明这样的图像目标处理十分必要;②与模型测试类似,模糊增强处理1~2次,能消除部分较强能量的背景信息,使河道边缘得到更好展示,有利于储层精细描述与油气预测。
5. 结语
对于频繁分叉改道、多期砂体叠置的复杂河流相储层,直接用振幅切片,识别能力有限。多道相干计算,利用河道的空间信息,河道的边缘特征得到了有效的提取。但由于背景噪声的影响,相干体上能量较弱的细小河道还是得不到很好成像。
利用图像处理中的灰度直方图均衡化技术,可以有效提高相干属性切片的对比度,凸显微小边缘特征。在此基础上进一步结合模糊增强,分离有效边缘和无效背景,可以进一步增强相干属性的边缘表征能力。
图像处理中,Canny、Prewitt、Sobel等算子处理,物理意义明确,操作简单、快捷,可用于河流相储层的边缘检测。
-
-
[1] MOLDOVEANU N, EGAN M. From narrow-azimuth to wide-azimuth acquisition in the Gulf of Mexico[J]. First Break, 2006, 24(12): 1−8.
[2] CORDSEN A, GALBRAITH M. Narrow-versus wide-azimuth land 3D seismic surveys[J]. The Leading Edge, 2002, 21(8): 764−770. doi: 10.1190/1.1503181
[3] 刘依谋, 印兴耀, 张三元, 等. 宽方位地震勘探技术新进展[J]. 石油地球物理勘探, 2014,49(3): 596−610, 420. LIU Y M, YIN X Y, ZHANG S Y, et al. Recent advances in wide-azimuth seismic exploration[J]. Oil Geophysical Prospecting, 2014, 49(3): 596−610, 420. (in Chinese).
[4] 黄福强, 张异彪, 李斌, 等. 高密度宽方位深拖地震采集技术在西湖凹陷的应用[J]. 海洋石油, 2019, 39(4): 1-5. HUANG F Q, ZHANG Y B, LI B, et al. Application of high density and wide azimuth seismic acquisition technology with deep-tow in Xihu Sag[J]. Offshore Oil, 2019, 39(4): 1-5. (in Chinese).
[5] KAPOOR J, MOLDOVEANU N, EGAN M, et al. Subsalt imaging: The RAZ-WAZ experience[J]. The Leading Edge, 2007, 26(11): 1414−1422. doi: 10.1190/1.2805764
[6] 王建花, 张金淼, 吴国忱. 宽方位杨氏模量反演和裂缝预测方法及应用—以渤中凹陷H构造潜山勘探为例[J]. 石油地球物理勘探, 2021,56(3): 593−602, 415. WANG J H, ZHANG J M, WU G C. Wide-azimuth Youngs modulus inversion and fracture prediction: An example of H structure in Bozhong sag[J]. Oil Geophysical Prospecting, 2021, 56(3): 593−602, 415. (in Chinese).
[7] 熊金红, 陈岑, 曹占元, 等. 基于叠前地震全方位各向异性预测裂缝发育−以普光气田须家河组为例[J]. 地质力学学报, 2017,23(2): 280−287. doi: 10.3969/j.issn.1006-6616.2017.02.011 XIONG J H, CHEN C, CAO Z Y, et al. Fracture development prediction based on azimuthal anisotropy analysis of prestack seismic: Take Xujiahe formation of Puguang filed for example[J]. Journal of Geomechanics, 2017, 23(2): 280−287. (in Chinese). doi: 10.3969/j.issn.1006-6616.2017.02.011
[8] 夏亚良, 魏小东, 王中凡, 等. OVT 域方位各向异性技术在中非花岗岩裂缝预测中的应用研究[J]. 石油物探, 2018,57(1): 140−147. doi: 10.3969/j.issn.1000-1441.2018.01.018 XIA Y L, WEI X D, WANG Z F, et al. Application of azimuthally anisotropy by OVT gather for granite fracture prediction in Central Africa[J]. Geophysical Prospecting for Petroleum, 2018, 57(1): 140−147. (in Chinese). doi: 10.3969/j.issn.1000-1441.2018.01.018
[9] 张宪旭. 基于宽方位地震数据的煤层裂隙预测方法研究[J]. 能源与环保, 2021,43(7): 114−119, 136. ZHANG X X. Study on coal seam fracture prediction method based on wide azimuth seismic data[J]. China Energy and Environmental Protection, 2021, 43(7): 114−119, 136. (in Chinese).
[10] 刘晓晖, 涂齐催, 姜雨, 等. 双宽地震资料在西湖凹陷 R 气田裂缝预测中的应用[J]. 海洋地质前沿, 2020,36(8): 57−64. LIU X H, TU Q C, JIANG Y, et al. Application of double-width seismic data to fracture prediction in R gasfield of Xihu sag[J]. Marine geology frontiers, 2020, 36(8): 57−64. (in Chinese).
[11] 王斌, 尹路, 陈永波, 等. 基于 OVT 域资料的低渗透砂砾岩储集层预测[J]. 新疆石油地质, 2018,39(1): 104−108. WANG B, YIN L, CHEN Y B, et al. OVT domain-based prediction technology for low-permeability sandy-conglomerate reservoirs[J]. Xinjiang Petroleum Geology, 2018, 39(1): 104−108. (in Chinese).
[12] 周华建. 基于叠前 OVT 域偏移的河道砂体预测方法[J]. 岩性油气藏, 2019,31(4): 112−120. doi: 10.12108/yxyqc.20190412 ZHOU H J. Prediction method of channel sand body based on prestack migration in OVT domain[J]. Lithologic Reservoirs, 2019, 31(4): 112−120. (in Chinese). doi: 10.12108/yxyqc.20190412
[13] 张春广, 丛培泓, 田军, 等. 宽方位角地震资料解释技术及应用[J]. 天然气工业, 2007,27(S1): 152−154. ZHANG C G, CONG P H, TIAN J, et al. Application of the wide azimuth seismic interpretation technology[J]. Gas Industry, 2007, 27(S1): 152−154. (in Chinese).
[14] 刘传虎. 宽方位地震技术与隐蔽油气藏勘探[J]. 石油物探, 2012,51(2): 138−145, 104. doi: 10.3969/j.issn.1000-1441.2012.02.005 LIU C H. Wide azimuth seismic technique and subtle hydrocarbon reservoir exploration[J]. Geophysical Prospecting for Petroleum, 2012, 51(2): 138−145, 104. (in Chinese). doi: 10.3969/j.issn.1000-1441.2012.02.005
[15] ZOEPPRITZ K. Erdbebenwellen VII b: Über reflexion und durchgang seismischer wellen durch unstetigkeitsflächen[J]. Gottinger Nachr, 1919, (1): 66−84.
[16] FATTI J L, SMITH G C, VAIL P J, et al. Detection of gas in sandstone reservoirs using AVO analysis: A 3D seismic case history using the geostack technique[J]. Geophysics, 1994, 59(9): 1362−1376. doi: 10.1190/1.1443695
[17] 方中于, 朱江梅, 李勇, 等. 叠前地震同步反演方法及应用[J]. 成都理工大学学报(自然科学版), 2014,41(2): 152−156. FANG Z Y, ZHU J M, LI Y, et al. Study on prestack seismic simultaneous inversion and its application[J]. Journal of Chengdu University of technology (Science & Technology Edition), 2014, 41(2): 152−156. (in Chinese).
[18] 张倩, 潘仁芳, 王湘君, 等. 叠前同步反演在港中油田的应用[J]. 长江大学学报(自然科学版), 2014,11(2): 92−95. ZHANG Q, PAN R F, WANG X J, et al. Application of prestack simultaneous inversion in Gangzhong oilfield[J]. Journal of Yangtze University (Science & Technology Edition), 2014, 11(2): 92−95. (in Chinese).
[19] 姜勇, 李宁, 涂齐催. 基于叠前同步反演的储层参数预测技术研究及应用[J]. 油气藏评价与开发, 2016,6(1): 9−13. doi: 10.3969/j.issn.2095-1426.2016.01.003 JIANG Y, LI N, TU Q C. Study on reservoir parameters prediction technique based on prestack simultaneous inversion and its application[J]. Reservoir Evaluation and Development, 2016, 6(1): 9−13. (in Chinese). doi: 10.3969/j.issn.2095-1426.2016.01.003
[20] 秦德文, 姜勇, 侯志强, 等. 叠前同步反演技术在西湖凹陷低孔渗储层“甜点”预测中的应用[J]. 油气藏评价与开发, 2015,5(6): 12−15. doi: 10.3969/j.issn.2095-1426.2015.06.003 QIN D W, JIANG Y, HOU Z Q, et al. Application of prestack synchronous inversion technology in “sweet spot”prediction of low porosity and permeability reservoir in Xihu sag[J]. Reservoir Evaluation and Development, 2015, 5(6): 12−15. (in Chinese). doi: 10.3969/j.issn.2095-1426.2015.06.003
[21] 詹仕凡, 陈茂山, 李磊, 等. OVT 域宽方位叠前地震属性分析方法[J]. 石油地球物理勘探, 2015,50(5): 956−966. ZHAN S F, CHEN M S, LI L, et al. OVT domain wide-azimuth prestack seismic attribute analysis[J]. Oil Geophysical Prospecting, 2015, 50(5): 956−966. (in Chinese).
[22] 李红星. 一种基于OVT域的地震储层预测方法[J]. 特种油气藏, 2019,26(4): 52−57. doi: 10.3969/j.issn.1006-6535.2019.04.009 LI H X. A seismic reservoir prediction based on OVT domain[J]. Special Oil & Gas Reservoirs, 2019, 26(4): 52−57. (in Chinese). doi: 10.3969/j.issn.1006-6535.2019.04.009
-
期刊类型引用(0)
其他类型引用(1)