石油物探  2021, Vol. 60 Issue (5): 805-815  DOI: 10.3969/j.issn.1000-1441.2021.05.011
0
文章快速检索     高级检索

引用本文 

郑笑雪, 唐金良. 基于多属性的断溶体特征提取与分带性自动识别[J]. 石油物探, 2021, 60(5): 805-815. DOI: 10.3969/j.issn.1000-1441.2021.05.011.
ZHENG Xiaoxue, TANG Jinliang. Feature extraction and automatic zonation of broken solutions based on multiple attributes[J]. Geophysical Prospecting for Petroleum, 2021, 60(5): 805-815. DOI: 10.3969/j.issn.1000-1441.2021.05.011.

基金项目

国家重大专项(2016ZX05014-001)、中国石化油气重大项目(P17049)和中国石化科技部项目(P18040)共同资助

第一作者简介

郑笑雪(1992—), 女, 硕士, 工程师, 主要从事反演与储层预测的研究工作。Email: zhengxx.swty@sinopec.com

文章历史

收稿日期:2020-09-11
改回日期:2021-07-08
基于多属性的断溶体特征提取与分带性自动识别
郑笑雪, 唐金良    
中国石油化工股份有限公司石油物探技术研究院, 江苏南京 211103
摘要:断溶体作为顺北油气田的主要储集体类型, 具有储层非均质性强, 地震波场特征复杂的特点, 如何精细描述断溶体空间形态和内部结构是断溶体油藏勘探开发的关键。基于二维经验模式分解(BEMD)和纹理分析方法开展了断溶体描述方法研究, 首先采用BEMD信号分解技术对地震数据进行预处理, 通过压制地层结构突出断溶体地震反射特征; 然后利用优化梯度结构张量属性和纹理属性算法, 提高断溶体空间形态检测精度; 最后基于多种优势属性聚类分析, 实现断溶体横向分带性的自动划分, 形成了满足实际生产应用需求的断溶体描述方法技术流程。将该方法技术应用于顺北油气田断溶体识别和描述, 得到的断溶体形态和分带性的认识与井资料吻合度较高。
关键词断溶体    属性    纹理    二维经验模式分解    聚类    降维    强反射    边界刻画    
Feature extraction and automatic zonation of broken solutions based on multiple attributes
ZHENG Xiaoxue, TANG Jinliang    
Sinopec Geophysical Research Institute, Nanjing 211103, China
Abstract: In the oil and gas fields in Shunbei, reservoir type identification methods often return broken solutions, which are indicative of strong reservoir heterogeneity and complex seismic wavefield characteristics.The description of the spatial morphology and internal structure of these reservoirs is essential for their exploration and development.To this aim, bidimensional empirical mode decomposition (BEMD) and texture analysis were used in this study.In particular, BEMD was used to preprocess the seismic data and highlight the seismic reflection characteristics of the broken solution by suppressing the formation structure.Then, the detection accuracy of the spatial morphology of the solution was improved by optimizing the gradient structure tensor attribute and the texture attribute algorithm.Finally, by clustering a number of advantageous attributes, a horizontal zonation of the broken solution was automatically achieved.The proposed method was applied for the identification and description of broken-solution reservoirs in the study area, providing results that were in good agreement with well data.
Keywords: broken solution    property    texture    bidimensional empirical mode decomposition    clustering    dimension reduction    strong reflection    boundary characterization    

鲁新便等[1]在塔河油田2013年的勘探开发中提出了“断溶体”的概念。致密碳酸盐岩或低孔、低渗的碳酸盐岩地层受多期次构造挤压作用, 沿断裂带形成了不同规模的破碎带, 多期岩溶水沿断裂下渗致使破碎带内断裂、裂缝被溶蚀改造, 从而形成板状溶蚀孔、洞储集体, 该储集体在上覆泥灰岩、泥岩等盖层封堵以及侧向致密灰岩遮挡下构成的圈闭类型, 称为“断溶体圈闭”。这类圈闭后期因油气沿断裂运移、充注后成藏, 形成“断溶体油(气)藏”。目前已探明的多个顺北油气富集区均沿大型断裂带分布, 断裂破碎带既是油气疏导通道, 又是成藏的有利空间。顺北特深断溶体埋深普遍超7000m, 构造应力差异导致断裂发育样式多, 岩溶作用导致储层非均质性强, 因其地震波场特征复杂, 故难以通过简单的地球物理方法精细描述和刻画断溶体储层特征, 因而制约了断溶体油藏高效勘探开发。因此, 有必要对断溶体进行针对性的地球物理方法技术和流程研究。

杨德彬[2]、焦方正[3]对顺北地区特深断溶体的油气富集规律进行了综合分析, 所总结的地震反射特征类型及规律对断溶体油藏后续处理解释具有一定的借鉴意义; 马乃拜等[4]研究了断溶体的地震反射特征, 建立的模型有效指导了断溶体的识别及开发; 徐红霞等[5]利用多属性融合技术预测断溶体的连通性; 李鹏飞等[6]利用Geoeast解释软件对断溶体进行了描述与刻画。上述技术均取得了良好的应用效果。

顺北奥陶系上部地层的强地震反射同相轴削弱了下部断溶体的地震响应特征, 为突出断溶体的地震异常特征, 使其便于识别和描述, 需采用相应的信号处理技术消除或减弱反映地层结构特征的连续强反射同相轴[7]。1998年, HUANG等[8]提出了经验模态分解(empirical mode decomposition, EMD)算法, 将信号分解为多个不同频率特性的分量和一个残余分量; 张义平等[9]将其用于地震信号处理, 陈东方等[10]利用经验模态分解算法压制信号中的噪声干扰; 2003年, 二维经验模式分解(bidimensional empirical mode decomposition, BEMD)[11-13]概念被提出并广泛应用于各领域; DAMERVAL等[14]采用三角剖分和分段三次多项式插值改进了EMD算法; 宋平舰等[15]采用BEMD算法对遥感图像进行了处理, 并取得了一定效果; 2005年, 刘忠轩等[16]将经验模态分解算法与纹理分析相结合; 姜晓宇等[17]利用基于相干技术与广义S变换频谱分解的不连续性检测技术对断溶体油藏进行了检测。采用BEMD信号处理技术和分频滤波等相关处理手段对断溶体地震响应特征进行增强预处理, 是实现断溶体精细描述和刻画的前提。

利用地震信息识别和描述断溶体是顺北断溶体油藏勘探和开发的关键环节之一。地震属性技术广泛应用于油气藏表征与储层预测等领域, 其中纹理属性和梯度结构张量属性在计算机视觉和图像分类等领域已经取得了较好的应用效果, 为将其应用领域拓展至地震数据处理, 许多学者进行了研究。GAO[18-20]通过计算纹理属性实现了断层识别; KNUTSSON等[21]提取了地震剖面的结构张量属性; RAO等[22]、FARAKLIOTI等[23]、RANDEN等[24]提取了各向异性介质中结构张量属性并进行了研究; 张军华等[25]和崔世凌等[26]分别将纹理属性和结构张量属性提取算法应用于不同工区完成了断层、河道等特殊异常体的识别; 龚屹等[27]利用纹理属性的聚类分析算法进行裂缝预测; 赵淑琴等[28]将纹理属性提取技术应用于火山岩边界刻画。

综合利用多种属性进行聚类分析, 得到构造、岩性和储层等多种信息, 可有效降低利用属性进行储层识别的多解性。SELIM等[29]发展了K-means方法并得到了广泛应用; MATOS等[30]利用属性进行聚类分析, 并划分了地震相; 杨培杰等[31]利用模糊C均值算法进行地震属性聚类分析; 蔡涵鹏等[32]利用纹理特征进行了半监督地震相分析; 王治国等[33]对地震多属性进行了优化分析; 聚类算法在石油工业中展现了广阔的应用前景, 具有重要的理论和实用价值[34]

本文提出的基于多属性的断溶体特征提取与分带性自动识别方法, 首先以BEMD信号分解理论为基础, 削弱反映地层结构的强反射同相轴, 突出断溶体异常响应; 然后利用优化梯度结构张量和地震纹理属性算法, 提高断溶体储层地震属性识别精度; 再利用优选的多种断溶体识别地震属性, 采用K-means算法进行降维化聚类; 最终实现断溶体空间形态和分带性的自动刻画。

1 断溶体边界刻画技术 1.1 基于BEMD信号处理的断溶体特征增强分析

采用BEMD[11-13]算法对断溶体剖面数据进行预处理, 假定地震数据满足相应条件[11-13], 具体分析过程如下。

1) 定义初始剖面参数, 将叠后地震剖面数据进行灰度化处理后, 作为输入数据:

$ r_{1,1}(x, y)=f(x, y) $ (1)

式中: f(x, y)为地震剖面; r1, 1(x, y)为初始输入剖面。

2) 计算输入剖面rp, k(x, y)的极大值和极小值, 其中p代表从地震剖面计算得到的第p个分解分量, k为筛分次数(P=1, 2, …, p; K=1, 2, …, k)。

3) 分别求取极大值点和极小值点的包络曲面ep, kmax(x, y)和ep, kmin(x, y), 计算包络曲面的均值:

$ m_{p, k}=\frac{e_{p, k}^{\max }(x, y)+e_{p, k}^{\min }(x, y)}{2} $ (2)

4) 将断溶体剖面rp, k(x, y)减去计算的曲面均值mp, k(x, y), 得到:

$ h_{p, k}(x, y)=r_{p, k}(x, y)-m_{p, k}(x, y) $ (3)

5) 重复上述步骤, 当hp, k(x, y)不能再被分解时, 即满足约束条件SD:

$ S_{D}=\sum\limits_{x=0}^{m} \sum\limits_{y=0}^{n} \frac{\left|r_{p, k}(x, y)-r_{p, k+1}(x, y)\right|^{2}}{r_{p, k}^{2}(x, y)} $ (4)

分解得到的第p个IMF分量为Cp(x, y), Cp(x, y)=hp, k(x, y)。

6) 计算残余量:

$ r_{p+1, k}(x, y)=r_{p, 1}(x, y)-C_{p}(x, y) $ (5)

7) 经过p次分解计算得到残余量rp(x, y)=rp+1, k(x, y), 目标信号为:

$ f(x, y)=\sum\limits_{p=1}^{L} C_{p}(x, y)+r_{p}(x, y) $ (6)

对如图 1a所示的地质模型进行BEMD分解测试。该模型通过在断裂带内添加不同大小和发育密度的缝洞体来模拟断溶体的分带性, 其中断溶体核部为大尺度缝洞体且缝洞体发育密集, 翼部缝洞体尺度相对较小且发育密度低, 断溶体上覆地层为泥岩层。从地震剖面可以看出, 断溶体发育的碳酸盐岩地层上部存在强能量地震同相轴, 影响了断溶体的特征识别, 此外, 断溶体弱破碎带(相当于断溶体翼部)地震响应呈现弱反射特征, 当破碎程度较强(相当于断溶体核部)时, 地震响应表现为杂乱异常反射特征(图 1b)。BEMD分解结果(图 1c)在一定程度上使异常更聚焦, 同时对背景数据进行了压制。对模拟地震数据进行处理, 通过BEMD分解得到单频或多频段内的不同组合剖面, 分解得到的剖面分别代表IMF1, IMF2, IMF3分量。分析其有效频段数据组合方式获取有效的BEMD提取结果(图 2), IMF1分量主要反映断溶体空间异常特征, IMF2分量主要反映地层结构信息, IMF3分量代表算法分解的残差部分, IMF1分量增强了断溶体地震响应特征信息, 为实际工区效果分析奠定了基础。

图 1 模型的BEMD分解 a地质模型; b地震剖面; c BEMD分解结果
图 2 模型BEMD提取结果 a地震剖面; b IMF1分量; c IMF2分量; d IMF3分量
1.2 基于地震纹理分析的断溶体识别研究

对剖面进行BEMD处理后, 通过计算纹理属性和结构张量属性, 探讨地震纹理分析在断溶体识别和预测中的应用。

1.2.1 基于地震纹理属性的断溶体识别

HARALICK等[34]利用灰度共生矩阵计算纹理属性, 并将其应用于地球物理领域。首先对地震剖面振幅数据进行灰度化计算, 然后统计一定距离和方向上地震数据的相似性, 得到地震振幅数据在不同方向、不同角度等方面的变化特征。进行灰度化计算时, 沿着某一方向计算距离为δ的任意两个剖面振幅满足一定条件的概率为:

$ \begin{gathered} p(i, j)=\sum\{g(x, y)=i, g(x+\Delta x, y+\Delta y)=j\} \\ i, j=0,1, \cdots, L-1, L \end{gathered} $ (7)

式中: L为地震剖面图像的灰度级数; x, y是地震剖面图像的像素坐标; g(x, y)是(x, y)处的灰度值; δ为数据在x, y方向的增量, $\delta =\sqrt{\Delta {{x}^{2}}+\Delta {{y}^{2}}} $p(i, j)经归一化处理后, 可表示为灰度共生矩阵形式, 灰度共生矩阵反映了剖面信息, 我们为了直观地获得剖面灰度, 即数据的变化情况, 计算由灰度共生矩阵衍生出的多种特征值。

1) 能量。

能量统计了纹理灰度的变化程度, 反映地震剖面振幅的分布均匀程度。

$ E_{\mathrm{ne}}=\sum\limits_{i=0}^{L-1} \sum\limits_{j=0}^{L-1} p_{i, j}^{2} $ (8)

2) 熵。

熵可以度量剖面中纹理的随机性。当共生矩阵的值大小都相等时, 熵为最大值; 反之, 其值较小。

$ E_{\mathrm{n}}=\sum\limits_{i=0}^{L-1} \sum\limits_{j=0}^{L-1} p_{i, j} \log p_{i, j} $ (9)

3) 对比度。

对比度为灰度共生矩阵主对角线附近的惯性矩, 可以反映剖面矩阵的清晰度和纹理的沟纹深浅。

$ C_{\mathrm{on}}=\sum\limits_{i=0}^{L-1} \sum\limits_{j=0}^{L-1}(i-j)^{2} p_{i, j} $ (10)

选取如图 1a所示模型进行纹理属性提取分析, 为了更好地识别断溶体发育特征, 提取不同方向、不同灰度级、不同窗口的纹理属性, 结果如图 3所示。其中对比度属性、熵属性、能量属性均可表征剖面上的断溶体, 优化后的纹理属性连续性得到了增强, 其中熵属性更为连贯, 说明断溶体灰度值变化较大。

图 3 纹理属性提取结果 a某商业软件的对比度属性; b某商业软件的能量属性; c某商业软件的熵属性; d采用本文方法得到的对比度属性; e采用本文方法得到的能量属性; f采用本文方法得到的熵属性
1.2.2 基于地震梯度结构张量属性的断溶体识别研究

梯度结构张量(GST)主要用于计算地震构造属性和预测裂缝发育区, 本文在现有算法的基础上研发优化了GST属性的提取方法, 具体包括以下步骤。

1) 在计算梯度结构张量之前利用BEMD进行断溶体特征增强处理, 然后采用特定滤波器算子对剖面进行处理, 增强断溶体边界。

2) 计算三维地震数据体每一点的梯度矢量:

$ \boldsymbol{g}=\nabla u(x, y, z)=\left[\begin{array}{l} \frac{\partial u(x, y, z)}{\partial x} \\ \frac{\partial u(x, y, z)}{\partial y} \\ \frac{\partial u(x, y, z)}{\partial z} \end{array}\right]=\left[\begin{array}{l} g_{1} \\ g_{2} \\ g_{3} \end{array}\right] $ (11)
$ g_{i}=u(x) \times \frac{\partial}{\partial x_{i}} G\left(x, \sigma_{g}\right), i \in\{1,2,3\} $ (12)

式中: u(x, y, z)为地震数据。

3) 构建梯度结构张量T[33]:

$ \boldsymbol{T}=\boldsymbol{g} \cdot \boldsymbol{g}^{\mathrm{T}}=\left[\begin{array}{c} g_{1}^{2} & g_{1} g_{2} & g_{1} g_{3} \\ g_{2} g_{1} & g_{2}^{2} & g_{2} g_{3} \\ g_{3} g_{1} & g_{3} g_{2} & g_{3}^{2} \end{array}\right] $ (13)

4) 采用不同的高斯核参数对梯度结构张量的每一个成分进行平滑测试, 选取最优结果:

$ \bar{T}_{i j}=T_{i j} \otimes G\left(x, \sigma_{\mathrm{T}}\right) $ (14)

式中: G(x, σT)是高斯核函数; σT是尺度参数。

$ \boldsymbol{T}=\boldsymbol{g} \cdot \boldsymbol{g}^{\mathrm{T}}=\left[\begin{array}{ccc} \bar{g}_{1}{ }^{2} & \bar{g}_{1} \bar{g}_{2} & \bar{g}_{1} \bar{g}_{3} \\ \bar{g}_{2} \bar{g}_{1} & \bar{g}_{2}{ }^{2} & \bar{g}_{2} \bar{g}_{3} \\ \bar{g}_{3} \bar{g}_{1} & \bar{g}_{3} \bar{g}_{2} & \bar{g}_{3}{ }^{2} \end{array}\right] $ (15)

5) 计算梯度结构张量的特征值和特征向量:

$ |\boldsymbol{T} \boldsymbol{\upsilon}-\lambda \boldsymbol{\upsilon}|=0 $ (16)

式中: λ为高斯核函数的相关参数。利用不同高斯核函数的相关参数对模型(图 1a)提取梯度结构张量属性, 结果如图 4所示, 从左至右所选用的参数λ分别为0.1, 0.5, 2.0, 2.5。不同高斯核参数下得到的结构张量属性λ1, λ2差异明显, λ3差异较小; λ1提取结果相对较差, 整体信息虽然较丰富, 但是背景噪声较严重; λ2相对较好, 但是仍存在部分模糊现象, 整体提取效果较理想且不受参数影响。对比剖面可知, 当参数λ为0.1和2.5时, 属性计算效果存在模糊性(图 4a); 参数λ为0.5时, 背景存在较多噪声; 参数λ为2时的效果较为理想。

图 4 模型梯度结构张量属性提取结果 a λ1分量; b λ2分量; c λ3分量
1.3 基于多属性聚类的断溶体分带性自动识别研究

属性解释的多解性造成储层的准确识别存在困难。属性提取结果是同一地质体的不同方面的响应, 地震属性聚类分析是对相互独立参数进行分析, 因此可以在无井约束的情况下完成储层分类预测。本文利用K-means原理将多个地震数据进行聚类, 在实现属性降维的同时, 实现断溶体分带性自动识别。本文将属性集分为c个类别, 算法描述如下。

1) 适当选择c个类的中心作为初始剖面聚类中心。

2) 在第k次迭代中, 对任意一个属性数据, 求其到c个中心剖面的距离, 将该样本归到距离最短的中心所在的类:

$ c^{i}=\arg \min \limits_{j}\left\|x^{i}-\mu_{j}\right\|^{2} $ (17)

式中: xi代表属性剖面; μj代表聚类中心剖面。

3) 利用均值等方法更新该类的中心值:

$ \mu_{j}=\frac{\sum\limits_{i=1}^{m} 1\left\{c^{i}=j\right\} x^{i}}{\sum\limits_{i=1}^{m} 1\left\{c^{i}=j\right\}} $ (18)

4) 对于所有的c个聚类中心, 如果利用步骤2)和3)迭代更新后, 中心值保持不变, 则迭代结束, 否则继续迭代。

采用上述方法计算相关属性, 模型优选的多种属性如图 5a所示, 从左到右分别是为能量性、熵属性、λ3属性、λ属性和对比度属性。聚类结果如图 5b所示。利用该方法可以实现断溶体的分带性预测, 聚类结果证实了方法的有效性。将属性聚类结果和原始剖面叠合显示(图 5c), 聚类结果与所设计模型的分带性高度吻合, 证明了算法的合理性。

图 5 模型属性聚类结果 a模型优选多种属性; b属性分带性聚类结果; c聚类结果与原剖面叠合显示

以断溶体地震响应特征增强预处理方法为基础, 开展以纹理属性和GST属性为核心的多属性断溶体识别分析, 结合空间描述和分带性识别技术, 提出了针对断溶体的一体化处理识别方法, 具体流程见图 6

图 6 基于多属性的断溶体特征提取与分带性自动识别算法流程
2 实际资料处理

对顺北油田B工区三维地震资料进行处理, 结果如图 7所示, 顺北油田发育一系列沿断裂分布的碳酸盐岩断溶体油藏, 其中典型断溶体地震响应特征如图 7中红框所示。对实际地震资料进行BEMD预处理, 结果如图 8所示。对预处理后的实际地震资料提取纹理属性(图 9a)和梯度结构张量属性(图 9b), 可见利用优化算法的纹理属性所识别断溶体形态更加连续, 梯度结构张量属性在一定程度上可以指示断溶体边界轮廓。将不同过井曲线的纹理属性、结构张量属性等作为输入数据进行断溶体分带性聚类分析, 结果如图 10所示, 断溶体分带性与钻井所揭示的储层信息特征相吻合。图 11为多属性聚类断溶体分带性自动识别的时间切片和三维可视化结果, 断溶体分带性刻画精度高, 与井资料吻合程度高, 证明了方法的有效性。

图 7 实际地震资料剖面(含断溶体油藏)
图 8 实际地震资料BEMD预处理结果 a IMF1分量; b IMF2分量; c IMF3分量
图 9 对预处理后的实际地震资料提取的纹理属性(a)和梯度结构张量属性(b)
图 10 断溶体分带性聚类分析结果 a过井1-5曲线; b过井1-1曲线
图 11 多属性聚类断溶体分带性聚类结果时间切片(a)和三维可视化结果(b)
3 结论

1) 断溶体的检测很大程度上依赖于地震剖面处理结果的质量。基于BEMD的地震信号分解技术有效削弱了与地层结构相关的地震强反射同相轴, 增强了断溶体的地震响应特征。

2) 通过优化纹理属性和GST属性算法, 有效强化了纹理属性识别断溶体的能力, 提高了属性识别断溶体的精度, 获得了突出的断溶体识别应用效果。综合利用多种地震属性, 充分挖掘单一属性内涵信息, 利用K值聚类算法可以实现属性的降维, 降低多属性预测储层的多解性, 进一步提高储层预测精度。

3) 地震属性处理结果, 关系到断溶体边界及内幕缝洞储层刻画精度, 随着越来越多智能化优化算法的出现以及混合算法的改进, 智能化自动识别断溶体描述和刻画算法具有广阔的应用前景, 将其应用于断溶体边界识别并结合内幕储层特征研究是下一步研究方向。

参考文献
[1]
鲁新便, 胡文革, 汪彦, 等. 塔河地区碳酸盐岩断溶体油藏特征与开发实践[J]. 石油与天然气地质, 2015, 36(3): 347-355.
LU X B, HU W G, WANG Y, et al. Characteristics and development practice of fault-karst carbonate reservoirs in Tahe area, Tarim Basin[J]. Oil & Gas Geology, 2015, 36(3): 347-355.
[2]
杨德彬. 塔河油田断溶体油藏地震反射特征类型及规律[J]. 石化技术, 2018, 25(4): 128-129.
YANG D B. Type and regularity of the seismic reflection characteristic of fault-karst carbonate reservoirs in the Tahe oilfield[J]. Petrochemical Technology, 2018, 25(4): 128-129. DOI:10.3969/j.issn.1006-0235.2018.04.103
[3]
焦方正. 塔里木盆地顺北特深碳酸盐岩断溶体油气藏发现意义与前景[J]. 石油与天然气地质, 2018, 39(2): 207-216.
JIAO F Z. Significance and prospect of ultra-deep carbonate fault-karst reservoirs in Shunbei Area, Tarim Basin[J]. Oil & Gas Geology, 2018, 39(2): 207-216.
[4]
马乃拜, 金圣林, 杨瑞召, 等. 塔里木盆地顺北地区断溶体地震反射特征与识别[J]. 石油地球物理勘探, 2019, 54(2): 398-403.
MA N B, JIN S L, YANG R Z, et al. Cross-well sand-body identification in the oilfield development late stage[J]. Oil Geophysical Prospecting, 2019, 54(2): 398-403.
[5]
徐红霞, 沈春光, 李斌, 等. 多属性分析技术在碳酸盐岩断溶体预测中的应用[J]. 石油地球物理勘探, 2017, 52(S2): 158-163.
XU H X, SHEN C G, LI B, et al. Fault-karst carbonate reservoir prediction with comprehensive multi-attribute analysis[J]. Oil Geophysical Prospecting, 2017, 52(S2): 158-163.
[6]
李鹏飞, 崔德育, 田浩男. 塔里木盆地塔北地区X区块断溶体刻画方法与效果[J]. 石油地球物理勘探, 2017, 52(S1): 189-194.
LI P F, CUI D Y, TIAN H N. Fault-karst carbonate reservoir description with Geoeast interpretation sub-system in the Tabei Area, Tarim Basin[J]. Oil Geophysical Prospecting, 2017, 52(S1): 189-194.
[7]
王震, 文欢, 邓光校, 等. 塔河油田碳酸盐岩断溶体刻画技术研究与应用[J]. 石油物探, 2019, 58(1): 153-158.
WANG Z, WEN H, DENG G X, et al. Fault-karst characterization technology in the Tahe Oilfield, China[J]. Geophysical Pospecting for Petroleum, 2019, 58(1): 153-158.
[8]
HUANG N E, SHEN Z, LONG S R, et al. The empirical mode decomposition and the Hilbert spectrum for nonlinear and nonstationary time series analysis[J]. Proceedings Mathematilal Physical & Engineering Science, 1998, 454(1971): 903-995.
[9]
张义平, 李夕兵, 赵国彦. 基于HHT方法的爆破地震信号分析[J]. 工程爆破, 2005, 11(1): 1-7.
ZHANG Y P, LI X B, ZHAO G Y. Analysis of blasting vibration signal based on HHT method[J]. Engineering Blasting, 2005, 11(1): 1-7. DOI:10.3969/j.issn.1006-7051.2005.01.001
[10]
陈东方, 吴先良. 采用EMD方法消除瞬态散射回波中的高斯白噪声干扰[J]. 电子学报, 2004, 32(3): 496-498.
CHEN D F, WU X L. Recovery of signal from transient scattered response contaminated by gaussian white noise based on EMD method[J]. Acta Electronica Sinica, 2004, 32(3): 496-498. DOI:10.3321/j.issn:0372-2112.2004.03.033
[11]
NUNES J C, BOUAOUNE Y, DELECHELLE E, et al. Image analysis by bidimensional empirical mode decomposition[J]. Image & Vision Computing, 2003, 21(12): 1019-1026.
[12]
NUNES J C, GUYOT S. Texture analysis based on local analysis of the bidimensional empirical mode decomposition[J]. Machine Vision and Applications, 2005, 16(3): 177-188. DOI:10.1007/s00138-004-0170-5
[13]
NUNES J C, BOUAOUNE Y, DELECHELLE E. Bidimensional empirical mode decomposition modified for texture analysis[J]. SCIA 2003 13th Scandinavian Conference on Image Analysis, 2003, 171-177.
[14]
DAMERVAL C, MEIGNEN S, PERRIER V. A fast algorithm for bidimensional EMD[J]. IEEE Signal Processing Letters, 2005, 12(10): 701-704. DOI:10.1109/LSP.2005.855548
[15]
宋平舰, 张杰. 二维经验模分解在海洋遥感图像信息分离中的应用[J]. 高技术通讯, 2001, 11(9): 62-67.
SONG P J, ZHANG J. The application of two-dimensional EMD to separating contents of oceanic remote sensing images[J]. High Technology Letters, 2001, 11(9): 62-67. DOI:10.3321/j.issn:1002-0470.2001.09.014
[16]
刘忠轩, 彭思龙. 方向EMD分解与其在纹理分割中的应用[J]. 中国科学E辑: 信息科学, 2005, 35(2): 113-123.
LIU Z X, PENG S L. Directional EMD decomposition and its application in texture segmentation[J]. Science in China series E: Information Science, 2005, 35(2): 113-123.
[17]
姜晓宇, 宋涛, 杜文辉, 等. 利用广义S变换频谱分解不连续性检测技术预测断溶体油藏[J]. 石油地球物理勘探, 2019, 54(6): 1324-1328.
JIANG X Y, SONG T, DU W H, et al. Fault-karst carbonate reservoir prediction with the spectral decomposition discontinuity detection based on generalized S transform[J]. Oil Geophysical Prospecting, 2019, 54(6): 1324-1328.
[18]
GAO D L. Structure-oriented texture model regression: Application to seismic structure visualization and interpretation[J]. Expanded Abstracts of 76th Annual SEG Internat Mtg, 2006, 1083-1086.
[19]
GAO D L. Application of seismic texture model regression to seismic facies characterization and interpretation[J]. The Leading Edge, 2008, 27(3): 394-397. DOI:10.1190/1.2896632
[20]
GAO D L. 3D seismic volume visualization and interpretation: An integrated workflow with case studies[J]. Geophysics, 2009, 74(1): W1-W12. DOI:10.1190/1.3002915
[21]
KNUTSSON H, WESTIN C F, ANDERSSON M. Representing local structure using tensors[J]. Proceedings of 6th Scandinavian Conference on Image Analysis, 1989, 619-622.
[22]
RAO A R, SCHUNCK B G. Computing oriented texture fields[J]. Cvgip Graphical Models & Image Processing, 1991, 53(2): 61-68.
[23]
FARAKLIOTI M, PETROU M. The use of structure tensors in the analysis of seismic data[M]. Mathematical methods and modelling in hydrocarbon exploration and production, Heidel berg: Springer, 2005: 978-981.
[24]
RANDEN T, PEDERSEN S I, SNNELAND L. Automatic detection and extraction of faults from three-dimensional seismic data[J]. Expanded Abstracts of 71st Annual SEG Internat Mtg, 2001, 410-412.
[25]
张军华, 董猛, 周振晓, 等. 基于GST的相干体方法研究及应用[J]. 天然气工业, 2007, 27(S1): 387-389.
ZHANG J H, DONG M, ZHOU Z X, et al. Research and application of coherent volume method based on GST[J]. Natural Gas Industry, 2007, 27(S1): 387-389.
[26]
崔世凌, 张军华, 王伟, 等. 地震纹理属性在JJD工区断层识别中的应用[J]. 物探化探计算技术, 2010, 32(3): 304-309.
CUI S L, ZHANG J H, WANG W, et al. The extraction of seismic texture attribute and its application in fault identification[J]. Computing Techniques for Geophysical and Geochemical Exploration, 2010, 32(3): 304-309.
[27]
龚屹, 桂志先, 王鹏, 等. 基于地震纹理属性聚类分析的裂缝分布预测[J]. 科学技术与工程, 2017, 17(30): 167-174.
GONG Y, GUI Z X, WANG P, et al. Prediction of fracture distribution based on clustering analysis of seismic texture attributes[J]. Science Technology and Engineering, 2017, 17(30): 167-174.
[28]
赵淑琴, 顾国忠, 韩宏伟, 等. 纹理属性在火山岩储层预测中的应用[J]. 石油地球物理勘探, 2017, 52(S1): 152-155, 168.
ZHAO S Q, GU G Z, HAN H W, et al. Volcanic reservoir prediction with texture attribute[J]. Oil Geophysical Prospecting, 2017, 52(S1): 152-155, 168.
[29]
SELIM S Z, ISMAIL M A. K-Means-type algorithms: A generalized convergence theorem and characterization of local optimality[J]. IEEE Transactions on Pattern Analysis and Machine Intelligence, 1984, 6(1): 81-87.
[30]
MATOS M C D, YENUGU M M, ANGELO S M, et al. Integrated seismic texture segmentation and cluster analysis applied to channel delineation and chert reservoir characterization[J]. Geophysics, 2011, 76(5): P11-P21.
[31]
杨培杰, 印兴耀, 张广智. 模糊C均值地震属性聚类分析[J]. 石油地球物理勘探, 2007, 42(3): 322-324.
YANG P J, YIN X Y, ZHANG G Z. Cluster analysis of seismic attributes by fuzzy C-mean algorithm[J]. Oil Geophysical Prospecting, 2007, 42(3): 322-324.
[32]
蔡涵鹏, 胡浩炀, 吴庆平, 等. 基于叠前地震纹理特征的半监督地震相分析[J]. 石油地球物理勘探, 2020, 55(3): 504-509.
CAI H P, HU H Y, WU Q P, et al. Semi-supervised seismic facies analysis based on prestack seismic texture[J]. Oil Geophysical Prospecting, 2020, 55(3): 504-509.
[33]
王治国, 尹成, 雷小兰, 等. 预测砂岩孔隙度的地震多属性优化模式对比[J]. 石油地球物理勘探, 2011, 46(3): 442-448.
WANG Z G, YIN C, LEI X L, et al. Multi-attribute optimization analysis for sandstone porosity prediction[J]. Oil Geophysical Prospecting, 2011, 46(3): 442-448.
[34]
HARALICK R M, SHANMUGAM K, DINSTEIN I. Textural features for image classification[J]. IEEE Transactions on Systems, Man, and Cybernetics, 1973, 3(6): 610-621.