CN109143353B - 一种基于深度卷积生成对抗网络的叠前地震波形分类方法 - Google Patents
一种基于深度卷积生成对抗网络的叠前地震波形分类方法 Download PDFInfo
- Publication number
- CN109143353B CN109143353B CN201810946436.8A CN201810946436A CN109143353B CN 109143353 B CN109143353 B CN 109143353B CN 201810946436 A CN201810946436 A CN 201810946436A CN 109143353 B CN109143353 B CN 109143353B
- Authority
- CN
- China
- Prior art keywords
- data
- depth convolution
- confrontation network
- label
- arbiter
- 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
- 238000012549 training Methods 0.000 claims abstract description 15
- 238000001914 filtration Methods 0.000 claims description 3
- 230000009467 reduction Effects 0.000 claims description 3
- 208000035126 Facies Diseases 0.000 abstract description 19
- 238000000034 method Methods 0.000 abstract description 12
- 238000013135 deep learning Methods 0.000 abstract description 2
- 230000006870 function Effects 0.000 abstract description 2
- 230000004913 activation Effects 0.000 description 3
- 238000013527 convolutional neural network Methods 0.000 description 3
- 238000000605 extraction Methods 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- ORILYTVJVMAKLC-UHFFFAOYSA-N Adamantane Natural products C1C(C2)CC3CC1CC2C3 ORILYTVJVMAKLC-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000005137 deposition process Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000009545 invasion Effects 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 239000013049 sediment Substances 0.000 description 1
- 238000004062 sedimentation Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/63—Seismic attributes, e.g. amplitude, polarity, instant phase
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种基于深度卷积生成对抗网络的叠前地震波形分类方法,属于地震相分析领域。本发明根据深度卷积生成对抗网络(DCGAN)对叠前地震波形进行半监督分类,先使用无标签样本让网络学习到叠前地震数据的特征,然后用少量有标签网络精调。本发明能从大量无标签数据中学习数据分布特性,具有很好的特征表示能力。相对于其他半监督方法需要使用多个分类器来增加训练样本,训练方法更简单。相对于其他深度学习特征提取方法,无须启发式损失函数,也能很好地表征图像。
Description
技术领域
本发明属于地震相分析领域,特别涉及一种基于深度卷积生成对抗网络的叠前地震波形分类方法。
背景技术
地震相分析的方法就是通过在划分地震层序的基础上,利用各种地震参数之间的差别以及参数之间的关系,将地震层序划分为不同的区域,然后再进行推断地质结构。地震相分析中应考虑的参数有:反射振幅、主反射频率、反射极性、层速度、反射连续性、反射结构、反射丰度、地震相单元几何、与其它单元的关系。地震数据就是地表检波器接收到的反射信号,然后,将地震信号的细微变化和地下结构信息进行映射,该操作可以通过信号分类技术来完成。地震相数据的解释可以是直接的,也可以是间接的。直接解释的目的是找出引起地震相单元地震特征的地质原因。所以,直接解释可能旨在预测岩性、孔隙度、流体含量,相对年龄,超压页岩、类型分层,对应的地震相单元及其地质背景地质体几何。间接解释的目的是得出一些关于沉积过程和环境、沉积物搬运方向和地质演化(海侵、消退、沉降、隆起、侵蚀)方面的结论。除了提供地震相分类,地震信号分类还可以通过同时评估瞬时属性,相似性及声阻抗的和AVO多属性分析相结合来更好地表达地下信息。地震相分析结果可在地震相剖面和地震相图上显示。根据该区现有的地震资料和地质条件,地震相图可能有不同的类型,如显示不同地震相单元分布的一般地震相图、砂泥岩比图、交错层理方向图和古迁移图等。
叠前地震波是不同方位地表角检波器接收到的原始反射信号,对于检波点都可以利用多重维度的数据来进行描述地下结构信息。叠前地震信号和叠后信号是密切相关的,叠后信号是通过叠加已经用速度模型“校正”或“迁移”的叠前信号获得的。速度模型是从地震时差中获得的,其中运动学可利用的叠前地震事件(通常是初级反射)的偏移。因此叠后数据量较小,数据维度也偏小,失去了原始的信息。如今,大数据技术的快速发展为叠前信号的处理提供了充分的技术支持,从而弥补了以往波形分类算法只能处理叠后信息的不足。
在石油勘探的初期,会产生大量的的叠前地震数据,这些数据可以通过无监督聚类技术来完成地震相分析,从而映射地下结构信息,进而预测以及选择测井的合理位置。而在获得一定数量测井属性后,可以结合测井数据、岩心等对地震相校准。通常使用机器学习中的有监督方法,自动根据测井信息对储层数据进行分类。但是由于测井数据相对地震数据是稀疏的,测井数据只能代表局部地质信息,在传统的有监督分类方法中,分类结果往往较差。
发明内容
为了解决现有技术中的问题,本发明提出了一种基于深度卷积生成对抗网络的叠前地震波形分类方法,基于深度学习技术,围绕地震叠前波形的去噪、特征提取、无监督学习和半监督学习等方面进行研究,研究出如何使用叠前地震波形更好地生成地震相图,有效帮助地质的解释工作。
一种基于深度卷积生成对抗网络的叠前地震波形分类方法,包括以下步骤:
步骤1,对叠前地震数据进行预处理,进行结构导向滤波后根据层位提取样本数据,根据测井位置选取测井邻域数据为有标签数据,其余数据为无标签数据;
步骤2,输入所述无标签数据至深度卷积生成对抗网络进行训练;
步骤3,将所述深度卷积生成对抗网络中判别器的最后一层替换为softmax分类器,构造分类网络模型;
步骤4,输入所述有标签数据至所述分类网络模型进行精调;
步骤5,输入地震工区数据至精调后的分类网络模型,得到所有样本的分类结果和地震相图。
进一步地,所述步骤1包括以下流程:
对叠前地震数据进行结构导向滤波降噪,根据层位提取样本数据;
根据测井位置,选取测井邻域数据为有标签样本,测井的类型为数据样本的标签,其余数据为无标签样本。
进一步地,所述步骤2包括以下流程:
输入所述无标签数据训练所述深度卷积生成对抗网络,所述深度卷积生成对抗网络包括生成器和判别器,生成器的输入为服从均匀分布的噪声矢量,输出为与所述无标签数据相同大小的地震数据;所述无标签数据为所述判别器的输入,所述判别器的输出为一个二分类器。
本发明的有益效果:本发明提供了一种基于深度卷积生成对抗网络的叠前地震波形分类方法,根据深度卷积生成对抗网络(DCGAN)对叠前地震波形进行半监督分类,先使用无标签样本让网络学习到叠前地震数据的特征,然后用少量有标签网络精调。本发明能从大量无标签数据中学习数据分布特性,具有很好的特征表示能力。相对于其他半监督方法需要使用多个分类器来增加训练样本,训练方法更简单。相对于其他深度学习特征提取方法,无须启发式损失函数,也能很好地表征图像。
附图说明
图1为本发明实施例的流程图。
图2为DCGAN的训练网络结构图。
图3为DCGAN的分类网络结构图。
具体实施方式
下面结合附图对本发明的实施例做进一步的说明。
请参阅图1,本发明提供了一种基于深度卷积生成对抗网络的叠前地震波形分类方法,通过以下流程实现:
步骤1,对叠前地震数据进行预处理,进行结构导向滤波后根据层位提取样本数据,根据测井位置选取测井邻域数据为有标签数据,其余数据为无标签数据。
本实施例中,对叠前地震数据进行结构导向滤波降噪,根据层位提取样本数据。
根据测井位置,选取测井邻域数据为有标签样本,测井的类型为数据样本的标签,其余数据为无标签样本。
步骤2,输入所述无标签数据至深度卷积生成对抗网络进行训练。
本实施例中,训练网络为深度卷积生成对抗网络(DCGAN),请参阅图2,其是在生成对抗网络(GAN)的基础上结合卷积神经网络(CNN)得到的一种网络模型,深度卷积生成对抗网络中的生成器和判别器为卷积神经网络,具体有:
用生成器的带步长卷积替换所有的池化层;
用判别器的微步幅卷积替换所有池化层;
在生成器和判别器上都使用批标准化(Batch-normalization),这个策略能有效地解决初始化不当引起训练崩溃的问题,但如果将批标准化应用于所有层又会引起模型的不稳定,所以采取的措施为在生成器的输出层和判别器的输入不适用批标准化;
删除深度网络中的全连接层;
生成器中输出层用Tanh激活函数,其他所有层用ReLU激活函数;
判别器中所有层的激活函数都用LeakyReLU。
本实施例中,对于生成器,输入为100维的均匀噪声,第一层为全连接层,将100维的向量投影成4×4大小的feature map,通道数为512。然后一次用四层步长为3×3的带步长卷积,这样使得每次卷积后图像尺寸加倍,通道数减半。最后转换成32×32的单通道图像,这些图像就是生成的假样本。对于判别器,输入为真实的样本和生成器生成的样本,均为32×32的单通道图像。判别器与生成器的各层的图像尺寸和通道数保持一致。判别器的前四层依次图像尺寸减半,通道数加倍,生成高级特征表示。最后一层为一个logistics回归二分类器,输出为一个标量,即样本为真的概率。训练中超参数设置如下:采用mini-batch进行训练,训练的batchsize为64,采用ADAM优化器进行训练,学习率设置为0.001;设置LeakyReLU的斜率为0.2。
本实施例中,输入无标签数据训练深度卷积生成对抗网络。生成器G的输入为服从均匀分布的噪声矢量,输出为与无标签数据相同大小的地震数据,生成器从随机分布的数据中生成“伪造”的叠前地震波形;判别器D的输入为无标签数据和生成器生成的伪造数据,输出为一个标量,即一个二分类器,输出是否为真实的训练样本,判别器判别输入的地震波形数据,是真实的还是伪造的。通过大量的无标签数据对训练网络进行训练。
步骤3,将所述深度卷积生成对抗网络中判别器的最后一层替换为softmax分类器,构造分类网络模型。
步骤4,输入所述有标签数据至所述分类网络模型进行精调。
本实施例中,输入有标签数据对分类网络模型训练分类器进行精调,如图3所示。
步骤5,输入地震工区数据至精调后的分类网络模型,得到所有样本的分类结果和地震相图。
本实施例中,将地震工区数据输入该分类网络,得到所有样本的分类结果,并绘制地震相图,结合实际地质信息和原理对结果加以分析。
通过两个网络进行博弈交替训练,使生成器生成的样本符合训练样本的概率分布。使用DCGAN进行半监督分类,分类网络会利用样本的分布信息,使得分类效果更好。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合,这些变形和组合仍然在本发明的保护范围内。
Claims (3)
1.一种基于深度卷积生成对抗网络的叠前地震波形分类方法,其特征在于,包括以下步骤:
步骤1,对叠前地震数据进行预处理,进行结构导向滤波后根据层位提取样本数据,根据测井位置选取测井邻域数据为有标签数据,其余数据为无标签数据;
步骤2,输入所述无标签数据至深度卷积生成对抗网络进行训练;
步骤3,将所述深度卷积生成对抗网络中判别器的最后一层替换为softmax分类器,构造分类网络模型;
步骤4,输入所述有标签数据至所述分类网络模型进行精调;
步骤5,输入地震工区数据至精调后的分类网络模型,得到所有样本的分类结果和地震相图。
2.如权利要求1所述的基于深度卷积生成对抗网络的叠前地震波形分类方法,其特征在于,所述步骤1包括以下流程:
对叠前地震数据进行结构导向滤波降噪,根据层位提取样本数据;
根据测井位置,选取测井邻域数据为有标签样本,测井的类型为数据样本的标签,其余数据为无标签样本。
3.如权利要求1所述的基于深度卷积生成对抗网络的叠前地震波形分类方法,其特征在于,所述步骤2包括以下流程:
输入所述无标签数据训练所述深度卷积生成对抗网络,所述深度卷积生成对抗网络包括生成器和判别器,生成器的输入为服从均匀分布的噪声矢量,输出为与所述无标签数据相同大小的地震数据;所述无标签数据为所述判别器的输入,所述判别器的输出为一个二分类器。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810946436.8A CN109143353B (zh) | 2018-08-20 | 2018-08-20 | 一种基于深度卷积生成对抗网络的叠前地震波形分类方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810946436.8A CN109143353B (zh) | 2018-08-20 | 2018-08-20 | 一种基于深度卷积生成对抗网络的叠前地震波形分类方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109143353A CN109143353A (zh) | 2019-01-04 |
CN109143353B true CN109143353B (zh) | 2019-10-01 |
Family
ID=64790255
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810946436.8A Active CN109143353B (zh) | 2018-08-20 | 2018-08-20 | 一种基于深度卷积生成对抗网络的叠前地震波形分类方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109143353B (zh) |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110007341B (zh) * | 2019-02-28 | 2020-10-20 | 长江大学 | 一种基于IfnoGAN和SSD模型的微地震有效信号的识别方法及系统 |
CN110009015A (zh) * | 2019-03-25 | 2019-07-12 | 西北工业大学 | 基于轻量化网络和半监督聚类的高光谱小样本分类方法 |
CN110032975B (zh) * | 2019-04-15 | 2021-09-07 | 禁核试北京国家数据中心 | 一种地震震相的拾取方法 |
CN110320557B (zh) * | 2019-06-10 | 2021-08-17 | 北京有隆科技服务有限公司 | 基于深度学习和进化学习的多尺度地质特征检测融合方法 |
CN112114360B (zh) * | 2019-06-21 | 2024-06-28 | 中国石油天然气集团有限公司 | 一种地震波形分析方法及装置 |
CN110609320B (zh) * | 2019-08-28 | 2021-03-16 | 电子科技大学 | 一种基于多尺度特征融合的叠前地震反射模式识别方法 |
CN110658557B (zh) * | 2019-09-03 | 2021-11-02 | 中国地质大学(北京) | 基于生成对抗网络的地震数据面波压制方法 |
CN110580682A (zh) * | 2019-09-16 | 2019-12-17 | 电子科技大学 | 一种基于优化生成对抗网络地震数据超分辨率重建方法 |
CN110703328B (zh) * | 2019-10-14 | 2021-09-24 | 中海石油(中国)有限公司上海分公司 | 一种超压界面识别方法、装置、设备及存储介质 |
CN112684497B (zh) * | 2019-10-17 | 2023-10-31 | 中国石油天然气集团有限公司 | 地震波形聚类方法和装置 |
CN111242201A (zh) * | 2020-01-07 | 2020-06-05 | 北京师范大学 | 基于对抗生成网络的恒星光谱小样本分类方法 |
CN111402266A (zh) * | 2020-03-13 | 2020-07-10 | 中国石油大学(华东) | 一种构建数字岩心的方法及系统 |
CN111983681B (zh) * | 2020-08-31 | 2021-09-14 | 电子科技大学 | 一种基于对抗学习的地震波阻抗反演方法 |
CN114428325B (zh) * | 2020-09-09 | 2024-08-20 | 中国石油化工股份有限公司 | 一种端到端的对抗网络流体相识别方法 |
CN116821800B (zh) * | 2023-08-31 | 2023-11-10 | 深圳市路桥建设集团有限公司 | 基于半监督生成对抗网络的结构状态分类方法及相关设备 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016073483A1 (en) * | 2014-11-05 | 2016-05-12 | Shell Oil Company | Systems and methods for multi-dimensional geophysical data visualization |
CN107563385A (zh) * | 2017-09-02 | 2018-01-09 | 西安电子科技大学 | 基于深度卷积生成式对抗网络的车牌字符识别方法 |
-
2018
- 2018-08-20 CN CN201810946436.8A patent/CN109143353B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016073483A1 (en) * | 2014-11-05 | 2016-05-12 | Shell Oil Company | Systems and methods for multi-dimensional geophysical data visualization |
CN107563385A (zh) * | 2017-09-02 | 2018-01-09 | 西安电子科技大学 | 基于深度卷积生成式对抗网络的车牌字符识别方法 |
Non-Patent Citations (1)
Title |
---|
unsupervised representation learning with deep convolutional generative adversarial networks;Alec Radford 等;《Computer Science》;20161231;第1-16页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109143353A (zh) | 2019-01-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109143353B (zh) | 一种基于深度卷积生成对抗网络的叠前地震波形分类方法 | |
US11313994B2 (en) | Geophysical deep learning | |
Zhu et al. | Intelligent logging lithological interpretation with convolution neural networks | |
AU2017343749B2 (en) | System and method for seismic facies identification using machine learning | |
US6560540B2 (en) | Method for mapping seismic attributes using neural networks | |
Marano et al. | Generative adversarial networks review in earthquake-related engineering fields | |
EP3894903A1 (en) | Automated reservoir modeling using deep generative networks | |
AU2023200300B2 (en) | Detecting fluid types using petrophysical inversion | |
WO2002029445A1 (en) | Method for seismic facies interpretation using textural analysis and neural networks | |
CN105259572A (zh) | 基于地震多属性参数非线性自动分类的地震相计算方法 | |
GB2592203A (en) | A System and method for improved geophysical data interpretation | |
CN109581487A (zh) | 一种沉积体波形空间结构分析与体控反演建模方法 | |
Ouadfeul et al. | Lithofacies classification using the multilayer perceptron and the self-organizing neural networks | |
CN111983683B (zh) | 一种少井条件下盐下湖相灰岩储层预测方法和系统 | |
Pintea et al. | Seismic inversion with deep learning: A proposal for litho-type classification | |
Xie et al. | First-break automatic picking with fully convolutional networks and transfer learning | |
Li et al. | Using GAN priors for ultrahigh resolution seismic inversion | |
Roy et al. | Application of 3D clustering analysis for deep marine seismic facies classification—An example from deep-water Northern Gulf of Mexico | |
Zahraa et al. | Characterizing geological facies using seismic waveform classification in sarawak basin | |
Bazulin et al. | VTI parameters determination from synthetic sonic logging data using a convolutional neural network | |
Zhang et al. | Self-organizing map (SOM) network for tracking horizons and classifying seismic traces | |
Zhao et al. | Different training sample selection strategies in unsupervised seismic facies analysis | |
Cai et al. | Deep learning for recognition of sedimentary microfacies with logging data | |
Heydarpour et al. | Applying deep learning method to develop a fracture modeling for a fractured carbonate reservoir using geologic, seismic and petrophysical data | |
Wang et al. | cigChannel: A massive-scale 3D seismic dataset with labeled paleochannels for advancing deep learning in seismic interpretation |
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 |