石油物探  2022, Vol. 61 Issue (2): 339-347, 363  DOI: 10.3969/j.issn.1000-1441.2022.02.015
0
文章快速检索     高级检索

引用本文 

张鹏云, 孙建孟, 邓志文, 等. 滩坝砂储集体测井产能等级划分与地震属性横向预测[J]. 石油物探, 2022, 61(2): 339-347. DOI: 10.3969/j.issn.1000-1441.2022.02.015.
ZHANG Pengyun, SUN Jianmeng, DENG Zhiwen, et al. Logging productivity grading and lateral prediction of seismic attributes in a beach-bar sandstone reservoir[J]. Geophysical Prospecting for Petroleum, 2022, 61(2): 339-347. DOI: 10.3969/j.issn.1000-1441.2022.02.015.

基金项目

国家自然科学基金项目(41874138)和国家重大科技专项(2016ZX05006002-004)共同资助

第一作者简介

张鹏云(1990—), 男, 博士在读, 主要研究方向为测井数据处理与解释。Email: pengchengwanlizpy@126.com

文章历史

收稿日期:2021-02-21
滩坝砂储集体测井产能等级划分与地震属性横向预测
张鹏云1,2, 孙建孟2, 邓志文3, 林磊2, 王军4    
1. 中海油田服务股份有限公司油田技术研究院, 北京 101149;
2. 中国石油大学(华东)地球科学与技术学院, 山东青岛 266580;
3. 中国石油集团东方地球物理勘探有限责任公司, 河北涿州 072750;
4. 中国石油化工股份有限公司胜利油田勘探开发研究院, 山东东营 257099
摘要:滩坝砂储层因其储层物性差、单层厚度薄等特点而导致产能预测较为困难, 提出了一种新的基于数据和模型双驱动的井震结合产能等级划分方法。通过机器学习中的降维算法结合聚类算法, 并利用聚类算法中的肘部法确定最佳聚类数, 进行储层测井产能等级自动划分。首次在平面径向流公式的基础上建立了测井产能等级指示模型, 为地震属性横向预测提供井点刻度。以胜利油田W工区为例, 利用最小二乘法推导出了测井产能等级指数公式, 计算了10口井目的层单井产能等级指数。对提取的多种地震属性在井点处的属性值与测井产能等级分类结果进行Pearson相关性分析, 优选了3种与产能等级指数显著相关且两两之间相互独立的地震属性。利用支持向量回归算法, 建立了地震属性融合的储层产能等级指数平面图, 经3口验证井检验, 产能等级指数平面图与实际产能吻合程度较好。研究结果表明, 基于测井产能等级划分的地震属性横向预测方法可以有效预测靶区储层产能等级。
关键词滩坝砂储层    测井产能等级划分    地震属性    横向预测    
Logging productivity grading and lateral prediction of seismic attributes in a beach-bar sandstone reservoir
ZHANG Pengyun1,2, SUN Jianmeng2, DENG Zhiwen3, LIN Lei2, WANG Jun4    
1. Well-tech R & D Institute, China Oilfield Services Limited, Beijing 101149, China;
2. School of Geosciences, China University of Petroleum, Qingdao 266580, China;
3. BGP Inc., China National Petroleum Corporation, Zhuozhou 072750, China;
4. Research Institute of Exploration and Development of Shengli Oilfield, Sinopec, Dongying 257099, China
Abstract: Predicting the productivity of beach-bar sandstone reservoirs is difficult owing to their poor physical properties and the low thickness of the individual layers.In this study, a new productivity grading method based on well logging and seismic data was proposed.Machine-learning dimensionality reduction and clustering algorithms were used to automatically evaluate the logging productivity grading.The optimal clustering number was determined using the elbow rule.For the first time, by using the plane radial flow formula, a model for the logging productivity grading index was established.This model can be used to provide a good one-point calibration for the lateral prediction of seismic attributes.Taking the W well area of the Shengli Oilfield as an example, a practical formula for the logging productivity grading index was derived using the least squares method.The single-well target layer logging productivity grading indices of 10 wells were calculated.Pearson's correlation analysis between the various attribute values extracted at the well point and the logging productivity classification results was conducted.Finally, three seismic attributes, which are significantly related to the productivity grading index and are independent of each other, were selected.Using the support vector regression algorithm, a plane graph of the reservoir productivity grading index combined with seismic attributes was generated.The plane graph of the reservoir productivity grading index was in good agreement with the actual production capacity according to data from three verification wells.The results of this study showed that the combined lateral prediction of seismic attributes based on the logging productivity grading can be effective in predicting the reservoir productivity grade in the target area.
Keywords: beach-bar sandstone reservoir    logging productivity grading    seismic attributes    lateral prediction    

胜利油田W工区目标储层属于典型的滩坝砂储层, 储层孔隙度基本小于10%, 空气渗透率基本小于1 mD(1 mD≈0.987×10-3 μm2), 单层厚度为1~5 m。随着勘探开发的深入, 这类面积相对较小、横向变化大、单层砂体厚度薄、纵向上多层互层的储集体日益受到重视[1-3]。如何在较为复杂的滩坝砂储层中进行产能预测进而找到“甜点区”(优质储层区域)是此类储层油气勘探中亟需解决的问题[4-6]

经过多年的研究和实践, 已有学者在测井产能预测方面取得了诸多进展, 其基本思想是综合测井和岩心分析等资料, 建立储层分类标准, 在此基础上提取储层产能测井敏感参数, 建立产能预测模型, 实现储层产能的定性或定量预测[7-10]。这些产能分级和预测方法证明利用测井资料进行储层产能预测是可行的, 但很少将地震数据考虑在内, 因此仅能获得井点处的产能信息, 无法进行储层产能的横向预测。测井相划分是指依据测井资料通过一些数学方法将地层划分为多个具有不同地质意义的纵向条带, 通常不同的测井相其岩性与物性具有较大差异[11-12]。地震属性是指由地震数据经过数学变换导出的有关地震波的几何形态、运动学、动力学和统计学特征, 目前在岩性预测、储层含油气性预测等研究中应用较为广泛[13-14]。由于测井曲线的分辨率远高于地震数据, 为了提高滩坝砂这类薄储层的产能预测精度, 本文在前人测井相研究的基础上提出测井产能等级指数, 并在数据和模型双核驱动下将测井产能等级分类结果与地震属性相结合进行产能预测。首先在降维和聚类的机器学习算法下进行单井的测井产能等级自动划分, 然后以平面径向流公式为基础建立测井产能等级指示模型, 以测井产能等级指数表征单井的产能大小, 最后与地震属性相结合, 采用支持向量回归算法将测井产能等级指数从单点外推到平面上, 实现了研究区滩坝砂储层的产能等级横向预测。

1 基于机器学习的测井产能等级自动划分

机器学习是数据分析的重要技术之一[15-16], 按处理任务的不同, 机器学习可分为回归、分类、降维和聚类四大类。测井产能等级划分是测井产能分析中至关重要的步骤, 测井产能等级划分结果的合理性直接影响后续分析的准确性。本文以测井数据为对象, 采用机器学习算法对地层测井产能等级进行自动划分。其主要步骤包括: 测井资料预处理、测井数据降维处理、测井数据降维后聚类分析、测井产能等级划分。

1.1 测井资料预处理

测井资料预处理包括环境校正、深度校正、测井曲线标准化等, 目的是使测井曲线尽可能真实地反映地层的岩性与孔隙流体等性质, 消除非地层因素的影响[17-18]。本文对测井资料进行了曲线质量检验、岩心深度归位和测井资料标准化等预处理。以声波测井资料标准化为例, 图 1给出了标准化前、后各井声波时差曲线的频率分布。由于不同的测井仪器以及人为操作等原因, 标准化前各井曲线频率的特征峰值一致性较差, 而经过标准化后特征峰值一致性较好, 消除了不同井间测井系列差异等非地层因素造成的误差。

图 1 标准化前(a)、后(b)声波时差曲线频率分布(1 ft≈0.304 8 m)
1.2 测井数据降维处理

为了尽可能准确地反映井下地层的特征, 往往使用多种测井曲线进行分析, 但不同的测井数据在反映地层信息方面存在重复性, 且将各种测井曲线都考虑在内可能导致样本的维度过大, 增加数据分析和处理的难度。主成分分析是机器学习中应用非常广泛的一种降维算法, 可以将m条测井曲线构造成n个主成分(m>n), 每个主成分都是原始测井曲线的线性组合, 这样既可以降低待处理数据的维度, 又能最大程度地保留原始数据包含的信息, 避免各个参数之间存在重复信息。从测井数据中选取自然伽马(GR)、自然电位(SP)、补偿中子(CNL)、补偿密度(DEN)、声波时差(AC)、深电阻率(RT)、浅电阻率(RS)共7条常规测井曲线作为原始输入数据进行主成分分析, 得到各主成分的特征向量和特征值以及累计贡献率(表 1)。从表 1中可以看出, 前4个主成分的累计贡献率达93.1%, 基本包含了原始7条曲线携带的地层信息, 同时也可以降低后续聚类分析的复杂度。

表 1 各主成分的特征向量和特征值及累计贡献率
1.3 降维后数据聚类分析

K-means聚类算法是机器学习中一种无监督学习算法, 本文将其应用于降维后测井数据的产能等级自动分类任务中。选择工区内的Li57, Li67, Li671和Li672 4口井作为测井产能等级划分建模标准井, 将4口井的主成分进行K-means聚类, 主要步骤为: ①随机选取k个点作为聚类中心; ②根据(1)式计算每个点分别到k个聚类中心的欧氏距离, 然后将该点分到最近的聚类中心, 这样就形成了k个簇; ③重新计算每个簇的质心(均值); ④重复步骤②和步骤③, 直到质心的位置不再发生变化或者达到设定的迭代次数。肘部法则对K-means算法的类数确定起指导作用。假设一共将n个样本点划分到k个类中(1≤k < n), 用Ci表示第i类测井产能等级, 且该类质心的位置记为ui, Ji为第i类的畸变程度((2)式), J为所有类的总畸变程度((3)式), 设定聚类类别数的最大值为m, 依次计算k=1, 2, 3, …, m时聚类后的总畸变程度。设定最大聚类数为19, 计算每个聚类数对应的总畸变程度。图 2为K-means聚类的类别数与总畸变程度的折线图。当聚类的类别数从1增大到4时, 畸变程度变化最大。超过4类后, 总畸变程度变化显著降低。因此确定肘部k=4, 将测井产能等级划分为4类。

$ \begin{array}{l} d\left( {x,y} \right) = \\ \sqrt {{{\left( {{x_1} - {y_1}} \right)}^2} + {{\left( {{x_2} - {y_2}} \right)}^2} + \cdots + {{\left( {{x_n} - {y_n}} \right)}^2}} \end{array} $ (1)
$ {J_i} = \sum\limits_{x \in {C_i}} {{{\left| {x - {u_i}} \right|}^2}} $ (2)
$ J = \sum\limits_{i = 1}^k {\sum\limits_{x \in {C_i}} {{{\left| {x - {u_i}} \right|}^2}} } $ (3)
图 2 聚类的类别数与总畸变程度折线显示结果

式中: d(x, y)为点到聚类中心的欧氏距离; xi为待分类点的第i个主成分; yi为聚类中心第i个主成分。

1.4 单井测井产能等级自动划分

对工区内的非标准井逐点计算其与K-means聚类后各个聚类中心的欧氏距离, 将当前点划入与其距离最近的类, 按照相同的方法和步骤依次完成工区内所有井的单井测井产能等级自动划分。图 3显示了Li673井的测井产能等级划分结果。工区内目的层测井产能由高到低划分为4个等级, 将取心井的岩心资料与测井产能等级划分结果进行统计分析(表 2), 可以看出, 不同测井产能等级的地层其岩性和物性差异较为明显。

图 3 Li673井测井产能等级划分结果(1 in≈2.54 cm)
表 2 不同测井产能等级地层的测井、岩性和物性特征
2 测井产能等级指示模型的建立

目前常规砂泥岩储层产能计算多采用平面径向流模型[19]:

$ q = \frac{{2{\rm{ \mathsf{ π} }}Kh\left( {{p_{\rm{e}}} - {p_{\rm{w}}}} \right)}}{{\mu B\left( {\ln \frac{{{r_{\rm{e}}}}}{{{r_{\rm{w}}}}} + S} \right)}} $ (4)

式中: q为地面原油产量; K为地层渗透率; h为地层厚度; μ为原油黏度; B为原油体积系数; S为表皮系数, 由试井解释结果获得; pe为油层压力; pw为井底流压; re为供油半径; rw为井筒半径。

平面径向流公式中参数较多, 在滩坝砂这类砂泥岩薄互型储集体的实际应用中较为困难。因此, 令

$ \frac{{2{\rm{ \mathsf{ π} }}K\left( {{p_{\rm{e}}} - {p_{\rm{w}}}} \right)}}{{\mu B\left( {\ln \frac{{{r_{\rm{e}}}}}{{{r_{\rm{w}}}}} + S} \right)}} = \alpha $ (5)

则(4)式简化为:

$ q = \alpha \times h $ (6)

式中: α是与原油性质、压力、井筒属性以及储层性质等有关的变量, 假设同一工区井中目的层的原油性质、压力、井筒属性等都相同, 则α值只与储层性质有关。

诸多学者研究认为, 测井相是指在地下具有相似的储层性质的一组测井数据集, 根据测井相可以将地层划分为多个具有不同地质意义的纵向条带[11-12, 20]。因此做出假设: 具有相同测井产能等级类型的地层具有相同的α值, 不同测井产能等级地层的α值不同。由于滩坝砂储层是典型的砂泥岩薄互型储层, 通常一段储层具有多种测井产能等级类型, 因此, 可将一段储层的产能大小看成是受厚度与测井产能等级共同影响的量。其中, α可看成不同类型测井产能等级的地层对储层产能的贡献率。基于此, 本文提出测井产能等级指示模型, 通过测井产能等级指数表征储层产能大小。测井产能等级指示模型表达式如下:

$ Q = {\alpha _1}{h_1} + {a_2}{h_2} + \cdots + {a_n}{h_n} = \sum\limits_{i = 1}^n {{\alpha _i}{h_i}} $ (7)

式中: Q为测井产能等级指数; h1, h2, …, hn分别为第1, 2, …, n级测井产能地层的累计厚度; α1, α2, …, αn分别为第1, 2, …, n级测井产能地层对产能的贡献率, 其值采用最小二乘法求得。多组试油层的日均产油量与各级测井产能地层厚度满足以下函数关系:

$ \left[ {\begin{array}{*{20}{c}} {{h_{11}}}&{{h_{12}}}& \cdots &{{h_{1n}}}\\ {{h_{21}}}&{{h_{22}}}& \cdots &{{h_{2n}}}\\ \vdots & \vdots &{}& \vdots \\ {{h_{m1}}}&{{h_{m2}}}& \cdots &{{h_{mn}}} \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {{\beta _1}}\\ {{\beta _2}}\\ \vdots \\ {{\beta _n}} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} {{y_1}}\\ {{y_2}}\\ \vdots \\ {{y_m}} \end{array}} \right] $ (8)

式中: hij为第i个试油层第j级测井产能地层的厚度; βj为第j级测井产能地层对产能的贡献程度; yi为第i个试油层的日均产油量。求解目标函数

$ {\rm{loss}} = \sum\limits_{i = 1}^m {{{\left( {{y_i} - \sum\limits_{j = 1}^n {{h_{ij}}{\beta _j}} } \right)}^2}} $ (9)

loss取最小值时的β值即为(8)式中β的最优估计值。此时, 不同测井产能等级的地层对产能的贡献率为:

$ {\alpha _j} = \frac{{{\beta _j}}}{{\sum\limits_{j = 1}^n {{\beta _j}} }} $ (10)

αj可由(8)式的解β1, β2, …βn经(10)式求得。在实际处理中试油层数通常大于测井产能等级数, 因此(8)式通常是一个超定方程组, 无唯一解, 通过最小二乘拟合求得该方程组的最小二乘解, 即(9)式中loss最小时对应的β值。求解出每一类测井产能等级的地层对应的α值就可以通过测井产能等级指数表征储层产能的大小。

研究区的测井产能等级被划分为4类, 因此其测井产能等级指数可表示为:

$ Q = {\alpha _1}{h_1} + {\alpha _2}{h_2} + {\alpha _3}{h_3} + {\alpha _4}{h_4} $ (11)

分别计算各等级产能对应地层的累计厚度, 将试油层测井产能甜点相剖面转化为累计厚度分别为h1, h2, h3, h4的等相厚度剖面。统计了研究区内8个试油层的测井产能等级地层厚度与日均产油量的数据(表 3)。由研究区内测井相与取心资料分析可知, 第4级产能储层的岩性主要为泥岩, 其通常为储层段中间的薄夹层, 对产能无贡献, 即β4=0。因此, 可以将前3级产能的储层对产能贡献率的求解问题转化为求超定方程组最优解的问题:

$ \left[ {\begin{array}{*{20}{c}} {{h_{11}}}&{{h_{12}}}&{{h_{13}}}\\ {{h_{21}}}&{{h_{22}}}&{{h_{23}}}\\ \vdots & \vdots & \vdots \\ {{h_{81}}}&{{h_{82}}}&{{h_{83}}} \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {{\beta _1}}\\ {{\beta _2}}\\ {{\beta _3}} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} {{y_1}}\\ {{y_2}}\\ \vdots \\ {{y_8}} \end{array}} \right] $ (12)
表 3 地层厚度与日均产油量关系

表 3中研究区内实际产能数据代入(12)式, 求得方程组的解为β1=5.211, β2=1.586, β3=0.696, 由(13)式求得α1=0.696, α2=0.212, α3=0.093。

$ {\alpha _j} = \frac{{{\beta _j}}}{{\sum\limits_{j = 1}^3 {{\beta _j}} }} $ (13)

由此可以导出研究区内测井产能等级指数的实际公式为:

$ Q = 0.696{h_1} + 0.212{h_2} + 0.093{h_3} $ (14)
3 储层产能等级地震属性横向预测

仅利用单井的测井产能等级指数信息难以进行储层产能横向预测, 结合地震资料可进行测井产能等级的平面展布研究。地震属性能描述和量化地震资料的特性, 是原始地震资料中所包含全部信息的子集[21-22], 其中包含了大量与岩性、物性、流体性质等相关的信息。因此, 本文利用机器学习中的支持向量回归算法搭建测井产能等级指数与地震属性之间的桥梁, 进行储层产能的横向预测, 具体步骤如下。

3.1 地震属性提取

地震属性提取是本文所提井震结合方法的重要组成部分, 提取的地震属性种类越多, 后续研究中发掘敏感属性的可能性就越大。根据研究区目的层顶、底界面的位置, 从地震数据体中提取了与振幅、相位等相关的20种地震属性, 并读取了井点处不同地震属性的值。

3.2 地震属性优选

地震属性优选应遵循如下原则: ①选择的地震属性与测井产能等级指数具有较好的相关性; ②选择的地震属性两两之间独立。本文采用Pearson相关性分析优选地震属性, 将提取的20种井点处的地震属性值与对应井的测井产能等级指数进行Pearson相关性分析, 结果如表 4所示。其中, 平均峰值振幅、最小振幅、零相位个数、振幅求和与测井产能等级分类结果具有较好的相关性, 且显著相关。为了确保优选的地震属性具有独立性, 对这4个属性两两之间进行Pearson相关性分析(表 5)。由表 5可知, 振幅求和属性与其它3个属性之间具有显著的相关性, 剔除该属性, 剩余的平均峰值振幅、最小振幅、零相位个数即为敏感地震属性。振幅是地震资料岩性解释和储层预测常用的动力学参数, 平均峰值振幅和最小振幅反映了地层的岩性、物性变化, 与储层的含油气性紧密相关[23-24]。相位属性反映地层的岩性变化, 根据测井产能等级指数与敏感地震属性的相关性, 零相位个数越多, 纵向上砂岩层越多的层段更有可能是有利的含油气区。

表 4 测井产能等级指数与地震属性相关性分析结果
表 5 地震属性相关性分析
3.3 基于支持向量回归的储层产能等级横向预测

为了避免单一属性对储层产能预测带来的多解性, 本文采用支持向量回归算法(SVR)建立敏感地震属性与测井产能等级指数之间的非线性映射关系, 将平均峰值振幅、最小振幅、零相位个数3个敏感地震属性(图 4)在井点处的属性值作为支持向量回归的输入, 测井产能等级指数作为支持向量回归的输出, 建立用3种敏感地震属性值预测测井产能等级指数的模型。图 5是利用支持向量回归算法计算的井点处测井产能等级指数与真实值的对比图。从图中可以看出预测值与真实值较为接近, 经统计, 3口未参与建模过程的验证井(Liang751, Liang756, Liang76)预测值与真实值的平均相对误差为6.8%, 预测结果较为准确。最后利用井点处建立的利用敏感地震属性预测测井产能等级指数的支持向量回归模型, 在非井点处预测储层产能等级。图 6是由3种地震属性预测的储层产能等级指数横向平面分布图。

图 4 目的层敏感属性切片 a平均峰值振幅; b最小振幅; c零相位个数
图 5 利用支持向量回归算法计算的井点处测井产能等级指数与真实值对比结果
图 6 由3种地震属性预测的储层产能等级指数横向平面分布
3.4 产能等级划分结果评价

统计了研究区内参与建模的Li67, Li671, Li672, Li673, Li674, Li675 6口井以及未参与建模的Liang751, Liang756, Liang76 3口验证井前6个月的产量数据(图 7图 8)。对于建模井来说, 处于储层产能等级指数平面图中高值区域的Li67, Li673井的产量明显高于处于高值边缘区域的Li671, Li675井的产量(图 7); 对于3口验证井来说, 处于产能等级指数平面图中高值区域内的Liang76井产量明显高于处于低值区域的Liang756, Liang751井的产量(图 8)。实际产能数据与计算得到的储层产能等级指数平面图吻合较好。该结果表明, 通过测井产能等级指数平面图可以较好地表征储层产能的横向变化。

图 7 建模井前6个月产量
图 8 验证井前6个月产量
4 结论

本文综合测井与地震资料, 采用机器学习中的降维、聚类和回归算法, 基于平面径向流公式建立了测井产能等级指示模型, 在数据和模型双核驱动下以井震结合的方式对滩坝砂储层的产能等级进行了预测, 得到以下3点认识:

1) 聚类分析中的肘部法可以解决测井产能等级分类中的最佳分类数问题, 在测井产能等级划分中无需与岩心数据对比人为确定聚类数, 避免了由主观因素造成的误差, 简化了测井产能等级划分的研究过程;

2) 研究过程中自动划分的各类测井产能等级地层具有不同的储层物性、含油气性等性质, 对储层产能的贡献率各不相同, 进一步在平面径向流公式基础上建立的测井产能等级指示模型可以表征单井产能的大小;

3) 通过建立测井产能等级指数与敏感地震属性之间的映射关系, 将测井产能等级指数从单点外推到平面实现了研究区储层的产能横向预测, 为滩坝砂砂泥岩薄互型储层的产能预测提供了一种新的基于数据和模型双驱动的井震结合解决方案。

参考文献
[1]
韩长城, 林承焰, 任丽华, 等. 地震波形指示反演在东营凹陷王家岗地区沙四上亚段滩坝砂的应用[J]. 中国石油大学学报(自然科学版), 2017, 41(2): 60-69.
HAN C C, LIN C Y, REN L H, et al. Application of seismic waveform inversion in Ess 4 beach-bar sandstone in Wangjiagang area, Dongying Depression[J]. Journal of China University of Petroleum(Edition of Natural Sciences), 2017, 41(2): 60-69. DOI:10.3969/j.issn.1673-5005.2017.02.007
[2]
朱卫星, 徐大年, 曹自强, 等. 基于扩展神经网络算法的弱信号分离技术在滩坝砂储层预测中的应用[J]. 天然气勘探与开发, 2015, 38(3): 39-42.
ZHU W X, XU D N, CAO Z Q, et al. Application of weak-signal separation technology based on extended neural network algorithm to predicting beach-bar sandstone reservoirs[J]. Natural Gas Exploration and Development, 2015, 38(3): 39-42. DOI:10.3969/j.issn.1673-3177.2015.03.009
[3]
袁红军. 东营凹陷博兴洼陷滨浅湖相滩坝砂岩储层预测[J]. 石油与天然气地质, 2007, 28(4): 497-503.
YUAN H J. Prediction of beach bar sand reservoirs of shore-shallow lake facies in Boxing subsag of Dongying sag[J]. Oil and Gas Geology, 2007, 28(4): 497-503. DOI:10.3321/j.issn:0253-9985.2007.04.009
[4]
关新, 潘树新, 曲永强, 等. 准噶尔盆地沙湾凹陷滩坝砂的发现及油气勘探潜力[J]. 岩性油气藏, 2021, 33(1): 90-98.
GUAN X, PAN S X, QU Y Q, et al. Discovery and hydrocarbon exploration potential of beach-bar sand in Shawan Sag, Junggar Basin[J]. Lithologic Reservoirs, 2021, 33(1): 90-98.
[5]
刘双莲, 李浩, 魏修平, 等. 含岩屑砂岩储层"甜点"测井评价技术难点与对策[J]. 石油物探, 2020, 59(1): 141-149.
LIU S L, LI H, WEI X P, et al. Difficulties and strategy of logging evaluation on "sweet spots" of lithic sandstone reservoirs[J]. Geophysical Prospecting for Petroleum, 2020, 59(1): 141-149. DOI:10.3969/j.issn.1000-1441.2020.01.016
[6]
刘力辉, 李建海, 刘玉霞. 地震物相分析方法与"甜点"预测[J]. 石油物探, 2013, 52(4): 432-437.
LIU L H, LI J H, LIU Y X. Seismic reservoir property facies analysis and sweet spot prediction[J]. Geophysical Prospecting for Petroleum, 2013, 52(4): 432-437. DOI:10.3969/j.issn.1000-1441.2013.04.014
[7]
冯进, 赵冰, 张占松, 等. 珠江口盆地惠州凹陷储层测井产能分级与识别方法[J]. 物探与化探, 2020, 44(1): 81-87.
FENG J, ZHAO B, ZHANG Z S, et al. Classification and identification method of reservoir logging capacity in Huizhou depression of Pearl River mouth basin[J]. Geophysical and Geochemical Exploration, 2020, 44(1): 81-87.
[8]
成家杰, 王文文, 侯振学, 等. 常规测井资料在致密砂岩储层产能评价中的应用[J]. 测井技术, 2017, 41(3): 325-330.
CHENG J J, WANG W W, HOU Z X, et al. Application of conventional log data to productivity evaluation of tight sandstone gas reservoirs[J]. Well Logging Technology, 2017, 41(3): 325-330.
[9]
冯春珍, 林伟川, 成志刚, 等. 低渗透储层测井分类和产能预测技术[J]. 测井技术, 2012, 36(3): 308-313.
FENG C Z, LIN W C, CHENG Z G, et al. Reservoir classification and productivity forecasting method for low porosity and permeability reservoir with log data[J]. Well Logging Technology, 2012, 36(3): 308-313. DOI:10.3969/j.issn.1004-1338.2012.03.019
[10]
李大军, 杨晓, 王小兰, 等. 四川盆地W地区龙马溪组页岩气压裂效果评估和产能预测研究[J]. 石油物探, 2017, 56(5): 735-745.
LI D J, YANG X, WANG X L, et al. Estimating the fracturing effect and production capacity of the Longmaxi Formation of the Lower Silurian in area W, Sichuan Basin[J]. Geophysical Prospecting for Petroleum, 2017, 56(5): 735-745. DOI:10.3969/j.issn.1000-1441.2017.05.014
[11]
杜启振, 杨少春, 王志欣, 等. 测井相模式识别自动分析[J]. 石油物探, 1997, 36(增刊1): 108-111.
DU Q Z, YANG S C, WANG Z X, et al. Automatic analysis of pattern recognition of logging facies[J]. Geophysical Prospecting for Petroleum, 1997, 36(S1): 108-111.
[12]
雍世和, 陈钢花. 自动测井相分析[J]. 石油学报, 1990, 11(4): 41-50.
YONG S H, CHEN G H. Automatic logging facies analysis[J]. Acta Petrolei Sinica, 1990, 11(4): 41-50. DOI:10.3321/j.issn:0253-2697.1990.04.009
[13]
杨占龙, 彭立才, 陈启林, 等. 地震属性分析与岩性油气藏勘探[J]. 石油物探, 2007, 46(2): 131-136.
YANG Z L, PENG L C, CHEN Q L, et al. Seismic attribute analysis and lithological reservoir exploration[J]. Geophysical Prospecting for Petroleum, 2007, 46(2): 131-136. DOI:10.3969/j.issn.1000-1441.2007.02.006
[14]
王永刚, 谢东, 乐友喜, 等. 地震属性分析技术在储层预测中的应用[J]. 中国石油大学学报(自然科学版), 2003, 27(3): 30-32.
WANG Y G, XIE D, YUE Y X, et al. Application of seismic attribute analysis technology in reservoir prediction[J]. Journal of China University of Petroleum(Edition of Natural Sciences), 2003, 27(3): 30-32. DOI:10.3321/j.issn:1000-5870.2003.03.009
[15]
ZHONG Z, SUN A Y, WANG Y, et al. Predicting field production rates for waterflooding using a machine learning-based proxy model[J]. Journal of Petroleum Science and Engineering, 2020, 194: 107574.
[16]
SHIRANGI M G. Applying machine learning algorithms to oil reservoir production optimization[R]. California: Stanford University, 2012
[17]
雍世和, 张超谟. 测井数据处理与综合解释[M]. 东营: 中国石油大学出版社, 2002: 90-91.
YONG S H, ZHANG C M. Logging data processing and comprehensive interpretation[M]. Dongying: China University of Petroleum Press, 2002: 90-91.
[18]
肖燕, 黄怡, 屈伟, 等. 测井资料预处理对井约束反演影响分析及应用[J]. 新疆地质, 2007, 25(4): 418-421.
XIAO Y, HUANG Y, QU W, et al. The influence analysis and application of well data pretreatment on well restrict inversion[J]. Xinjiang Geology, 2007, 25(4): 418-421. DOI:10.3969/j.issn.1000-8845.2007.04.017
[19]
孙建孟, 运华云, 冯春珍. 测井产能预测方法与实例[J]. 测井技术, 2012, 36(6): 628-634.
SUN J M, YUN H Y, FEN C Z. The methods and instances of productivity prediction using log data[J]. Well Logging Technology, 2012, 36(6): 628-634. DOI:10.3969/j.issn.1004-1338.2012.06.017
[20]
DELFINER P, PEYRET O, SERRA O. Automatic determination of lithology from well logs[J]. SPE Formation Evaluation, 1987, 2(3): 303-310.
[21]
王开燕, 徐清彦, 张桂芳, 等. 地震属性分析技术综述[J]. 地球物理学进展, 2013, 28(2): 815-823.
WANG K Y, XU Q Y, ZHANG G F, et al. Summary of seismic attribute analysis[J]. Progress in Geophysics, 2013, 28(2): 815-823.
[22]
郭华军, 刘庆成. 地震属性技术的历史、现状及发展趋势[J]. 物探与化探, 2008, 32(1): 19-22.
GUO H J, LIU Q C. The discussion of earthquake attribute technology's history, present situation and development tendency[J]. Geophysical and Geochemical Exploration, 2008, 32(1): 19-22.
[23]
邹才能, 张颖. 油气勘探开发实用地震新技术[M]. 北京: 石油工业出版社, 2002: 202-218.
ZOU C N, ZHANG Y. The new practical seismic techniques in hydrocarbon exploration and development[M]. Beijing: Petroleum Industry Press, 2002: 202-218.
[24]
李坦, 殷小舟. 地震属性的地质意义分析[J]. 复杂油气藏, 2009, 2(3): 25-28.
LI T, YIN X Z. Analysis on the geological significance of seismic attributes[J]. Complex Hydrocarbon Reservoirs, 2009, 2(3): 25-28.