CN107688201B - 基于rbm地震叠前信号聚类方法 - Google Patents
基于rbm地震叠前信号聚类方法 Download PDFInfo
- Publication number
- CN107688201B CN107688201B CN201710728113.7A CN201710728113A CN107688201B CN 107688201 B CN107688201 B CN 107688201B CN 201710728113 A CN201710728113 A CN 201710728113A CN 107688201 B CN107688201 B CN 107688201B
- Authority
- CN
- China
- Prior art keywords
- data
- rbm
- training
- layer
- seismic
- 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
- 238000000034 method Methods 0.000 title claims abstract description 28
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 25
- 230000009467 reduction Effects 0.000 claims abstract description 15
- 238000013528 artificial neural network Methods 0.000 claims abstract description 14
- 238000012549 training Methods 0.000 claims description 31
- 210000002569 neuron Anatomy 0.000 claims description 14
- 230000004913 activation Effects 0.000 claims description 11
- 239000013598 vector Substances 0.000 claims description 10
- 238000010587 phase diagram Methods 0.000 claims description 4
- 230000000007 visual effect Effects 0.000 claims description 4
- 230000003044 adaptive effect Effects 0.000 claims description 2
- 239000011159 matrix material Substances 0.000 claims description 2
- 210000002856 peripheral neuron Anatomy 0.000 claims description 2
- 208000035126 Facies Diseases 0.000 abstract description 12
- 238000000605 extraction Methods 0.000 abstract description 3
- 238000004458 analytical method Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000012880 independent component analysis Methods 0.000 description 4
- 238000009826 distribution Methods 0.000 description 3
- 238000001914 filtration Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000012300 Sequence Analysis Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000007635 classification algorithm Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000000513 principal component analysis Methods 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
- 230000002123 temporal effect 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/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
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
本发明公开了一种基于RBM地震叠前信号聚类方法,本发明的聚类方法针对地震叠前信号维度高的特点,引入基于深度受限玻尔兹曼机网络的降维方法,深度受限玻尔兹曼机具有较好的抽象特征提取能力,降维后的输出能够好的表达原始信号;针对地震信号往往噪声多且难以较好去除的特点,引入模糊自组织神经网络聚类算法,使得地震相分类由隶属度表达,模糊分类结果提供的信息更加丰富合理,使得聚类结果更稳定,而且还能对分类结果进一步探索。本发明的方法有着较低的计算复杂度,能够在大规模地震叠前数据推广。
Description
技术领域
本发明属于地震技术领域,涉及地震信号分类技术,具体涉及一种地震叠前信号聚类方法。
背景技术
通过地震信号分类技术生成地震相图来确定地下储层对指导石油勘探有着重要作用,现有地震信号分类技术主要针对叠后信号,叠后信号是叠前信号的横向求和,导致很多地震信息丢失。地震叠前信号相对于叠后信号显著特点就是维度变高,直接使用现有方法会带来维度灾难,分类算法不适应等问题。同时数据量变大,不能使用计算复杂度较高的算法。在获得测井信息后,测井可以作为有标签数据帮助后续测井的选择。现有方法通常将此作为有监督学习进行处理,但由于测井信息即有标签数据较少,有监督分类结果可能不合实际。
自组织神经网络(Self-organizing Map,SOM)聚类是一种非监督的地震相分析技术,能够自动确定地震相的数量。SOM主要针对地震信号波形形状,也可以选取其他处理过的数据作为聚类对象从而关注某种地质结构,在文献“Steeghs P,DrijkoningenG.Seismic sequence analysis and attribute extraction using quadratic time-frequency representations.Geophysics,2001,66(6):1947-1959”提出使用联合时频分析来表征储层,因为频率特征对反射信息的细微变化敏感。在文献“de Matos M C,OsorioP L,Johann P R.Unsupervised seismic facies analysis using wavelet transformand self-organizing maps.Geophysics,2006,72(1):P9-P21”提出连续小波变换(Continuous Wavelet Transform,CWT)和离散小波变换(Discrete Wavelet Transform,DWT)无时间抽取可以应用于检测奇异性。SOM相对于K均值的主要优点是数据空间中变形的簇可以直接映射到规则网格化的潜在空间。当然SOM算法也存在缺陷,首先,每次迭代的邻域函数的理论基础不足,只能依靠经验选择,选择不同可能最终的聚类结果相差很大;第二,缺少定量误差测量方法,难以判断收敛是否达到最优;第三,算法发现最可能的聚类结果,但没有对相关分类置信度的定量标度,不能说明向量可以被所在簇代表。
高斯混合模型(Gaussian Mixture Model,GMM)是概率分布的参数模型,其可以比传统的无监督聚类算法在建模中提供更大的灵活性和精度。这些GMM岩性用于校准三维地震叠加反演结果,以产生三维岩石性质模型。
在聚类算法中往往需要结合特征降维算法,因为高维的数据会使聚类算法不适应和难以解释。这里简单介绍地质领域的常用降维算法如下:
生成拓扑映射(Generative Topographic Mapping,GTM)是一种基于概率的非线性降维技术,位于较低维度上的每个节点对数据提供一些数学支持并且被认为在某种程度上“负责”数据向量。GTM理论基于概率统计学,可以通过投影多属性数据的特定窗口的平均后验概率,例如测井信息来扩展在地震勘探中的应用。
独立分量分析(Independent Component Analysis,ICA)与基于高斯统计的主成分分析不同,ICA试图将数据投影到子空间上,从而产生非高斯分布,然后更易于分离和可视化。ICA通常用于减少冗余属性集以形成较小的独立元属性集。
不管使用何种地震相分析方法,都需要对储层区域的时间和空间仔细划分。地震相分析对地震波采集误差敏感特别是在储层区域,对结果解释的置信程度取决于地质***的复杂性、地震数据质量和解释者的经验,进而造成与实际情况产生误差。
发明内容
为解决现有技术存在的上述问题,本发明提出了一种基于RBM地震叠前信号聚类方法。
本发明的具体技术方案为:一种基于RBM地震叠前信号聚类方法,具体包括如下步骤:
步骤S1.去除地震叠前数据的噪声,保留结构特征,并将目标层位段数据组织成Xi=[x1;x2;...;xs],其中,x为单个方位角数据,维度为d,s为地震数据的方位角个数;
步骤S2.对n个位置点的数据U=[X1,X2,...,Xn]加入不同的噪声生成m份数据集U1,U2,...,Um;
步骤S3.按深度受限玻尔兹曼机网络对数据提取特征和维度规约;
步骤S4.将降维后的数据作为模糊自组织神经网络的输入,进行无监督训练,得到最终分类结果。
本发明的有益效果:本发明的地震叠前信号聚类方法针对地震叠前信号维度高的特点,引入基于深度受限玻尔兹曼机网络的降维方法,深度受限玻尔兹曼机具有较好的抽象特征提取能力,降维后的输出能够好的表达原始信号;针对地震信号往往噪声多且难以较好去除的特点,引入模糊自组织神经网络聚类算法,使得地震相分类由隶属度表达,模糊分类结果提供的信息更加丰富合理,使得聚类结果更稳定,而且还能对分类结果进一步探索。本发明的方法有着较低的计算复杂度,能够在大规模地震叠前数据推广。
附图说明
图1为本发明实施例的地震叠前信号聚类方法流程示意图。
图2为本发明实施例的深度RBM结构示意图。
图3为现有的方法生成的结果图。
图4为本发明实施例的方法生成的结果图。
图5为本发明实施例的方法生成的隶属度梯度示意图。
具体实施方式
下面结合附图对本发明的实施例做进一步的说明。
本发明实施例的地震叠前信号聚类方法流程如图1所示,具体包括如下步骤:
步骤S1.结构导向滤波算法去除噪声,保留结构特征,并将目标层位段数据组织成后续算法需要的格式。这里通常采用结构导向滤波算法去除噪声,这样能尽量减少噪声对后续处理的影响;
叠前信号每一个位置有多个方位角,即如果单个方位角数据维度d,对于共有s个方位角的地震数据而言,每个位置的特征数据维度为s·d。进行聚类算法之前,需要对地震叠前信号用结构导向滤波算法去除噪声,保留结构特征,并将目标层位段数据组织成后续算法需要的格式Xi=[x1;x2;...;xs],其中,x为单个方位角数据,维度为d。
步骤S2.在叠前数据预处理完成后,对n个位置点的数据U=[X1,X2,...,Xn]加入不同的噪声生成m份数据集U1,U2,...,Um,以便聚类算法更加关注数据的主要特征。
步骤S3.按深度受限玻尔兹曼机网络(Restricted Boltzmann Machines,RBM)对数据提取特征和维度规约,后续聚类算法对于抽象特征分类效果更好。
如图2所示,每层高斯RBM都单独训练,但是上一层的隐藏层直接作为下一层的可视层,逐层贪婪训练,这样就得到深度RBM网络。将神经元数逐层减少,最顶层的隐藏层即为降维后的数据编码,其输入为地震叠前数据训练样本,输出为训练完毕的深度RBM网络,降维后的地震叠前数据,深度RBM网络训练的具体流程如下:
S31:选取加入了噪声的数据集,将每个点的叠前数据Xi由矩阵转换成向量输入RBM网络,新的样本数据代表这个点处的所有信息,也是后续降维、分类的训练样本数据;
S32:将所有训练样本均值归一化到[0,1]范围内,这里数据预处理完成,保存为样本数据;
S33:因为临近位置样本总是相近的,直接输入易使RBM陷入局部最小或欠拟合,所以需要将所有训练样本随机排列。然后将数据分成多批,一方面使用批梯度上升可以降低算法复杂度,另一方面可以满足程序并发处理需求;
S34:开始训练第一层RBM,预处理过的训练数据作为可视层输入,训练完毕后得到RBM的各项参数与隐藏层输出,保存各项参数;
S35:开始训练第二层RBM,将上一层隐藏层的激活概率作为本层RBM可视层输入,训练完毕后得到RBM的各项参数与隐藏层输出,保存各项参数;
S36:按照步骤S35逐层训练RBM,直至算法训练完毕后,得到适应本训练样本的深度RBM网络;
S37:取S32中保存的样本数据作为输入,计算隐藏层的激活概率:这里计算激活概率可以采用现有技术,不再详细展开说明。
S38:传递上一层的隐藏层的激活概率,继续计算,直至计算到顶层节点激活概率,输出顶层隐藏层激活概率即为训练样本降维结果。
步骤S4.将降维后的数据作为模糊自组织神经网络的输入,进行无监督训练,得到最终分类结果。
对于模糊自组织神经网络的输出,每个样本对所有簇都有一个隶属度值,既可以选取最大隶属度值的簇作为归属类别,也可以对每个隶属度值做进一步探索。
由于实际地质结构往往很复杂,难以用一个精确的分类结果描述,或者说精确的描述本身就是不合理的。为了平滑不同簇之间的界限,本发明引入模糊理论将其与SOM结合,即模糊自组织神经网络(FSOM),这样能够对更好地解释地震相,更丰富的展示地震相信息。
对比标准自组织神经网络,模糊自组织神经网络只是将归属于某个神经元的簇改为对每个神经元的隶属度。所以算法的性能并不会降低,算法复杂度也没有改变。模糊自组织神经网络算法流程如下:
输入:数据集D=[z1,z2,...,zn],其中,zi是原始地震叠前数据Xi经过RBM降维后的向量;簇的数目k。
输出:输入层与输出层节点之间的权值向量Wj(t),j=1,2,...,p,其中,t是迭代次数,p是输出层节点个数;第i个样本对第j个簇的隶属度uij。
S41:初始化两层神经元之间的权值向量,随机选择输入数据的特征向量进行初始化;
S42:对任意样本初始化uij=1/k,表示初始每个样本对所有模糊簇的隶属度是相同的;
S43:依次将输入数据zi输入网络,计算dij=||zi-Wj||,表示zi到Wj的距离;
S44:求最大相似度对应的节点,即dij最小的节点,记为获胜神经元Wc(t);
S45:利用侧向反馈规则,对获胜神经元以及其周围的神经元Wj进行更新,具体的更新规则如下:
Wj(t+1)=Wj(t)+σ(t)hcj(t)[zi-Wc(t)] (2)
其中,σ(t)代表学习率,是一个调节参数用来保证神经网络最终的收敛,一般随着学习迭代次数的增多而下降,0<σ(t)<1,hcj(t)是侧向反馈距离,是获胜神经元与其他未获胜神经元之间的距离,这使得每一次迭代中,大部分神经元都能获得调整;
S46:更新所有样本对每个模糊簇的隶属度uij;
S47:i=i+1,继续输入数据,直到每个数据依次输入完毕,或者满足其他截止条件,例如当数据较少时,一般会考虑重复输入整个数据集合,直到达到一定的迭代次数。
FSOM使得每个样本不再固定被分为同类或不同类,给后续的地震相成图提供了更丰富的信息。
实验工区选取某实际地震工区部分工区,工区大小为600×700,6个方位角数据,选取目标段数据100个采样点。即样本数据量为420000,样本数据维度600。作为对比首先使用PCA将数据降维,然后SOM聚类成16类得到如图3所示的结果图;使用本发明的方法处理本工区,分成16类,结果图如图4所示。
可以明显看出,现有的方法生成的地震相图已可以大致反映出地下地质结构的分布,但精细度不够,没有清晰地描述复杂小裂缝,难以满足更高的勘探需求,不利于实际应用。
由于模糊聚类算法使得每个样本对每个簇都有一个隶属度值,对此提出一种的展示地震相图的方法。在每个类别上计算二维梯度,对于地震相图每个点赋予坐标(i,j),计算此位置对每个簇的隶属度记为uk(i,j),假设分为n类则k={1,2,...,n}。对于每个点计算其二维梯度如(3)式。
如图5所示,相对于普通地震相图隶属度梯度能够更好地展示地震相变化的边缘地带,对于大尺度裂缝的展示更加清晰。
Claims (3)
1.一种基于RBM地震叠前信号聚类方法,具体包括如下步骤:
步骤S1.去除地震叠前数据的噪声,保留结构特征,并将目标层位段数据组织成Xi=[x1;x2;...;xs],其中,x为单个方位角数据,维度为d,s为地震数据的方位角个数;
步骤S2.对n个位置点的数据U=[X1,X2,...,Xn]加入不同的噪声生成m份数据集U1,U2,...,Um;
步骤S3.按深度受限玻尔兹曼机网络对数据提取特征和维度规约;
步骤S4.将降维后的数据作为模糊自组织神经网络的输入,进行无监督训练,得到最终分类结果;模糊自组织神经网络的输出,每个样本对所有簇都有一个隶属度值,对此提出一种的展示地震相图的方法,在每个类别上计算二维梯度,对于地震相图每个点赋予坐标(i,j),计算此位置对每个簇的隶属度记为uk(i,j),假设分为n类则k={1,2,...,n},对于每个点计算其二维梯度如(3)式:
2.根据权利要求1所述的基于RBM地震叠前信号聚类方法,其特征在于,所述步骤S3的具体包括如下分步骤:
S31:选取加入了噪声的数据集,将每个点的叠前数据Xi由矩阵转换成向量输入RBM网络;
S32:将所有训练样本均值归一化到[0,1]范围内,保存为样本数据;
S33:将所有训练样本随机排列,然后将数据分成多批;
S34:开始训练第一层RBM,预处理过的训练数据作为可视层输入,训练完毕后得到RBM的各项参数与隐藏层输出,保存各项参数;
S35:开始训练第二层RBM,将上一层隐藏层的激活概率作为本层RBM可视层输入,训练完毕后得到RBM的各项参数与隐藏层输出,保存各项参数;
S36:按照步骤S35逐层训练RBM,直至算法训练完毕后,得到适应本训练样本的深度RBM网络;
S37:取S32中保存的样本数据作为输入,计算隐藏层的激活概率;
S38:传递上一层的隐藏层的激活概率,继续计算,直至计算到顶层节点激活概率,输出顶层隐藏层激活概率即为训练样本降维结果。
3.根据权利要求2所述的基于RBM地震叠前信号聚类方法,其特征在于,所述步骤S4的具体包括如下分步骤:
S41:初始化两层神经元之间的权值向量,随机选择输入数据的特征向量进行初始化;
S42:对任意样本初始化uik=1/k,表示初始每个样本对所有模糊簇的隶属度是相同的,k是簇的数目;
S43:依次将输入数据zi输入网络,计算dij=||zi-Wj||,表示zi到Wj的距离;
S44:求最大相似度对应的节点,即dij最小的节点,记为获胜神经元Wc(t);
S45:利用侧向反馈规则,对获胜神经元以及其周围的神经元Wj进行更新,具体的更新规则如下:
Wj(t+1)=Wj(t)+σ(t)hcj(t)[zi-Wc(t)] (2)
其中,σ(t)代表学习率,用来保证神经网络最终的收敛,0<σ(t)<1,hcj(t)是侧向反馈距离,是获胜神经元与其他未获胜神经元之间的距离;
S46:更新所有样本对每个模糊簇的隶属度uik;
S47:i=i+1,继续输入数据,直到每个数据依次输入完毕。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710728113.7A CN107688201B (zh) | 2017-08-23 | 2017-08-23 | 基于rbm地震叠前信号聚类方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710728113.7A CN107688201B (zh) | 2017-08-23 | 2017-08-23 | 基于rbm地震叠前信号聚类方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107688201A CN107688201A (zh) | 2018-02-13 |
CN107688201B true CN107688201B (zh) | 2019-12-31 |
Family
ID=61153639
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710728113.7A Active CN107688201B (zh) | 2017-08-23 | 2017-08-23 | 基于rbm地震叠前信号聚类方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107688201B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108732620B (zh) * | 2018-03-09 | 2019-12-27 | 山东科技大学 | 一种非监督与监督学习下的多波地震油气储层预测方法 |
CN108596630B (zh) * | 2018-04-28 | 2022-03-01 | 招商银行股份有限公司 | 基于深度学习的欺诈交易识别方法、***及存储介质 |
CN108629593B (zh) * | 2018-04-28 | 2022-03-01 | 招商银行股份有限公司 | 基于深度学习的欺诈交易识别方法、***及存储介质 |
GB2590555B (en) * | 2018-06-08 | 2022-06-08 | Schlumberger Technology Bv | Methods for characterizing and evaluating well integrity using unsupervised machine learning of acoustic data |
CN109272029B (zh) * | 2018-08-31 | 2022-02-22 | 电子科技大学 | 井控稀疏表征大规模谱聚类地震相划分方法 |
CN112444868B (zh) * | 2019-08-30 | 2024-04-09 | 中国石油化工股份有限公司 | 一种基于改进的k均值算法的地震相分析方法 |
AU2020101108A4 (en) * | 2020-06-24 | 2020-07-30 | Institute Of Geology And Geophysics, Chinese Academy Of Sciences | A method for geophysical observation information fusion |
CN112764110A (zh) * | 2020-07-09 | 2021-05-07 | 五季数据科技(北京)有限公司 | 一种基于限制波尔兹曼机特征编码的聚类地震相分析方法 |
CN114114410B (zh) * | 2020-08-28 | 2022-09-09 | 中国科学技术大学 | 一种识别异常地震信号的方法 |
CN114624763A (zh) * | 2020-12-10 | 2022-06-14 | 中国石油天然气股份有限公司 | 地震相图生成方法及装置 |
CN114137610B (zh) * | 2021-11-24 | 2023-05-02 | 同济大学 | 联合监督和非监督学习的低勘探区地层和岩性地震评价方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101430386A (zh) * | 2007-11-08 | 2009-05-13 | 中国石油天然气股份有限公司 | 一种地震多参数融合气藏检测方法 |
US8270733B2 (en) * | 2009-08-31 | 2012-09-18 | Behavioral Recognition Systems, Inc. | Identifying anomalous object types during classification |
CN103257360A (zh) * | 2013-04-25 | 2013-08-21 | 刘立峰 | 基于模糊c均值聚类的碳酸盐岩流体识别方法 |
CN104280771A (zh) * | 2014-10-27 | 2015-01-14 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种基于em算法的三维地震数据波形的半监督聚类方法 |
CN106097339A (zh) * | 2016-06-12 | 2016-11-09 | 成都甄识科技有限公司 | 一种基于模糊集的图像边缘检测方法 |
CN106226816A (zh) * | 2016-09-12 | 2016-12-14 | 电子科技大学 | 一种叠前地震信号波形分类方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6763303B2 (en) * | 2002-03-25 | 2004-07-13 | Council Of Scientific And Industrial Research | System for classifying seafloor roughness |
CN106405640B (zh) * | 2016-08-26 | 2018-07-10 | 中国矿业大学(北京) | 基于深度信念神经网络的微震信号到时自动拾取方法 |
-
2017
- 2017-08-23 CN CN201710728113.7A patent/CN107688201B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101430386A (zh) * | 2007-11-08 | 2009-05-13 | 中国石油天然气股份有限公司 | 一种地震多参数融合气藏检测方法 |
US8270733B2 (en) * | 2009-08-31 | 2012-09-18 | Behavioral Recognition Systems, Inc. | Identifying anomalous object types during classification |
CN103257360A (zh) * | 2013-04-25 | 2013-08-21 | 刘立峰 | 基于模糊c均值聚类的碳酸盐岩流体识别方法 |
CN104280771A (zh) * | 2014-10-27 | 2015-01-14 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种基于em算法的三维地震数据波形的半监督聚类方法 |
CN106097339A (zh) * | 2016-06-12 | 2016-11-09 | 成都甄识科技有限公司 | 一种基于模糊集的图像边缘检测方法 |
CN106226816A (zh) * | 2016-09-12 | 2016-12-14 | 电子科技大学 | 一种叠前地震信号波形分类方法 |
Non-Patent Citations (1)
Title |
---|
模糊综合评判法评价尕斯库勒N1-N2油藏北区块储层;胡利华;《青海石油》;20120930;第30卷(第3期);第56-61页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107688201A (zh) | 2018-02-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107688201B (zh) | 基于rbm地震叠前信号聚类方法 | |
Zhu et al. | Intelligent logging lithological interpretation with convolution neural networks | |
CN109709603B (zh) | 地震层位识别与追踪方法、*** | |
AlRegib et al. | Subsurface structure analysis using computational interpretation and learning: A visual signal processing perspective | |
AU2013378058B2 (en) | Detecting subsurface structures | |
Qi et al. | Seismic attribute selection for machine-learning-based facies analysis | |
CN110609320B (zh) | 一种基于多尺度特征融合的叠前地震反射模式识别方法 | |
CN109272029B (zh) | 井控稀疏表征大规模谱聚类地震相划分方法 | |
EA004551B1 (ru) | Способ интерпретации сейсмических фаций с использованием текстурного анализа и нейронных сетей | |
CN109581487B (zh) | 一种沉积体波形空间结构分析与体控反演建模方法 | |
CN112731522B (zh) | 地震地层智能识别方法、装置、设备及存储介质 | |
CN117546051A (zh) | 使用s波速度模型和机器学习进行地震成像的方法和*** | |
Xiao et al. | Conditioning of deep-learning surrogate models to image data with application to reservoir characterization | |
WO2022187685A1 (en) | Method and system for faster seismic imaging using machine learning | |
WO2017152123A1 (en) | Context based bounded hydrocarbon formation identification | |
Guo et al. | Reconstruction of missing shear wave slowness in shale gas reservoir based on deep learning algorithm | |
CN108226997A (zh) | 一种基于叠前地震数据的地震相划分方法 | |
CN115639605B (zh) | 基于深度学习的高分辨率断层的自动识别方法和装置 | |
CN116068644B (zh) | 一种利用生成对抗网络提升地震数据分辨率和降噪的方法 | |
CN111273346B (zh) | 去除沉积背景的方法、装置、计算机设备及可读存储介质 | |
CN106997060A (zh) | 一种基于Shearlet‑fastICA的地震多属性融合方法 | |
Koochak et al. | A variability aware GAN for improving spatial representativeness of discrete geobodies | |
CN113419278B (zh) | 一种基于状态空间模型与支持向量回归的井震联合多目标同时反演方法 | |
Hong et al. | A DTW distance-based seismic waveform clustering method for layers of varying thickness | |
Liu et al. | Robust prestack seismic facies analysis using shearlet transform-based deep learning |
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 |