CN110806602B - 基于深度学习的智能化地震数据随机噪声压制方法 - Google Patents
基于深度学习的智能化地震数据随机噪声压制方法 Download PDFInfo
- Publication number
- CN110806602B CN110806602B CN201810886433.XA CN201810886433A CN110806602B CN 110806602 B CN110806602 B CN 110806602B CN 201810886433 A CN201810886433 A CN 201810886433A CN 110806602 B CN110806602 B CN 110806602B
- Authority
- CN
- China
- Prior art keywords
- training
- model
- self
- seismic
- seismic data
- 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
- 230000001629 suppression Effects 0.000 title claims abstract description 36
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000013135 deep learning Methods 0.000 title claims abstract description 23
- 238000012549 training Methods 0.000 claims abstract description 132
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 36
- 238000013527 convolutional neural network Methods 0.000 claims abstract description 36
- 230000006870 function Effects 0.000 claims description 27
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000013528 artificial neural network Methods 0.000 claims description 5
- 230000005284 excitation Effects 0.000 claims description 5
- 230000004913 activation Effects 0.000 claims description 4
- 230000001419 dependent effect Effects 0.000 claims description 4
- 238000005422 blasting Methods 0.000 claims description 3
- 230000001537 neural effect Effects 0.000 claims 1
- 238000005516 engineering process Methods 0.000 description 21
- 238000007418 data mining Methods 0.000 description 10
- 238000012545 processing Methods 0.000 description 10
- 230000000694 effects Effects 0.000 description 9
- 238000005070 sampling Methods 0.000 description 6
- 238000013473 artificial intelligence Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 230000009466 transformation Effects 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000013499 data model Methods 0.000 description 2
- 238000000354 decomposition reaction Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005065 mining Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 230000003213 activating effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 238000003062 neural network model Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 229910052704 radon Inorganic materials 0.000 description 1
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Images
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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/30—Noise handling
- G01V2210/32—Noise reduction
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
本发明提供了一种基于深度学习的智能化地震数据随机噪声压制方法,包括:生成用于地震资料正演的多个速度模型;利用声波方程,针对每个速度模型正演得到相对应的地震单炮记录作为训练输出标签,并加入噪声形成训练数据输入集;以及利用卷积神经网络算法并结合所述训练数据输入集,构建自编码训练模型、完成模型训练,并将其用于地震资料随机噪声压制。
Description
技术领域
本发明涉及地球物理勘探技术领域。具体地,涉及一种利用深度学习算法进行地震资料去噪的方法,可应用于石油地球物理勘探中的地震资料去噪阶段。
背景技术
地震信号在产生、传输和记录过程中,经常会受到各种噪声的干扰,噪声问题一直是地震资料处理的主要问题。随着勘探程度的不断加深,“高信噪比、高分辨率、高保真度”已经成为地震资料处理的目标,实现“三高”的地震处理目标的关键在于高质量的去噪。如何更好地压制噪声,保护有效信号的能量和频率不受损失是目前地震资料去噪技术的关键。
目前去噪的主流方法依旧是通过多域变换,将时间域信号变换到不同域,利用噪声与有效信号的性质差异进行去噪,例如基于傅里叶变换的F-K滤波,它通过噪声和有效信号在频率上的差异进行噪声压制;基于Radon变换的去噪技术,它通过将信号变换到τ-p域实现噪声和有效信号的分离;基于希尔伯特-黄变换的噪声压制技术,该技术通过二维经验模式分解对信号进行多尺度分解,实现噪声与有效信号的分离,再通过数据重构达到去噪的效果。这些技术虽然能够去除大部分噪声,但在一定程度上对有效信号造成破坏,且算法复杂,甚至会产生虚假有效信号。另外,这些方法的另一主要问题在于引入了处理人员的主观干预,可能将有效信号当噪声,或将噪声识别为有效信号,从而影响去噪效果。
随着地震采集技术的进步和计算机技术的飞速发展,勘探地震领域中积累了大量的数据,对这些数据进行智能化分析以发掘数据中蕴含的有用信息,成为时代发展的需求,尤其是利用深度学习开展智能化地震数据去噪技术研究。基于数据挖掘及深度学习的人工智能技术在新兴行业领域和其他传统行业领域已经取得较大进展,无论是算法层面还是技术层面都已相对成熟。在地震勘探领域,数据挖掘技术和基于深度学习的人工智能技术正在大规模开展广泛研究,但针对地震勘探资料的处理阶段目前尚未找到有效切入点。
发明内容
面向地震资料噪声信息的数据挖掘目的在于利用数据挖掘和地震数据处理中的方法,自动提取叠前地震道集的特征信息,并进行抽象描述以达到最小化模拟数据与实际数据的特征信息差别,达到去除噪声、识别有效信号的目的,实现基于深度学习的智能化去噪技术。
针对地震数据处理中占用了大量人力、时间的噪声压制问题,将目前功能强大的人工智能技术引入到勘探地震中,开展基于数据挖掘和深度学习的智能化随机噪声压制技术研究,充分融合地球物理专家经验,最大力度识别及去除噪声,不伤害有效信号,形成针对勘探地震的数据挖掘技术和基于深度学习的智能化去噪技术,有效提高地震资料信噪比,降低地震勘探成本,缩短勘探周期,减少人工主观干预,为地震勘探中的噪声问题提供新的解决方案。
根据本发明的一方面,提供了一种基于深度学习的智能化地震数据随机噪声压制方法,包括:
生成用于地震资料正演的多个速度模型;
利用声波方程,针对每个速度模型正演得到相对应的地震单炮记录作为训练输出标签,并加入噪声形成训练数据输入集;以及
利用卷积神经网络算法并结合所述训练数据输入集,构建自编码训练模型、完成模型训练,并将其用于地震资料随机噪声压制。
在实施例中,生成用于地震资料正演的多个速度模型,包括:
每个速度模型中的地层层位数量、地层的倾斜角度,以及各层的层内速度均随机生成。
具体地,在优选实施例中,每个速度模型中的地层层位数量在2-12之间随机生成,地层的倾斜角度在0-80度之间随机生成,以及各层的层内速度在500m/s-7000m/s之间随机生成。
在实施例中,利用声波方程,针对每个速度模型正演得到相对应的地震单炮记录作为训练输出标签,并加入噪声形成训练数据输入集,包括:
利用二维声波方程,模拟地震波在每个所述速度模型中的传播状态,即激发方式采用中间放炮,两侧接收的观测系统,正演得到相对应的地震单炮记录作为训练输出标签,并加入噪声形成训练数据输入集,使得每个速度模型可以正演模拟得到一个单炮记录。
具体地,在优选实施例中,该观测系统可以左右两侧各放置64个检波器,共128道接收,纵向采样点为1000,采样间隔为4ms,正演得到相对应的地震单炮记录作为训练输出标签,并加入噪声形成训练数据输入集,使得每个速度模型可以正演模拟得到一个单炮记录。
在实施例中,将正演得到的地震单炮记录加入预定信噪比的随机噪声,以作为所述训练数据输入集。
在优选实施例中,所述预定信噪比为10。
在实施例中,利用卷积神经网络算法,构建自编码训练模型、完成模型训练,并将其用于地震资料随机噪声压制,包括:
利用如下的三层卷积神经网络算法公式,提取地震数据特征,
其中,上述公式表示第j层卷积神经网络算法,hj W,b(xj)表示第j层经过卷积神经网络计算后的结果,其中j=1,2,3,x1表示所述训练数据输入集;x2即为h1 W,b(x1);x3表示h2 W,b(x2);W表示权重,即卷积核;b为数据扰动项,i表示第i个单元,n表示单元总数;f表示激活函数,在此采用ReLU函数,其定义为:
其中m表示函数自变量,f表示函数因变量。
在实施例中,利用卷积神经网络算法并结合所述训练数据输入集,构建自编码训练模型、完成模型训练,并将其用于地震资料随机噪声压制,还包括:
在经过三层卷积核为预定尺寸的卷积神经网络计算之后,利用最邻近插值算法,恢复地震数据原始尺寸。
在实施例中,利用卷积神经网络算法并结合所述训练数据输入集,构建自编码训练模型、完成模型训练,并将其用于地震资料随机噪声压制,还包括:
利用如下的最小均方差(MSE)函数,确定所述自编码训练模型的输出与输入之间的残差,
其中y'表示所述训练数据输入集经过所述自编码训练模型计算得到的输出值,y表示所述训练输出标签。
在实施例中,利用卷积神经网络算法并结合所述训练数据输入集,构建自编码训练模型、完成模型训练,并将其用于地震资料随机噪声压制,还包括:
通过将残差E对权重W和扰动项b求偏导,得到所述自编码训练模型中每层各单元的权重W和扰动项b:
其中h是训练步长,k是反馈迭代次数。
在实施例中,利用卷积神经网络算法并结合所述训练数据输入集,构建自编码训练模型、完成模型训练,并将其用于地震资料随机噪声压制,还包括:
利用多个所述地震单炮记录,完成所述自编码训练模型的训练,以用于地震资料的随机噪声压制。
根据本发明的另一方面,还提供了一种存储介质,其中存储有计算机可执行的程序,所述程序在被执行时适于实施基于深度学习的智能化地震数据随机噪声压制方法,所述方法包括:
生成用于地震资料正演的多个速度模型;
利用声波方程,针对每个速度模型正演得到相对应的地震单炮记录作为训练输出标签,并加入噪声形成训练数据输入集;以及
利用卷积神经网络算法并结合所述训练数据输入集,构建自编码训练模型、完成模型训练,并将其用于地震资料随机噪声压制。
本发明的其它特征和优点将在随后的说明书中阐述,并且部分地从说明书中变得显而易见,或者通过实施本发明而了解。本发明的目的和其他优点可通过在说明书、权利要求书以及说明书附图中所特别指出的结构来实现和获得。
附图说明
附图用来提供对本发明的进一步理解,并且构成说明书的一部分,与本发明的实施例共同用于解释本发明,并不构成对本发明的限制。在附图中:
图1示出根据本发明实施例的基于三层卷积神经网络的自编码器模型的示意图;
图2示出根据本发明实施例的随机生成的某个速度场;
图3示出根据本发明实施例的图1中的速度场对应生成的单炮地震记录;
图4示出根据本发明实施例的加入信噪比为10随机噪声后的单炮记录;
图5示出根据本发明实施例的自编码训练模型在训练阶段的收敛曲线,横坐标为迭代次数,纵坐标为训练误差(%);
图6示出根据本发明实施例的利用自编码训练模型进行噪声压制后的结果;
图7示出根据本发明实施例的无噪数据、含噪数据及去噪数据的频谱对比,横坐标表示频率(Hz),纵坐标表示振幅;
图8示出根据本发明实施例的无噪数据与去噪数据的单道对比,其中横坐标表示时间(s),纵坐标表示振幅;
图9示出根据本发明实施例的无噪数据(a)、含噪数据(b)及去噪数据(c)的局部细节对比。
具体实施方式
以下将结合附图及实施例来详细说明本发明的实施方式,借此对本发明如何应用技术手段来解决技术问题,并达成技术效果的实现过程能充分理解并据以实施。需要说明的是,只要不构成冲突,本发明中的各个实施例以及各实施例中的各个特征可以相互结合,所形成的技术方案均在本发明的保护范围之内。
同时,在以下说明中,出于解释的目的而阐述了许多具体细节,以提供对本发明实施例的彻底理解。然而,对本领域的技术人员来说显而易见的是,本发明可以不必用这里的具体细节或者所描述的特定方式来实施。
如上所述,本发明针对地震数据处理中占用了大量人力、时间的噪声压制问题,将功能强大的人工智能技术引入到勘探地震中,开展基于数据挖掘和深度学习的智能化随机噪声压制技术研究,充分融合地球物理专家经验,最大力度识别及去除噪声,不伤害有效信号,形成针对勘探地震的数据挖掘技术和基于深度学习的智能化去噪技术。
总体上,本发明通过:首先生成大量的速度模型用于地震资料正演;然后采用声波方程,针对每个速度模型正演得到相对应的地震单炮记录作为训练输出标签,并加入噪声形成数据训练输入集;最后采用卷积神经网络算法并结合所述训练数据输入集,构建自编码训练模型,实现基于数据挖掘的无监督训练模块,通过最小化输入-输出残差完成模型训练,并将其用于地震资料随机噪声压制。
本发明的方法能够有效提高地震资料信噪比,降低地震勘探成本,缩短勘探周期,减少人工主观干预,为地震勘探中的噪声问题提供新的解决方案。
以下,具体描述本发明的基于深度学习的智能化地震数据随机噪声压制方法的技术方案:
1)随机生成多个例如10000个速度模型,其中,各速度模型中地层层位数量在2-12之间随机生成、地层的倾斜角度在0-80度之间随机生成,各层的层内速度在500m/s-7000m/s之间随机生成,以保证样本的随机性;以上数值及其范围是非限制性的,可以根据实际需要选择不同数值及其范围;
2)在实施例中,设定速度模型数据体大小为128x1000,该尺寸可根据实际资料的需求调整;
3)采用二维声波方程,针对每个速度模型对应生成一个单炮地震记录,激发方式采用中间放炮,两边接收的观测系统,左右两边各放置64个检波器,共128道接收,纵向采样点为1000,采样间隔为4ms,该观测系统可根据实际资料调整;
4)对正演得到的单炮地震记录加入预定信噪比,例如10的随机噪声作为训练数据输入集,将正演得到的原始炮集作为训练输出标签;
5)利用三层卷积神经网络提取地震数据特征,卷积神经网络算法公式如下:
其中上述公式代表第j层卷积神经网络算法,hj W,b(xj)表示第j层经过卷积神经网络计算后的结果,进入下一层后hj W,b(xj)将作为输入x计算hj+1 W,b(xj+1),这样建立三层卷积神经网络算法用于数据特征的提取,即在实施例中j=1,2,3。x1表示输入地震数据单元,即第四步加噪后的训练数据输入集;x2即为h1 W,b(x1);x3表示h2 W,b(x2)。W表示权重,即卷积核,本发明中采用2x2的卷积核。b为数据扰动项,i表示第i个单元,n表示单元总数。f表示激活函数,这里采用ReLU函数,其定义为:
其中m表示函数自变量,f表示函数因变量,激活函数的目的是为了增加神经网络模型的表达能力,提高模型的复杂度。
6)经过三层卷积核为2x2的卷积神经网络计算之后,128x1000的数据将被压缩为16x125,这种操作可以保持原有数据特征的同时减少权系数W和数据扰动项b的数量,提高训练效率,最后再利用最邻近插值算法(resize_nearest_neighbor)恢复数据原始尺寸,具体的做法是:为了保证数据特征不出现模糊现象,每次插值纵、横向仅扩大一倍,那么16x125的数据经过一次插值之后就变味32x250,这样经过三次插值即可实现数据模型大小的复原。
7)更新权值W与扰动项b:通过上面的模型搭建可以看出,实际上模型训练是想通过最小化输出与期望标签的残差来确定权值W和扰动项b,即权值W和扰动项b才是本发明中的未知数;要更新权值和扰动项,首先要确定价值函数(cost function),通常采用的MSE(最小均方差)函数,即:期望标签值(即第四步中没有加入噪声的原始炮集)减去实际输出值(利用第四步含噪数据作为输入经过模型计算后的值)后求均方差。其中。公式如下:
其中y'表示第4步的训练数据输入集经过所搭建模型计算得到的输出值,y表示训练输出标签,即第四步中的原始炮集,i表示第i个单元序列,n表示数据单元总数。
8)确定价值函数以后,通过残差E的大小实现数据反馈,从而对W和b不断调整,直到E足够小为止,停止计算。残差E足够小实际上代表了含噪数据经过模型计算后达到了去噪的效果,所得结果若与无噪数据非常接近,这样训练后的模型就可以实现去噪的效果。具体的反馈算法如下:通过E对W和b求偏导,得到自编码训练模型中每层各单元的权重W和扰动项b:
其中h是训练步长,k是反馈迭代次数。
9)通过上述步骤搭建智能化去噪模型,利用10000个正演记录进行训练,训练结束后所得模型即为针对地震资料去噪的自编码训练模型。
图1示出根据本发明实施例的基于三层卷积神经网络的自编码器模型的示意图。如图1所示,数据体大小为128x1000的输入层,在经过三层卷积核为2x2的卷积神经网络计算之后,128x1000的数据将被压缩为16x125,再经过三次插值,即可实现数据模型大小的复原。通过构建三层卷积神经网络的自编码训练模型,最小化目标单元和输入单元的差别来进行优化,即让输出层(利用含噪数据作为输入通过模型计算后的结果)尽可能地去学习有效波信号的特征信息,从而达到去噪目的。
以下,结合图2-图9描述根据本发明实施例的基于深度学习的智能化地震数据随机噪声压制方法的正确性和有效性。
图2示出根据本发明实施例的随机生成的某个速度场。如图2所示,这是随机生成10000个速度模型之一,其中速度模型中的层位数量、倾角以及层内速度均为随机生成。
图3示出根据本发明实施例的图1中的速度场对应生成的单炮地震记录。如图3所示,采用声波方程正演,对每个速度模型生成单个地震单炮记录,激发方式采用中间放炮,两边接收的观测系统,左右两边各放置64个检波器,共128道接收,纵向采样点为1000,采样间隔为4ms。
图4示出根据本发明实施例的加入信噪比为10随机噪声后的单炮记录。如图4所示,对生成的单炮地震记录加入信噪比为10的随机噪声。
图5示出根据本发明实施例的自编码训练模型在训练阶段的收敛曲线,横坐标为迭代次数,纵坐标为训练误差(%)。如图5所示,利用10000个样本组成数据训练集对模型进行迭代训练更新,该图是训练对应的收敛曲线,可以看出,该模型收敛迅速。
图6示出根据本发明实施例的利用自编码训练模型进行噪声压制后的结果。如图6所示,这是利用训练完成的卷积自编码器模型得到的去噪效果图,可以看出,噪声得到了较好的压制,有效信号得以保存。
图7示出根据本发明实施例的无噪数据、含噪数据及去噪数据的频谱对比,横坐标表示频率(Hz),纵坐标表示振幅。如图7所示,分别对同一个样本的无噪数据、含噪数据及去噪数据进行频谱分析,验证去噪效果,从结果中可以看出,去噪前后数据频谱基本一致,也证明了该方法去噪能力。
图8示出根据本发明实施例的无噪数据与去噪数据的单道对比,其中横坐标表示时间(s),纵坐标表示振幅。如图8所示,提取出去噪前后的单道地震记录(此处以第一道为例)进行对比,通过单道对比,可以更直观的看出该模型的去噪效果较好,数据重构能力强。。
图9示出根据本发明实施例的无噪数据(a)、含噪数据(b)及去噪数据(c)的局部细节对比。如图9所示,对去噪前后的数据进行局部放大,这里以浅层近偏移距数据为例,通过对比验证其去噪效果。其中图9(a)表示无噪数据,图9(b)表示含噪数据,图9(c)表示去噪后数据。
综上所述,本发明提供了一种基于深度学习的智能化地震数据随机噪声压制方法。与现有技术相比,本发明的有益效果是:利用数据挖掘技术,自动挖掘有效波与噪声的特征信息,形成适合深度学习的规则输入数据,搭建卷积自编码器模型,进而开展基于深度学习的智能化地震数据去噪技术研究,研发针对地震数据噪声压制损失函数,建立地震资料去噪的人工智能神经网络,通过大规模数值训练优化神经网络,形成智能化噪声压制技术,实现保幅去噪,在提高信噪比的同时保证有效信号不受损失,特别是保护低频信号方面。
1)为了确保训练数据集的随机性,生成的单炮记录对应的速度模型在层位数量、层速度及层位倾角均为随机产生;
2)将深度学习引入地震资料去噪处理环节,通过卷积神经网络算法实现地震资料的特征挖掘,自动区分有效波信号及噪声;
3)通过卷积神经网络算法,对原始地震资料进行低秩降维处理,得到适于机器学习的规范化地震数据;
4)搭建三层卷积神经网络的自编码器模型,通过去最小化目标单元和输入单元的差别来进行优化,即让输出层(利用含噪数据作为输入通过模型计算后的结果)尽可能地去学习有效波信号的特征信息,达到去噪目的。
根据本发明的另一方面,还提供了一种存储介质,其中存储有计算机可执行的程序,所述程序在被执行时适于实施基于深度学习的智能化地震数据随机噪声压制方法,所述方法包括:生成用于地震资料正演的多个速度模型;利用声波方程,针对每个速度模型正演得到相对应的地震单炮记录作为训练输出标签,并加入噪声形成训练数据输入集;以及利用卷积神经网络算法并结合所述训练数据输入集,构建自编码训练模型、完成模型训练,并将其用于地震资料随机噪声压制。
应该理解的是,本发明所公开的实施例不限于这里所公开的特定处理步骤或材料,而应当延伸到相关领域的普通技术人员所理解的这些特征的等同替代。还应当理解的是,在此使用的术语仅用于描述特定实施例的目的,而并不意味着限制。
说明书中提到的“实施例”意指结合实施例描述的特定特征、或特性包括在本发明的至少一个实施例中。因此,说明书通篇各个地方出现的短语“实施例”并不一定均指同一个实施例。
本领域的技术人员应该明白,结合本文中所公开的实施例描述的各示例的单元及算法步骤,能够以电子硬件、计算机软件或者二者的结合来实现,为了清楚地说明硬件和软件的可互换性,在上述说明中已经按照功能一般性地描述了各示例的组成及步骤。这些功能究竟以硬件还是软件方式来执行,取决于技术方案的特定应用和设计约束条件。本领域的技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能,但是这种实现不应认为超出本发明的范围。
结合本文中所公开的实施例描述的方法或算法的步骤可以直接用硬件、处理器执行的软件模块,或者二者的结合来实施。软件模块可以置于随机存储器(RAM)、内存、只读存储器(ROM)、电可编程ROM、电可擦除可编程ROM、寄存器、硬盘、可移动磁盘、CD-ROM或技术领域内所公知的任意其它形式的存储介质中。
虽然本发明所公开的实施方式如上,但所述的内容只是为了便于理解本发明而采用的实施方式,并非用以限定本发明。任何本发明所属技术领域内的技术人员,在不脱离本发明所公开的精神和范围的前提下,可以在实施的形式上及细节上作任何的修改与变化,但本发明的保护范围,仍须以所附的权利要求书所界定的范围为准。
Claims (5)
1.一种基于深度学习的智能化地震数据随机噪声压制方法,包括:
生成用于地震资料正演的多个速度模型;
利用二维声波方程,模拟地震波在每个所述速度模型中的传播状态,即激发方式采用中间放炮,两侧接收的观测系统,正演得到相对应的地震单炮记录作为训练输出标签,并加入噪声形成训练数据输入集,使得每个速度模型可以正演模拟得到一个单炮记录;
利用卷积神经网络算法并结合所述训练数据输入集,构建自编码训练模型、完成模型训练,并将其用于地震资料随机噪声压制;
所述利用卷积神经网络算法并结合所述训练数据输入集,构建自编码训练模型、完成模型训练,并将其用于地震资料随机噪声压制,包括:
利用如下的三层卷积神经网络算法公式,提取地震数据特征,
其中,上述公式表示第j层卷积神经网络算法,hj W,b(xj)表示第j层经过卷积神经网络计算后的结果,其中j=1,2,3,x1表示所述训练数据输入集;x2即为h1 W,b(x1);x3表示h2 W,b(x2);W表示权重,即卷积核;b为数据扰动项,i表示第i个单元,n表示单元总数;f表示激活函数,在此采用ReLU函数,其定义为:
其中m表示函数自变量,f表示函数因变量;
在经过三层卷积核为预定尺寸的卷积神经网络计算之后,利用最邻近插值算法,恢复地震数据原始尺寸;
利用如下的最小均方差(MSE)函数,确定所述自编码训练模型的输出与输入之间的残差,
其中y'表示所述训练数据输入集经过所述自编码训练模型计算得到的输出值,y表示所述训练输出标签;
通过将残差E对权重W和扰动项b求偏导,得到所述自编码训练模型中每层各单元的权重W和扰动项b:
δh=hj W,b(xj)-hj-1 W,b(xj-1),
δy=y'-y
其中δh是训练步长,k是反馈迭代次数。
2.根据权利要求1所述的方法,其中,生成用于地震资料正演的多个速度模型,包括:
每个速度模型中的地层层位数量、地层的倾斜角度,以及各层的层内速度均随机生成。
3.根据权利要求1所述的方法,其中,将正演得到的地震单炮记录加入预定信噪比的随机噪声,以作为所述训练数据输入集。
4.根据权利要求1所述的方法,其中,利用卷积神经网络算法并结合所述训练数据输入集,构建自编码训练模型、完成模型训练,并将其用于地震资料随机噪声压制,还包括:
利用多个所述地震单炮记录,完成所述自编码训练模型的训练,以用于地震资料的随机噪声压制。
5.一种存储介质,其中存储有计算机可执行的程序,所述程序在被执行时适于实施基于深度学习的智能化地震数据随机噪声压制方法,所述方法包括:
生成用于地震资料正演的多个速度模型;
利用二维声波方程,模拟地震波在每个所述速度模型中的传播状态,即激发方式采用中间放炮,两侧接收的观测系统,正演得到相对应的地震单炮记录作为训练输出标签,并加入噪声形成训练数据输入集,使得每个速度模型可以正演模拟得到一个单炮记录;
利用卷积神经网络算法并结合所述训练数据输入集,构建自编码训练模型、完成模型训练,并将其用于地震资料随机噪声压制;
所述利用卷积神经网络算法并结合所述训练数据输入集,构建自编码训练模型、完成模型训练,并将其用于地震资料随机噪声压制,包括:
利用如下的三层卷积神经网络算法公式,提取地震数据特征,
其中,上述公式表示第j层卷积神经网络算法,hj W,b(xj)表示第j层经过卷积神经网络计算后的结果,其中j=1,2,3,x1表示所述训练数据输入集;x2即为h1 W,b(x1);x3表示h2 W,b(x2);W表示权重,即卷积核;b为数据扰动项,i表示第i个单元,n表示单元总数;f表示激活函数,在此采用ReLU函数,其定义为:
其中m表示函数自变量,f表示函数因变量;
在经过三层卷积核为预定尺寸的卷积神经网络计算之后,利用最邻近插值算法,恢复地震数据原始尺寸;
利用如下的最小均方差(MSE)函数,确定所述自编码训练模型的输出与输入之间的残差,
其中y'表示所述训练数据输入集经过所述自编码训练模型计算得到的输出值,y表示所述训练输出标签;
通过将残差E对权重W和扰动项b求偏导,得到所述自编码训练模型中每层各单元的权重W和扰动项b:
δh=hj W,b(xj)-hj-1 W,b(xj-1),
δy=y'-y
其中δh是训练步长,k是反馈迭代次数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810886433.XA CN110806602B (zh) | 2018-08-06 | 2018-08-06 | 基于深度学习的智能化地震数据随机噪声压制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810886433.XA CN110806602B (zh) | 2018-08-06 | 2018-08-06 | 基于深度学习的智能化地震数据随机噪声压制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110806602A CN110806602A (zh) | 2020-02-18 |
CN110806602B true CN110806602B (zh) | 2022-04-22 |
Family
ID=69487254
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810886433.XA Active CN110806602B (zh) | 2018-08-06 | 2018-08-06 | 基于深度学习的智能化地震数据随机噪声压制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110806602B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111580161A (zh) * | 2020-05-21 | 2020-08-25 | 长江大学 | 基于多尺度卷积自编码神经网络的地震随机噪音压制方法 |
AU2020101108A4 (en) * | 2020-06-24 | 2020-07-30 | Institute Of Geology And Geophysics, Chinese Academy Of Sciences | A method for geophysical observation information fusion |
CN113900146A (zh) * | 2020-07-06 | 2022-01-07 | 中国石油天然气股份有限公司 | 面波压制方法及系统 |
CN114167500A (zh) * | 2020-09-10 | 2022-03-11 | 中国石油天然气股份有限公司 | 基于卷积神经网络的邻炮干扰分离方法及装置 |
CN112419166B (zh) * | 2020-09-24 | 2024-01-05 | 南京晓庄学院 | 一种基于局部区域分割与scn结合的图像去雾方法 |
CN112946749B (zh) * | 2021-02-05 | 2022-05-20 | 北京大学 | 基于数据增广训练深度神经网络压制地震多次波的方法 |
CN113743193B (zh) * | 2021-07-22 | 2022-12-09 | 西安交通大学 | 一种叠前地震资料线性干扰压制方法及系统 |
CN113495291A (zh) * | 2021-07-27 | 2021-10-12 | 成都爱为贝思科技有限公司 | 基于深度学习的偏前道集保幅性智能量化评价方法 |
CN114114421B (zh) * | 2021-11-05 | 2023-09-29 | 中国石油大学(华东) | 基于深度学习的导向自学习地震数据去噪方法及装置 |
CN114035226B (zh) * | 2021-11-12 | 2023-10-20 | 东北石油大学 | 基于噪声语义相关性深度迁移去噪自编码网络的微地震资料噪声压制方法及系统 |
CN115221912A (zh) * | 2022-01-05 | 2022-10-21 | 中国地质大学(北京) | 一种基于自监督学习的地震数据随机噪音压制方法及系统 |
CN116561481B (zh) * | 2023-05-10 | 2023-12-19 | 广东海洋大学 | 一种针对avo效应的速度分析方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106291701A (zh) * | 2016-10-13 | 2017-01-04 | 成都理工大学 | 储层检测方法及装置 |
WO2018106805A1 (en) * | 2016-12-09 | 2018-06-14 | William Marsh Rice University | Signal recovery via deep convolutional networks |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106204468B (zh) * | 2016-06-27 | 2019-04-26 | 深圳市未来媒体技术研究院 | 一种基于ReLU卷积神经网络的图像去噪方法 |
CN106873036A (zh) * | 2017-04-28 | 2017-06-20 | 中国石油集团川庆钻探工程有限公司 | 一种基于井震结合的去噪方法 |
CN108198154B (zh) * | 2018-03-19 | 2020-06-26 | 中山大学 | 图像去噪方法、装置、设备及存储介质 |
-
2018
- 2018-08-06 CN CN201810886433.XA patent/CN110806602B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106291701A (zh) * | 2016-10-13 | 2017-01-04 | 成都理工大学 | 储层检测方法及装置 |
WO2018106805A1 (en) * | 2016-12-09 | 2018-06-14 | William Marsh Rice University | Signal recovery via deep convolutional networks |
Non-Patent Citations (5)
Title |
---|
Beyond a Gaussian Denoiser: Residual Learning of Deep CNN for Image Denoising;Kai Zhang et al.;《IEEE TRANSACTIONS ON IMAGE PROCESSING》;20170731;第26卷(第7期);第3142-3155页 * |
Deep Feature Learning for Medical Image Analysis with Convolutional Autoencoder Neural Network;Min Chen et al.;《2016 IEEE》;20161231;第1-10页 * |
Random Noise Attenuation Using Convolutional Neural Networks;J. Liu et al;《80th RAGR Conference & Rxhibition 2018》;20180630;第2-3页 * |
卷积神经网络研究综述;李彦冬等;《计算机应用》;20161231;第2508-2515页 * |
基于深度学习的地震噪声压制;于四伟等;《CPS/SEG北京2018国际地球物理会议暨展览电子论文集》;20180430;第437-439页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110806602A (zh) | 2020-02-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110806602B (zh) | 基于深度学习的智能化地震数据随机噪声压制方法 | |
CN111190227B (zh) | 基于残差卷积生成对抗模型的低信噪比地震数据消噪方法 | |
CN105467444A (zh) | 一种弹性波全波形反演方法及装置 | |
CN110658557B (zh) | 基于生成对抗网络的地震数据面波压制方法 | |
CN109633753B (zh) | 地表噪声压制方法及装置 | |
CN105700020A (zh) | 一种地震数据随机噪声压制方法及装置 | |
Meng et al. | Self-supervised learning for seismic data reconstruction and denoising | |
Liu et al. | Unsupervised deep learning for random noise attenuation of seismic data | |
CN113962244A (zh) | 瑞雷波地震数据噪声去除方法、存储介质及电子设备 | |
CN115905805A (zh) | 一种基于全局信息判别gan的das数据多尺度噪声消减方法 | |
Li et al. | An attention‐guided convolution neural network for denoising of distributed acoustic sensing–vertical seismic profile data | |
CN112926232A (zh) | 一种基于分层融合的地震低频分量恢复方法 | |
US11460594B1 (en) | Mitigating wireless channel impairments in seismic data transmission using deep neural networks | |
CN107255833B (zh) | 地震叠加剖面的确定方法和装置 | |
Fang et al. | Seismic Random Noise Suppression Model based on Downsampling and Super-resolution | |
CN112578439B (zh) | 一种基于空间约束的地震反演方法 | |
CN113009564B (zh) | 地震数据处理方法和装置 | |
CN108037533A (zh) | 一种基于块重组期望对数似然的地震勘探噪声压制方法 | |
CN113567129A (zh) | 一种列车轴承振动信号基于ceemd的降噪方法 | |
CN110008633B (zh) | 基于人工智能深度神经网络的公路噪音压制方法及系统 | |
Liao et al. | A twice denoising autoencoder framework for random seismic noise attenuation | |
Cheng et al. | An effective self-supervised learning method for various seismic noise attenuation | |
CN112363217A (zh) | 一种地震数据随机噪声压制方法及系统 | |
CN114460633B (zh) | 基于加权框架变换-低维流形模型的地震去噪与插值方法 | |
CN116610937B (zh) | 在隐式空间上进行低频信息延拓的方法、装置及电子设备 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |