CN113671594B - 基于bp神经网络的碳酸盐岩高频层序自动化识别方法 - Google Patents
基于bp神经网络的碳酸盐岩高频层序自动化识别方法 Download PDFInfo
- Publication number
- CN113671594B CN113671594B CN202110684299.7A CN202110684299A CN113671594B CN 113671594 B CN113671594 B CN 113671594B CN 202110684299 A CN202110684299 A CN 202110684299A CN 113671594 B CN113671594 B CN 113671594B
- Authority
- CN
- China
- Prior art keywords
- neural network
- sequence
- frequency
- division
- neurons
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
- G06N3/084—Backpropagation, e.g. using gradient descent
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Computing Systems (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- Biomedical Technology (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- General Engineering & Computer Science (AREA)
- Artificial Intelligence (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Health & Medical Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种基于BP神经网络的碳酸盐岩高频层序自动化识别方法,包括以下步骤:S1、建立高频层序识别划分标准数据库;S2、编写并调试基于python的BP神经网络模型。过对具有丰富层序划分经验的地质人员的层序划分标准进行学习,建立具有强迁移能力的BP神经网络,有效解决了层序划分难度大、周期长的问题,实现了高频层序的自动化划分,能够作为认识海相地层沉积环境演化、辅助油气勘探的重要手段。
Description
技术领域
本发明属于油气勘探技术领域,特别涉及一种基于BP神经网络的碳酸盐岩高频层序自动化识别方法。
背景技术
层序(Sequence):“层序”这一概念对油气勘探具有重要指导意义,自20世纪80年代以来,四大层序地层学派先后兴起。“层序”概念最早由Sloss等学者提出,1977年Vail等人在此基础上与海平面升降相结合进一步建立了经典层序地层学,此后,以Galloway等学者为代表的成因层序地层学派及以Cross为代表的T-R层序地层学派先后兴起。通常情况下,高频层序(High Frequency Sequence)通常被认为是四级及其后续更高级别的旋回。
传统的层序地层学是以不整合面或与之对应的整合面为界限进行层序界面识别,在陆相层序地层研究中,层序边界往往是不整合面,其识别主要是依据地震反射特征(如削截、上超、下超、顶超等),能够在全盆内进行对比进而建立区域上的层序地层格架。常规层序地层划分通常针对碎屑岩颗粒沉积物,沉积过程稳定,容易找到明显的不整合面,因此根据地震反射特征是无法进行细粒沉积岩层序界面识别(图1)。另一方面,在常规油气勘探中,经典层序地层学是建立在湖扩体系或海侵体系,往往被认为是一个完整的沉积过程,很少考虑碳酸盐岩地层沉积体系,利用测井资料难以进行精确的高频层序界面识别。
调研显示,当前对高频层序的划分以地震剖面的人工划分为主,通过寻找地震剖面上的不整合界面,测井、录井数据上的电性、岩性突变处,进一步识别划分体系域,划分层序。目前层序地层划分的指导方法主要分为两种,分别是Vail学派和Cross学派。
Vail层序地层学理论是以地层不整合或与此不整合可以对比的整合面为层序边界,以地震、测井资料等为基础,通过识别地震接触关系和地层叠加样式变化等,利用井震结合的方式,建立盆地或区域等时层序地层格架。
Vail学派层序划分方法可分为以下步骤:
(1)测录井资料综合分析。测录井资料上的响应特征,直接反映了岩性的显著性变化,岩性变化又直接受海平面变化控制,依据测录井资料响应特征变化即可划定层序界面。
(2)地震资料分析。层序界面在地震资料上表现为反射结构和终止关系变化组合特征,如:中—高连续的中—强反射结构特征,区域上地层终止关系主要表现为顶超或削截特征;各层序界面之下的顶超特征均较普遍,表现为前积反射层终止于低角度面上。
(3)多级层序划分。在层序界面识别和追踪对比的基础上,结合研究区垂向上储、盖配置好的特征,考虑到界面的可操作性和可对比性,以海泛面作为四级层序的边界,对研究区四级层序进行了划分。通过识别追踪对比测井资料的海泛面,井震结合,划分出多级层序。
Vail层序地层学理论通过地震、测井技术剖面人工识别、划分体系域,难度较高,人为因素影响较大。具体缺点总结如下:
(1)体系复杂。Vail学派总结的体系域放分方案包含低位、水进、高位体系域,体系域之下地层又分为进积、退积、加积等多种组合样式,分类复杂。
(2)人为干扰因素多。地震剖面的解释划分完全需要人工进行,不同人员进行解释时往往会产生较大差距。此外,在地震剖面上划分层序界面及后续多级层序时,完全依靠解释人员的经验,对于缺乏经验的解释人员通常不知如何下手,工作难度大,影响层序划分的精度。
Cross高分辨率层序地层研究方法是以基准面旋回为理论依据,强调基准面变化对层序发育的控制作用,将基准面“升—降”旋回期间形成的地层作为一套层序,将基准面旋回变化的转换点作为层序边界;以岩心、三维露头、测井和高分辨率地震剖面等为基础,以基准面旋回与过程-响应原理为理论依据,运用精细层序划分对比技术,建立区域、油田乃至油藏级储层的高分辨率层序地层对比格架。
Cross学派层序划分方法可分为以下步骤:
(1)首先利用钻井岩心,确定研究区超短期基准面旋回类型及特征。超短期基准面旋回类型及特征主要为向上变深的非对称型(Ⅰ型)和向上变浅的非对称型(Ⅱ型)2种超短期旋回,进一步可分为低可容空间向上变浅的非对称型Ⅱ-1和高可容空间向上变浅的非对称型Ⅱ-2型。
(2)然后用钻井标定测井,根据测井曲线反映的地层叠加样式,确定研究区的短期旋回。短期基准面旋回类型及特征主要为对称型层序。此处水下分流河道、河口坝主要发育在基准面对称类型中的下降半旋回中,自然伽马曲线分别呈箱形和漏斗形。
(3)确定中期基准面旋回特征。依据短期旋回的叠加样式,根据基准面旋回叠加理论,可进一步确定中期基准面旋回特征。
Vail和Cross的层序划分方案虽不相同,但操作过程均需要大量人工操作,因此,Cross学派方案存在与Vail学派方案相同的缺点。
发明内容
本发明的目的在于克服现有技术的不足,提供一种通过对具有丰富层序划分经验的地质人员的层序划分标准进行学习,建立具有强迁移能力的BP神经网络,实现了高频层序的自动化划分的基于BP神经网络的碳酸盐岩高频层序自动化识别方法。
本发明的目的是通过以下技术方案来实现的:基于BP神经网络的碳酸盐岩高频层序自动化识别方法,包括以下步骤:
S1、建立高频层序识别划分标准数据库:通过在少量生产井的地层综合柱状图上利用岩性、测井曲线、钻井测试信息,建立高频旋回的人工划分方案;在人工划分的旋回基础上,利用One-Hot编码的方案,将正旋回编码为(0,1),将反旋回编码为(1,0),将高频旋回编码结果与测井曲线一起建立高频层序识别划分标准数据库;
S2、编写并调试基于python的BP神经网络模型,具体步骤如下:
S21、数据预处理:利用sklearn中preprocessing对测井数据进行归一化处理,其计算方法为:
Pi为归一化后的测井数据,Li为原始测井数据,Lmax和Lmin分别为测井数据的最大值和最小值;
S22、利用shuffle函数将测井数据打乱;
S23、建立神经网络模型,神经网络模型包含5个Dense全连接层:第一个Dense全连接层包含1000个神经元,激活函数为elu;第二个Dense全连接层包含500个神经元,激活函数为selu;第三个Dense全连接层包含250个神经元,激活函数为tanh;第四个Dense全连接层包含125个神经元,激活函数为softsign;第五个Dense全连接层包含2个神经元,激活函数为softmax;神经网络模型的损失函数采用交叉熵损失函数,具体公式如下:
其中,f(n)表示交叉熵,n表示所有类别的个数,本发明对应正反两种旋回;y是每种旋回内训练的数据点个数,a是预测的准确率;
S24、采用牛顿梯度下降法对神经元的权重进行动态调整。
本发明的有益效果是:本发明通过对具有丰富层序划分经验的地质人员的层序划分标准进行学习,建立具有强迁移能力的BP神经网络,有效解决了层序划分难度大、周期长的问题,实现了高频层序的自动化划分,能够作为认识海相地层沉积环境演化、辅助油气勘探的重要手段。
附图说明
图1为利用地震剖面识别划分层序示意图;
图2为本发明的碳酸盐岩高频层序自动化识别方法的流程图。
具体实施方式
下面结合附图进一步说明本发明的技术方案。
如图2所示,本发明的一种基于BP神经网络的碳酸盐岩高频层序自动化识别方法,包括以下步骤:
S1、建立高频层序识别划分标准数据库:标准数据库根据少量井的生产数据建立,通过在少量生产井的地层综合柱状图上利用岩性、测井曲线、钻井测试信息,建立高频旋回的人工划分方案;在人工划分的旋回基础上,利用One-Hot编码的方案,将正旋回编码为(0,1),将反旋回编码为(1,0),将高频旋回编码结果与测井曲线一起建立高频层序识别划分标准数据库;
本实施例层序识别划分标准数据库建立依托中国石油天然气集团勘探开发研究院提供的高频层序划分结果,部分数据如表1所示。
表1
顶深m | 底深m | 旋回 | DT | GR | RHOB |
1542.796 | 1562.3 | 反旋回 | 98 | 20.411 | 2.302 |
1562.3 | 1586.84 | 正旋回 | 99 | 20.788 | 2.294 |
1586.84 | 1604.12 | 反旋回 | 99 | 21.164 | 2.287 |
1604.12 | 1617.844 | 正旋回 | 100 | 21.457 | 2.282 |
1617.844 | 1704.875 | 反旋回 | 100 | 21.75 | 2.277 |
1704.875 | 1756.91 | 正旋回 | 101 | 22.044 | 2.273 |
1756.91 | 1772.875 | 反旋回 | 101 | 22.337 | 2.268 |
1772.875 | 1888.232 | 正旋回 | 101 | 22.528 | 2.265 |
S2、编写并调试基于python的BP神经网络模型,具体步骤如下:
S21、数据预处理:利用sklearn中preprocessing对测井数据进行归一化处理,其计算方法为:
Pi为归一化后的测井数据,Li为原始测井数据,Lmax和Lmin分别为测井数据的最大值和最小值;
S22、利用shuffle函数将测井数据打乱;Shuffle函数是对建立的训练集进行打乱,打乱的目的是减少神经网络训练的过拟合。本实施例中采用基于25%采样率随机切片,即将75%的样本划分为训练集,将25%的样本划分为验证集。
S23、建立神经网络模型,神经网络模型包含5个Dense全连接层:第一个Dense全连接层包含1000个神经元,激活函数为elu;第二个Dense全连接层包含500个神经元,激活函数为selu;第三个Dense全连接层包含250个神经元,激活函数为tanh;第四个Dense全连接层包含125个神经元,激活函数为softsign;第五个Dense全连接层包含2个神经元,激活函数为softmax;Softmax是对分类函数,判定类别对应Mishrif组地层所有的岩性类别。
数据经过5个全连接层后进入损失函数,本发明的神经网络模型的损失函数采用交叉熵损失函数(categorical_crossentropy),交叉熵是用来评估当前训练得到的概率分布与真实分布的差异情况,刻画了实际输出(概率)与期望输出(概率)的距离,也就是交叉熵的值越小,两个概率分布就越接近;具体公式如下:
其中,f(n)表示交叉熵,n表示所有类别的个数,本发明对应正反两种旋回;y是每种旋回内训练的数据点个数,a是预测的准确率;
S24、采用牛顿梯度下降法对神经元的权重进行动态调整。牛顿法是一种在实数域和复数域上近似求解方程的方法,使用损失函数的泰勒级数的前面几项来寻找损失函数方程=0的根。首先是所有神经元权重均为1时,计算此时的交叉熵函数f(1),计算此时损失函数的斜率。再结算以上一步计算的斜率,穿过(1,f(1))的点与损失函数的交点x1,将x1带入损失函数,重复上述步骤。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合,这些变形和组合仍然在本发明的保护范围内。
Claims (1)
1.基于BP神经网络的碳酸盐岩高频层序自动化识别方法,其特征在于,包括以下步骤:
S1、建立高频层序识别划分标准数据库:通过在少量生产井的地层综合柱状图上利用岩性、测井曲线、钻井测试信息,建立高频旋回的人工划分方案;在人工划分的旋回基础上,利用One-Hot编码的方案,将正旋回编码为(0,1),将反旋回编码为(1,0),将高频旋回编码结果与测井曲线一起建立高频层序识别划分标准数据库;
S2、编写并调试基于python的BP神经网络模型,具体步骤如下:
S21、数据预处理:利用sklearn中preprocessing对测井数据进行归一化处理,其计算方法为:
Pi为归一化后的测井数据,Li为原始测井数据,Lmax和Lmin分别为测井数据的最大值和最小值;
S22、利用shuffle函数将测井数据打乱;
S23、建立神经网络模型,神经网络模型包含5个Dense全连接层:第一个Dense全连接层包含1000个神经元,激活函数为elu;第二个Dense全连接层包含500个神经元,激活函数为selu;第三个Dense全连接层包含250个神经元,激活函数为tanh;第四个Dense全连接层包含125个神经元,激活函数为softsign;第五个Dense全连接层包含2个神经元,激活函数为softmax;
神经网络模型的损失函数采用交叉熵损失函数,具体公式如下:
其中,f(n)表示交叉熵,n表示所有类别的个数,对应正反两种旋回;y是每种旋回内训练的数据点个数,a是预测的准确率;
S24、采用牛顿梯度下降法对神经元的权重进行动态调整。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110684299.7A CN113671594B (zh) | 2021-06-21 | 2021-06-21 | 基于bp神经网络的碳酸盐岩高频层序自动化识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110684299.7A CN113671594B (zh) | 2021-06-21 | 2021-06-21 | 基于bp神经网络的碳酸盐岩高频层序自动化识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113671594A CN113671594A (zh) | 2021-11-19 |
CN113671594B true CN113671594B (zh) | 2023-06-20 |
Family
ID=78538230
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110684299.7A Active CN113671594B (zh) | 2021-06-21 | 2021-06-21 | 基于bp神经网络的碳酸盐岩高频层序自动化识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113671594B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117991335A (zh) * | 2022-10-27 | 2024-05-07 | 中国石油天然气股份有限公司 | 基于常规测井资料的高频旋回碳酸盐岩识别方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2004008389A1 (en) * | 2002-07-12 | 2004-01-22 | Chroma Energy, Inc. | Pattern recognition applied to oil exploration and production |
CN1975462A (zh) * | 2006-09-05 | 2007-06-06 | 孟召平 | 基于地震属性的煤层厚度分析方法 |
WO2013002233A1 (ja) * | 2011-06-30 | 2013-01-03 | 株式会社 日立メディコ | 磁気共鳴イメージング装置および傾斜磁場波形推定方法 |
CN104035120A (zh) * | 2014-07-02 | 2014-09-10 | 成都理工大学 | 粉末环境样品α放射性核素的测量方法及测量装置 |
CA3075770A1 (en) * | 2017-09-12 | 2019-03-21 | Schlumberger Canada Limited | Seismic image data interpretation system |
CN110187391A (zh) * | 2019-06-20 | 2019-08-30 | 克拉玛依亚新石油科技有限公司 | 一种利用声波时差确定地层层序界面的方法 |
AU2020101809A4 (en) * | 2020-06-09 | 2020-09-17 | China University Of Geosciences, Beijing | A simulation system for depositional sequence formation and evolution |
CN112712025A (zh) * | 2020-12-29 | 2021-04-27 | 东北石油大学 | 一种基于长短期记忆神经网络的复杂岩性识别方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10877178B2 (en) * | 2019-05-03 | 2020-12-29 | Institute Of Geology And Geophysics, Chinese Academy Of Sciences | Method for suppressing airborne transient electromagnetic in-band vibration noise |
-
2021
- 2021-06-21 CN CN202110684299.7A patent/CN113671594B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2004008389A1 (en) * | 2002-07-12 | 2004-01-22 | Chroma Energy, Inc. | Pattern recognition applied to oil exploration and production |
CN1975462A (zh) * | 2006-09-05 | 2007-06-06 | 孟召平 | 基于地震属性的煤层厚度分析方法 |
WO2013002233A1 (ja) * | 2011-06-30 | 2013-01-03 | 株式会社 日立メディコ | 磁気共鳴イメージング装置および傾斜磁場波形推定方法 |
CN104035120A (zh) * | 2014-07-02 | 2014-09-10 | 成都理工大学 | 粉末环境样品α放射性核素的测量方法及测量装置 |
CA3075770A1 (en) * | 2017-09-12 | 2019-03-21 | Schlumberger Canada Limited | Seismic image data interpretation system |
CN110187391A (zh) * | 2019-06-20 | 2019-08-30 | 克拉玛依亚新石油科技有限公司 | 一种利用声波时差确定地层层序界面的方法 |
AU2020101809A4 (en) * | 2020-06-09 | 2020-09-17 | China University Of Geosciences, Beijing | A simulation system for depositional sequence formation and evolution |
CN112712025A (zh) * | 2020-12-29 | 2021-04-27 | 东北石油大学 | 一种基于长短期记忆神经网络的复杂岩性识别方法 |
Non-Patent Citations (1)
Title |
---|
层序地层学研究进展及发展趋势;雷国明;郝洁;陈静静;;科技咨询导报(13);第162和164页 * |
Also Published As
Publication number | Publication date |
---|---|
CN113671594A (zh) | 2021-11-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106886043B (zh) | 基于地震数据深度学习的储层检测方法 | |
CN105259572A (zh) | 基于地震多属性参数非线性自动分类的地震相计算方法 | |
CN111983691B (zh) | 一种多模型融合的储层预测方法及软件系统 | |
CN103472484A (zh) | 基于rs三维敏感地震属性分析的水平井轨迹优化方法 | |
CN113610945B (zh) | 一种基于混合神经网络的地应力曲线预测方法 | |
CN106545337A (zh) | 一种基于支持向量机的沉积微相测井识别方法 | |
CN110275210A (zh) | 一种碳酸盐岩高频层序格架的沉积微相相模式的识别方法 | |
CN110412662A (zh) | 基于地震多属性深度学习的薄互层储层预测方法 | |
CN111983683B (zh) | 一种少井条件下盐下湖相灰岩储层预测方法和系统 | |
CN113671594B (zh) | 基于bp神经网络的碳酸盐岩高频层序自动化识别方法 | |
CN115220101A (zh) | 一种深层碎屑岩砂体沉积结构的建模方法 | |
Liu* et al. | Seismic lithofacies computation method based on deep learning | |
CN114482995B (zh) | 一种细粒沉积物泥质含量精细确定方法 | |
CN115964667A (zh) | 基于深度学习和重采样的河流-湖泊岩相测井识别方法 | |
CN104834934A (zh) | 一种用于识别储层流体的核体俘获法 | |
CN103048688A (zh) | 一种基于三步法的地震属性优选方法 | |
CN112379435B (zh) | 相控岩溶型缝洞集合体刻画方法及装置 | |
CN109765621B (zh) | 一种沉积相精细刻画的成图方法 | |
CN112052591A (zh) | 一种储层构型约束下的隔夹层精细刻画及嵌入式建模方法 | |
CN112147676A (zh) | 一种煤层及夹矸厚度预测方法 | |
CN110443801A (zh) | 一种基于改进AlexNet的盐丘识别方法 | |
CN117665933B (zh) | 一种低分辨率地震资料限制下稀井网区的沉积微相刻画方法 | |
He et al. | Shale seismic facies recognition technology based on sparse autoencoder | |
CN110320567A (zh) | 一种快速恢复碳酸盐岩台地古水深的方法 | |
CN112287532B (zh) | 一种缘控地震驱动建模方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |