CN117434604A - 地震数据处理方法、装置、存储介质及电子设备 - Google Patents
地震数据处理方法、装置、存储介质及电子设备 Download PDFInfo
- Publication number
- CN117434604A CN117434604A CN202310183922.XA CN202310183922A CN117434604A CN 117434604 A CN117434604 A CN 117434604A CN 202310183922 A CN202310183922 A CN 202310183922A CN 117434604 A CN117434604 A CN 117434604A
- Authority
- CN
- China
- Prior art keywords
- seismic data
- seismic
- initial
- bandwidth
- determining
- 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.)
- Granted
Links
- 238000003672 processing method Methods 0.000 title claims abstract description 18
- 238000000034 method Methods 0.000 claims abstract description 56
- 238000012545 processing Methods 0.000 claims abstract description 38
- 230000006870 function Effects 0.000 claims description 28
- 238000004364 calculation method Methods 0.000 claims description 19
- 238000002922 simulated annealing Methods 0.000 claims description 11
- 238000004590 computer program Methods 0.000 claims description 8
- 238000005516 engineering process Methods 0.000 abstract description 7
- 238000010586 diagram Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 239000002131 composite material Substances 0.000 description 4
- 238000005070 sampling Methods 0.000 description 4
- 238000005553 drilling Methods 0.000 description 3
- 230000003068 static effect Effects 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 239000010410 layer Substances 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000000985 reflectance spectrum Methods 0.000 description 2
- KLDZYURQCUYZBL-UHFFFAOYSA-N 2-[3-[(2-hydroxyphenyl)methylideneamino]propyliminomethyl]phenol Chemical compound OC1=CC=CC=C1C=NCCCN=CC1=CC=CC=C1O KLDZYURQCUYZBL-UHFFFAOYSA-N 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 201000001098 delayed sleep phase syndrome Diseases 0.000 description 1
- 208000033921 delayed sleep phase type circadian rhythm sleep disease Diseases 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000011229 interlayer Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000013517 stratification Methods 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
-
- 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
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
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
技术领域
本发明涉及地震资料数据处理技术领域，具体涉及一种地震数据处理方法、装置、存储介质及电子设备。
背景技术
地震数据的高分辨率处理，主要是通过反褶积来实现。按照数学算法划分，反褶积方法可分为两类：第一类为线性反褶积，如最小平方反褶积，其特点是不能产生高频，只在有效信号占绝大部分的频带中，通过抬高其中的高频达到提高分辨率的目的；第二类是非线性反褶积，如稀疏约束反褶积，其特点是通过加入各种稀疏约束项，使反褶积变成非线性问题，由于所加的各种约束项具备不同程度创造高频的能力，从而克服了线性反褶积带限的缺陷，能大幅度的提高地震资料的分辨率。
目前常用的一类非线性反褶积方法是稀疏约束反演，这类方法在地层反射系数是稀疏的前提下，通过在反演中加入各种稀疏约束项，达到提高地震资料分辨率的目的，对于相对简单的地质目标而言，具有较强的适用性；但是，在实际应用中，地层反射系数不是严格稀疏的，由于稀疏约束反演出的高频信息不受地震资料本身所控制，而且各种稀疏约束是定性约束，因此，当面对复杂的地质目标时，当前基于稀疏约束反演的地震数据处理技术容易造成反演结果与地震资料的差异较大而无法判断真伪的问题。
故，当前地震数据处理技术存在稀疏约束反演方式导致反演结果不准确的技术问题。
发明内容
为缓解当前地震数据处理技术存在的反演结果不准确的技术问题，本发明实施例提供一种地震数据处理方法、装置、存储介质及电子设备。
第一方面，本发明实施例提供一种地震数据处理方法，所述方法包括：
获取初始地震数据对应的地震子波；
确定所述初始地震数据的初始带宽；
根据预设的噪声约束条件以及所述初始带宽，确定拓宽处理后的地震反射系数；
根据所述地震反射系数与所述地震子波，得到高分辨率地震数据。
在一些实施方式中，所述根据预设的噪声约束条件以及所述初始带宽，确定拓宽处理后的地震反射系数的步骤，包括：
确定反射系数计算方式及地震信号计算方式；
根据所述噪声约束条件，确定高频噪声定量约束函数；
根据所述高频噪声定量约束函数，对所述初始带宽进行扩宽处理得到扩宽带宽；
获取所述扩宽带宽对应的所述地震反射系数。
在一些实施方式中，所述根据所述高频噪声定量约束函数，对所述初始带宽进行扩宽处理得到扩宽带宽的步骤，包括：
采用模拟退火方式求解所述高频噪声定量约束函数，得到中间地震数据；
采用迭代方式，基于所述中间地震数据确定所述扩宽带宽。
在一些实施方式中，所述采用模拟退火方式求解所述高频噪声定量约束函数，得到中间地震数据的步骤，包括：
根据所述初始带宽得到所述高频噪声定量约束函数的初始值；
基于所述初始带宽得到中间解对应的中间值；
基于所述初始值、所述中间值、迭代层数以及约束条件，得到所述中间地震数据。
在一些实施方式中，所述采用迭代方式，基于所述中间地震数据确定所述扩宽带宽的步骤，包括：
根据所述初始地震数据与所述中间地震数据，得到残差数据体；
根据所述残差数据体中的有效数据体，确定所述扩宽带宽。
在一些实施方式中，在所述根据所述残差数据体中的有效数据体，确定所述扩宽带宽的步骤之前，还包括：
根据所述有效数据体和所述中间数据体，得到候选地震数据；
在所述候选地震数据满足分辨率要求时，输出所述地震反射系数；
在所述候选地震数据不满足分辨率要求时，根据所述残差数据体中的有效数据体，确定所述扩宽带宽。
在一些实施方式中，所述确定所述初始地震数据的初始带宽的步骤，包括：
对所述初始地震数据进行分频扫描，得到所述初始地震数据的信噪比；
根据信噪比大于预设值的频率带宽，确定所述初始带宽。
第二方面，本发明实施例提供一种地震数据处理装置，所述装置包括：
获取模块，用于获取初始地震数据对应的地震子波；
第一确定模块，用于确定所述初始地震数据的初始带宽；
第二确定模块，用于根据预设的噪声约束条件以及所述初始带宽，确定拓宽处理后的地震反射系数；
实现模块，用于根据所述地震反射系数与所述地震子波，得到高分辨率地震数据。
第三方面，本发明实施例提供一种存储介质，所述存储介质上存储有计算机程序，所述计算机程序被一个或多个处理器执行时，实现如第一方面所述的方法。
第四方面，本发明实施例提供一种电子设备，包括存储器和处理器，所述存储器上存储有计算机程序，所述计算机程序被所述处理器执行时实现如第一方面所述的方法。
与现有技术相比，本发明的一个或多个实施例至少能带来如下有益效果：
本发明实施例提供了一种地震数据处理方法、装置、存储介质及电子设备；所述方法包括获取初始地震数据对应的地震子波，确定所述初始地震数据的初始带宽，根据预设的噪声约束条件以及所述初始带宽，确定拓宽处理后的地震反射系数，根据所述地震反射系数与所述地震子波，得到高分辨率地震数据。在本方案提供的方法中，基于初始地震数据及预设的噪声约束条件，采用多次迭代方法，逐步拓宽地震数据的频带宽度，摆脱了当前方法需要在“有限带宽”内提高地震数据分辨率的束缚，实现了在通放带内提高地震数据分辨率的目的，缓解了当前基于稀疏约束反演的地震数据处理技术容易造成反演结果与地震资料的差异较大而无法判断真伪的技术问题；同时，提高分辨率处理后的地震数据的信噪比高、信号保真度高，能够很好的完成地质任务。
附图说明
为了更清楚地说明本发明实施例的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本发明的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。
图1为本申请实施例提供的地震数据处理方法的第一种流程示意图；
图2为本申请实施例提供的地震数据处理方法的第二种流程示意图；
图3为本申请实施例提供的地震反射系数分解示意图；
图4为本申请实施例提供的模拟退火法计算示意图；
图5为本申请实施例提供的迭代法求最终反射系数流程图；
图6为本申请实施例提供的CR工区初次迭代与最终迭代(第三次迭代)地震剖面对比图；
图7为本申请实施例提供的MM地区S2井测井计算反射系数与最终反演的地震反射系数对比图；
图8为本申请实施例提供的ZZ地区页岩反演前后对比剖面图；
图9为本申请实施例提供的TD地区A1-A2-A3连井剖面反演前后储层分布对比图；
图10为本申请实施例提供的KD地区KD893井-KD894井连井剖面反演前后对比图；
图11为本申请实施例提供的地震数据处理装置的结构示意图。
在附图中，相同的部件使用相同的附图标记，附图并未按照实际的比例绘制。
具体实施方式
以下将结合附图及实施例来详细说明本申请的实施方式，借此对本申请如何应用技术手段来解决技术问题，并达到相应技术效果的实现过程能充分理解并据以实施。本申请实施例以及实施例中的各个特征，在不相冲突前提下可以相互结合，所形成的技术方案均在本申请的保护范围之内。
同时，在以下说明中，出于解释的目的而阐述了许多具体细节，以提供对本发明实施例的彻底理解。然而，对本领域的技术人员来说显而易见的是，本发明可以不用这里的具体细节或者所描述的特定方式来实施。
实施例一
图1为本申请实施例提供的地震数据处理方法的第一种流程示意图，请参阅图1，本实施例提供的地震数据处理方法包括：
步骤S110：获取初始地震数据对应的地震子波。
在本申请中，初始地震数据包括叠后地震数据；从叠后地震数据中精细提取时变地震子波W(t)，即完成获取初始地震数据对应的地震子波的步骤。
步骤S120：确定所述初始地震数据的初始带宽。
在一些实施方式中，本步骤包括：对所述初始地震数据进行分频扫描，得到所述初始地震数据的信噪比；根据信噪比大于预设值的频率带宽，确定所述初始带宽。
在本申请中，初始带宽为初始地震数据的信噪比大于预设值的频率带宽，也即有效频宽；此处，预设值可以为1。
具体的，对初始地震数据进行分频扫描，计算信噪比的公式为：
式中：SNR为信噪比，M为分析时窗的采样点数，N为地震道数，hij是地震记录振幅；根据信噪比大于1的频率带宽来确定“有效频宽”(p1 p2)。
步骤S130：根据预设的噪声约束条件以及所述初始带宽，确定拓宽处理后的地震反射系数。
在一些实施方式中，本步骤包括：确定反射系数计算方式及地震信号计算方式；根据所述噪声约束条件，确定高频噪声定量约束函数；根据所述高频噪声定量约束函数，对所述初始带宽进行扩宽处理得到扩宽带宽；获取所述扩宽带宽对应的所述地震反射系数。
在本申请中，反射系数计算方式可以选取奇偶反射系数的计算公式：
如图3所示，设某一地震反射道的ti时间位置的反射系数为R(ti)，则：
R(ti)＝αi+βi；
式中αi、βi分别代表ti时间位置的奇偶反射系数，根据傅里叶变换可以得到：
在本申请中，地震信号计算公式为：
根据地震数据处理原理，在频率域内，地震信号可以表示为：
H(f)＝W(f)R(f)+N(f) (3)
式中：H(f)为地震记录，W(f)为地震子波，R(f)为反射系数，N(f)为噪声，f表示地震道的频率。根据公式(2)可以得到：
式中：αi为奇反射系数，βi为偶反射系数，ti为产生反射系数的道所对应的时间。
假设地震子波W(f)已知，频率用fc表示，则反射系数谱方程表示为：
式中：分别表示反射系数和噪声的反褶积结果，fc＝(f1，f2，f3……fc)，即C表示频率个数，则其向量表达式为：
A＝Qg+N (6)
式中：A＝[H(f1)/W(f1)，H(f2)/W(f2)，H(f3)/W(f3)…，H(fc)/W(fc)]T
在本申请中，高频噪声定量约束函数具体为：
设高频噪声定量约束函数由F(p1-2有效带宽，t，α，β)表达，目标函数采用最小二乘法使噪声达到最小值。其表达式如下：
式中，k表示时窗内的采样点，η表示地震道数，p为初始有效带宽(p1、p2)。
在一些实施方式中，根据所述高频噪声定量约束函数，对所述初始带宽进行扩宽处理得到扩宽带宽的步骤包括：采用模拟退火方式求解所述高频噪声定量约束函数，得到中间地震数据；采用迭代方式，基于所述中间地震数据确定所述扩宽带宽。
在一些实施方式中，采用模拟退火方式求解所述高频噪声定量约束函数，得到中间地震数据的步骤包括：根据所述初始带宽得到所述高频噪声定量约束函数的初始值；基于所述初始带宽得到中间解对应的中间值；基于所述初始值、所述中间值、迭代层数以及约束条件，得到所述中间地震数据。
在本申请中，中间地震数据是指合成数据体(Si)。
具体的，如图4所示，在采用模拟退火方式求解高频噪声定量约束函数方程的计算中，给出初始有效带宽(p1、p2)，初始ti及对应的初始反射系数αi、βi；按照模拟退火法进行逐一计算，最终得到反射系数体(计算流程请参阅图4)，反射系数体与子波(地震子波)褶积得到合成数据体(Si)。
在一些实施方式中，采用迭代方式，基于所述中间地震数据确定所述扩宽带宽的步骤包括：根据所述初始地震数据与所述中间地震数据，得到残差数据体；根据所述残差数据体中的有效数据体，确定所述扩宽带宽。
在本申请中，有效数据体是指对残差数据体进行去噪后的数据体。
具体的，利用原始数据体(即初始地震数据)与合成数据体(Si)(即中间地震数据)相减得到最小残差Sfe Error1数据体。用下式表示：
Sfe Error1＝Sfe Residual1+N1 fe (8)
式中：Sfe Error1表示最小残差Sfe Error1数据体，Sfe Residual1表示最小残差中剩余的有效信号，N1 fe表示最小残差中的噪声信号。
对Sfe Error1数据体进行去噪，对去噪后的数据体(即有效数据体)进行分频扫描，将高频扩展到p3(即扩宽带宽)。
在一些实施方式中，在所述根据所述残差数据体中的有效数据体，确定所述扩宽带宽的步骤之前，还包括：根据所述有效数据体和所述中间数据体，得到候选地震数据；在所述候选地震数据满足分辨率要求时，输出所述地震反射系数；在所述候选地震数据不满足分辨率要求时，根据所述残差数据体中的有效数据体，确定所述扩宽带宽。
在本申请中，候选地震数据是指去噪后的地震数据Si+1，该候选地震数据是利用有效数据体与中间数据体相加得到的。
具体的，如图5所示，在所述候选地震数据满足分辨率要求时，输出所述地震反射系数；在所述候选地震数据不满足分辨率要求时，根据所述残差数据体中的剩余有效数据体，保持低频端p1不变，利用p3作为初始高频pi新进行新一轮迭代(频率从最初的p1-p2拓宽致p1-pm)。
即频带宽度变化为(p1-p2；p1-p3；p1-p4……p1-pm)；直到残差结果最小时，确定最终t所对应α、β取值，用于计算最终地震反射系数(具体流程请参阅图5)。
步骤S140：根据所述地震反射系数与所述地震子波，得到高分辨率地震数据。
如图6所示，图6是CR工区初次迭代(上)与最终迭代(第三次迭代S3)后XX60-XX49-XX69测线地震剖面对比图，从上、下地震剖面的对比分析可以看出，地震资料的分辨率有了极大的提高，其主频在45Hz左右，频带宽度从5-100Hz拓宽到达到5Hz-140Hz，能很好的完成地质任务。
图2为本申请实施例提供的地震数据处理方法的第二种流程示意图，请参阅图2，本实施例提供的地震数据处理方法包括：
S1：提取地震子波W(t)。
在本申请中，初始地震数据包括叠后地震数据；从叠后地震数据中精细提取时变地震子波W(t)。
S2：确定地震资料(对应上文中的初始地震数据)有效带宽(对应上文中的初始带宽)。
对叠后处理的地震资料进行分频扫描，计算信噪比：
式中：SNR为信噪比，M为分析时窗的采样点数，N为地震道数，hij是地震记录振幅；根据信噪比大于1的频率带宽来确定“有效频宽”(p1 p2)。
S3：设计高频噪声定量约束最大概率准则及约束函数，求取最终拓频后的地震反射系数；
设计最大概率准则中包括地震资料的有效频带宽度、地震反射时间及其对应的反射系数的门槛值等相关控制参数；
S31：设计计算公式：
其中，奇偶反射系数的计算公式如下：
设某一地震反射道的ti时间位置的反射系数为R(ti)，则：
R(ti)＝αi+βi(请参见图3)；
式中αi、βi分别代表ti时间位置的奇偶反射系数，根据傅里叶变换：
其中，地震信号的表达公式如下：
根据地震数据处理原理，在频率域内，地震信号可以表示为：
H(f)＝W(f)R(f)+N(f)
式中：H(f)为地震记录，W(f)为地震子波，R(f)为反射系数，N(f)为噪声，f表示地震道的频率。根据公式可以得到：
式中：式中：αi为奇反射系数，βi为偶反射系数，ti为产生反射系数的道所对应的时间。
假设地震子波W(f)已知，频率用fc表示，则反射系数谱方程表示为：
式中：分别表示反射系数和噪声的反褶积结果，fc＝(f1，f2，f3，……fc)，即C表示频率个数，则其向量表达式为：
A＝Qg+N
式中：A＝[H(f1)/W(f1)，H(f2)/W(f2)，H(f3)/W(f3)…，H(fc)/W(fc)]T
S32：设计高频噪声定量约束函数。
设高频噪声定量约束函数由F(p1-2有效带宽，t，α，β)表达，目标函数采用最小二乘法使噪声达到最小值。其表达式如下：
式中，k表示时窗内的采样点，η表示地震道数，p为初始有效带宽(p1、p2)；
S33：采用模拟退火法求解高频噪声定量约束函数方程。
计算中，给出初始有效带宽(p1、p2)，初始ti及对应的初始反射系数αi、βi；按照模拟退火法进行逐一计算，最终得到反射系数体(计算流程请参见图4)，反射系数体与子波褶积得到合成数据体(Si)(对应上文中的中间地震数据)。
S34：逐步迭代，求取最终拓频后的反射系数(请参见图5)；
①利用原始数据体与合成数据体(Si体)相减得到最小残差Sfe Error1数据体。用下式表示：
Sfe Error1＝Sfe Residual1+N1 fe
式中：Sfe Error1表示最小残差Sfe Error1数据体，Sfe Residual1表示最小残差中剩余的有效信号，N1 fe表示最小残差中的噪声信号。
②对Sfe Error1数据体进行去噪，对去噪后的数据体(对应上文中的有效数据体)进行分频扫描，将高频扩展到p3；利用去噪后的结果与合成数据Si相加得到去噪后的地震数据Si+1(对应上文中的候选地震数据)。
③低频端p1不变，利用p3作为初始高频pi新进行新一轮迭代(频率从最初的p1-p2拓宽致p1-pm)。
即频带宽度变化为(p1-p2；p1-p3；p1-p4……p1-pm)；直到残差结果最小时，确定最终t所对应α、β取值，用于计算最终反射系数。
S4：用高频噪声定量约束函数计算得到的最终地震反射系数与地震子波褶积形成高分辨率数据体。
请参阅图6，图6是CR工区初次迭代(上)与最终迭代(第三次迭代S3)后XX60-XX49-XX69测线地震剖面对比图，如图6所示，从上、下地震剖面的对比分析可以看出，地震资料的分辨率有了极大的提高，其主频在45Hz左右，频带宽度从5-100Hz拓宽到达到5Hz-140Hz，能很好的完成地质任务。
综上，本实施例基于初始地震数据及预设的噪声约束条件，采用多次迭代方法，逐步拓宽地震数据的频带宽度，摆脱了当前方法需要在“有限带宽”内提高地震数据分辨率的束缚，实现了在通放带内提高地震数据分辨率的目的，缓解了当前基于稀疏约束反演的地震数据处理技术容易造成反演结果与地震资料的差异较大而无法判断真伪的技术问题；同时，提高分辨率处理后的地震数据的信噪比高、信号保真度高，能够很好的完成地质任务。
实施例二
图11为本申请实施例提供的地震数据处理装置的结构示意图，请参阅图11，本实施例提供的地震数据处理装置包括：
获取模块1110，用于获取初始地震数据对应的地震子波；
第一确定模块1120，用于确定所述初始地震数据的初始带宽；
第二确定模块1130，用于根据预设的噪声约束条件以及所述初始带宽，确定拓宽处理后的地震反射系数；
实现模块1140，用于根据所述地震反射系数与所述地震子波，得到高分辨率地震数据。
基于上述各模块执行方法的具体实施例已在实施例一中详述，此处不再赘述。
现基于前述的记载，结合具体场景对本申请进行详细的说明。
利用本申请的地震数据处理方法对多个工区叠后三维数据体开展高分辨率地震数据处理，结果显示，利用测井资料计算的地层反射系数与经过高分辨率数据处理后抽取的井旁地震道地震反射系数对应关系良好。几个地震工区原始地震资料与叠后高分辨率地震数据处理结果相比，后者不仅地震频带有效展宽，分辨率得到不同程度的提高，而且，地质现象更加清楚，解决地质问题的能力得到有效提升。
请参阅图7，图7为MM地区S2井测井计算反射系数与反演成果对比图。如图7所示，曲线a代表利用测井资料计算的地层反射系数，曲线b为本发明反演的地层反射系数。从图中可以看出，利用S2井实际测井资料计算的地震反射系数与利用发明抽取的井旁地震道反射系数的相对关系良好，表明反演结果真是可靠。
请参阅图8，图8为ZZ地区页岩反演前后对比剖面图。如图8所示，从图中可以看出，原始地震剖面上页岩整体为一套强反射特征(上)，利用本发明处理的高分辨率地震资料，可以将原来的一套强反射分解为上、下2套强反射(下)。
请参阅图9，图9为TD地区A1-A2-A3连井剖面反演前后储层分布对比图。如图9所示，从图中可以看出，原始地震剖面上三口井均钻遇1套连续发布的厚储层(虚线c所示)，而如曲线d(3个椭圆形曲线)所示，实际钻井资料揭示为2套薄储层被1套泥岩夹层隔开；利用本发明处理的高分辨率地震资料显示，2套薄储层不仅具有横向不连续性，而且分布不均，与实际钻井资料吻合良好。
请参阅图10，图10为KD地区KD893-KD894井2口开发井的连井剖面反演前后对比图，钻井后，地震地质分层有很大的争议。利用本发明处理的高分辨率地震资料显示，2口井的地震地质层位清楚，如画圈区域所示，由于地震分辨率提高，地质现象清楚。
实施例三
本实施例还提供一种计算机可读存储介质，如闪存、硬盘、多媒体卡、卡型存储器(例如，SD或DX存储器等)、随机访问存储器(RAM)、静态随机访问存储器(SRAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、可编程只读存储器(PROM)、磁性存储器、磁盘、光盘、服务器、App应用商城等等，该存储介质上存储有计算机程序，所述计算机程序被处理器执行时可以实现如下方法步骤：
获取初始地震数据对应的地震子波；
确定所述初始地震数据的初始带宽；
根据预设的噪声约束条件以及所述初始带宽，确定拓宽处理后的地震反射系数；
根据所述地震反射系数与所述地震子波，得到高分辨率地震数据。
上述方法步骤的具体实施例过程可参见实施例一，本实施例在此不再重复赘述。
实施例四
本实施例提供一种电子设备，该电子设备可以是手机、电脑或平板电脑等作为本申请中的分享服务器，包括存储器和处理器，所述存储器上存储有计算器程序，该计算机程序被处理器执行时实现如实施例一中所述的地震数据处理方法。可以理解，电子设备还可以包括，输入/输出(I/O)接口，以及通信组件。
其中，处理器用于执行如实施例一中的地震数据处理方法中的全部或部分步骤。存储器用于存储各种类型的数据，这些数据例如可以包括终端设备中的任何应用程序或方法的指令，以及应用程序相关的数据。
所述处理器可以是专用集成电路(Application Specific Integrated Circuit，简称ASIC)、数字信号处理器(Digital Signal Processor，简称DSP)、数字信号处理设备(Digital Signal Processing Device，简称DSPD)、可编程逻辑器件(Programmable LogicDevice，简称PLD)、现场可编程门阵列(Field Programmable Gate Array，简称FPGA)、控制器、微控制器、微处理器或其他电子元件实现，用于执行上述实施例一中的地震数据处理方法。
所述存储器可以由任何类型的易失性或非易失性存储设备或者它们的组合实现，例如静态随机存取存储器(Static Random Access Memory，简称SRAM)，电可擦除可编程只读存储器(Electrically Erasable Programmable Read-Only Memory，简称EEPROM)，可擦除可编程只读存储器(Erasable Programmable Read-Only Memory，简称EPROM)，可编程只读存储器(Programmable Read-Only Memory，简称PROM)，只读存储器(Read-Only Memory，简称ROM)，磁存储器，快闪存储器，磁盘或光盘。
综上，本申请提供一种地震数据处理方法、装置、存储介质及电子设备；所述方法包括获取初始地震数据对应的地震子波，确定所述初始地震数据的初始带宽，根据预设的噪声约束条件以及所述初始带宽，确定拓宽处理后的地震反射系数，根据所述地震反射系数与所述地震子波，得到高分辨率地震数据。在本方案提供的方法中，基于初始地震数据及预设的噪声约束条件，采用多次迭代方法，逐步拓宽地震数据的频带宽度，摆脱了当前方法需要在“有限带宽”内提高地震数据分辨率的束缚，实现了在通放带内提高地震数据分辨率的目的，缓解了当前基于稀疏约束反演的地震数据处理技术容易造成反演结果与地震资料的差异较大而无法判断真伪的技术问题；同时，提高分辨率处理后的地震数据的信噪比高、信号保真度高，能够很好的完成地质任务。
需要说明的是，在本文中，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
虽然本申请所揭露的实施方式如上，但所述的内容只是为了便于理解本申请而采用的实施方式，并非用以限定本申请。任何本申请所属技术领域内的技术人员，在不脱离本申请所揭露的精神和范围的前提下，可以在实施的形式上及细节上作任何的修改与变化，但本申请的专利保护范围，仍须以所附的权利要求书所界定的范围为准。
Claims (10)
1.一种地震数据处理方法，其特征在于，包括：
获取初始地震数据对应的地震子波；
确定所述初始地震数据的初始带宽；
根据预设的噪声约束条件以及所述初始带宽，确定拓宽处理后的地震反射系数；
根据所述地震反射系数与所述地震子波，得到高分辨率地震数据。
2.根据权利要求1所述的地震数据处理方法，其特征在于，所述根据预设的噪声约束条件以及所述初始带宽，确定拓宽处理后的地震反射系数的步骤，包括：
确定反射系数计算方式及地震信号计算方式；
根据所述噪声约束条件，确定高频噪声定量约束函数；
根据所述高频噪声定量约束函数，对所述初始带宽进行扩宽处理得到扩宽带宽；
获取所述扩宽带宽对应的所述地震反射系数。
3.根据权利要求2所述的地震数据处理方法，其特征在于，所述根据所述高频噪声定量约束函数，对所述初始带宽进行扩宽处理得到扩宽带宽的步骤，包括：
采用模拟退火方式求解所述高频噪声定量约束函数，得到中间地震数据；
采用迭代方式，基于所述中间地震数据确定所述扩宽带宽。
4.根据权利要求3所述的地震数据处理方法，其特征在于，所述采用模拟退火方式求解所述高频噪声定量约束函数，得到中间地震数据的步骤，包括：
根据所述初始带宽得到所述高频噪声定量约束函数的初始值；
基于所述初始带宽得到中间解对应的中间值；
基于所述初始值、所述中间值、迭代层数以及约束条件，得到所述中间地震数据。
5.根据权利要求4所述的地震数据处理方法，其特征在于，所述采用迭代方式，基于所述中间地震数据确定所述扩宽带宽的步骤，包括：
根据所述初始地震数据与所述中间地震数据，得到残差数据体；
根据所述残差数据体中的有效数据体，确定所述扩宽带宽。
6.根据权利要求5所述的地震数据处理方法，其特征在于，在所述根据所述残差数据体中的有效数据体，确定所述扩宽带宽的步骤之前，还包括：
根据所述有效数据体和所述中间数据体，得到候选地震数据；
在所述候选地震数据满足分辨率要求时，输出所述地震反射系数；
在所述候选地震数据不满足分辨率要求时，根据所述残差数据体中的有效数据体，确定所述扩宽带宽。
7.根据权利要求1至6任一项所述的地震数据处理方法，其特征在于，所述确定所述初始地震数据的初始带宽的步骤，包括：
对所述初始地震数据进行分频扫描，得到所述初始地震数据的信噪比；
根据信噪比大于预设值的频率带宽，确定所述初始带宽。
8.一种地震数据处理装置，其特征在于，包括：
获取模块，用于获取初始地震数据对应的地震子波；
第一确定模块，用于确定所述初始地震数据的初始带宽；
第二确定模块，用于根据预设的噪声约束条件以及所述初始带宽，确定拓宽处理后的地震反射系数；
实现模块，用于根据所述地震反射系数与所述地震子波，得到高分辨率地震数据。
9.一种存储介质，其特征在于，所述存储介质上存储有计算机程序，所述计算机程序被一个或多个处理器执行时，实现如权利要求1至7中任一项所述的方法。
10.一种电子设备，其特征在于，包括存储器和处理器，所述存储器上存储有计算机程序，所述计算机程序被所述处理器执行时实现如权利要求1至7中任一项所述的方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310183922.XA CN117434604B (zh) | 2023-02-20 | 2023-02-20 | 地震数据处理方法、装置、存储介质及电子设备 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310183922.XA CN117434604B (zh) | 2023-02-20 | 2023-02-20 | 地震数据处理方法、装置、存储介质及电子设备 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117434604A true CN117434604A (zh) | 2024-01-23 |
| CN117434604B CN117434604B (zh) | 2024-08-13 |
Family
ID=89552248
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310183922.XA Active CN117434604B (zh) | 2023-02-20 | 2023-02-20 | 地震数据处理方法、装置、存储介质及电子设备 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117434604B (zh) |
Citations (20)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6393366B1 (en) * | 2000-07-07 | 2002-05-21 | Saudi Arabian Oil Company | Deconvolution of seismic data based on fractionally integrated noise |
| CN101609161A (zh) * | 2009-07-17 | 2009-12-23 | 中国石化集团胜利石油管理局 | 基于地震层序体理论多尺度资料联合频带拓展方法 |
| US20140102694A1 (en) * | 2012-10-12 | 2014-04-17 | Rock Solid Images Inc | Geophysical surveying |
| CN104122588A (zh) * | 2014-07-30 | 2014-10-29 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 基于谱分解提高叠后地震资料分辨率的方法 |
| CN104181589A (zh) * | 2014-08-20 | 2014-12-03 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种非线性反褶积方法 |
| CN104237945A (zh) * | 2014-09-18 | 2014-12-24 | 中国石油集团东方地球物理勘探有限责任公司 | 一种地震资料自适应高分辨处理方法 |
| CN107436450A (zh) * | 2017-07-26 | 2017-12-05 | 西安交通大学 | 一种基于连续小波变换的地震信号带宽拓展方法 |
| CN109061764A (zh) * | 2018-09-07 | 2018-12-21 | 中国石油化工股份有限公司 | 一种分频融合波阻抗反演方法 |
| CN110174702A (zh) * | 2018-09-30 | 2019-08-27 | 中海油海南能源有限公司 | 一种海上地震数据低频弱信号恢复的方法和系统 |
| CN111025395A (zh) * | 2020-01-06 | 2020-04-17 | 中国石油化工股份有限公司 | 基于压缩感知的高斯频率域提高薄互层分辨率的方法 |
| CN111766624A (zh) * | 2020-06-22 | 2020-10-13 | 南京加宝囤信息科技有限公司 | 一种地震数据拓频处理方法、装置、储存介质及电子设备 |
| CN112213773A (zh) * | 2019-07-12 | 2021-01-12 | 中国石油化工股份有限公司 | 一种地震分辨率提高方法及电子设备 |
| CN112363222A (zh) * | 2020-10-28 | 2021-02-12 | 中国石油天然气集团有限公司 | 叠后自适应宽带约束波阻抗反演方法及装置 |
| CN112444867A (zh) * | 2019-08-30 | 2021-03-05 | 中国石油化工股份有限公司 | 一种基于压缩感知的地震资料频带拓展处理方法 |
| CN112558161A (zh) * | 2020-12-24 | 2021-03-26 | 中国石油大学(华东) | 基于压缩感知的井约束地震低频恢复方法 |
| CN112882099A (zh) * | 2021-02-25 | 2021-06-01 | 中海石油深海开发有限公司 | 一种地震频带拓宽方法、装置、介质及电子设备 |
| CN112925013A (zh) * | 2021-01-28 | 2021-06-08 | 中国石油化工股份有限公司 | 基于全频带延拓保真的地震数据高分辨率处理方法 |
| CN114428291A (zh) * | 2020-09-22 | 2022-05-03 | 中国石油化工股份有限公司 | 多道的反射系数脉冲确定方法、存储介质及电子设备 |
| CN114428280A (zh) * | 2021-12-20 | 2022-05-03 | 中国石油化工股份有限公司 | 一种地震数据低频信息补偿方法及其应用 |
| CN114966856A (zh) * | 2022-08-02 | 2022-08-30 | 中国科学院地质与地球物理研究所 | 基于多频带地震资料的碳封存场址优选方法、系统和设备 |
-
2023
- 2023-02-20 CN CN202310183922.XA patent/CN117434604B/zh active Active
Patent Citations (20)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6393366B1 (en) * | 2000-07-07 | 2002-05-21 | Saudi Arabian Oil Company | Deconvolution of seismic data based on fractionally integrated noise |
| CN101609161A (zh) * | 2009-07-17 | 2009-12-23 | 中国石化集团胜利石油管理局 | 基于地震层序体理论多尺度资料联合频带拓展方法 |
| US20140102694A1 (en) * | 2012-10-12 | 2014-04-17 | Rock Solid Images Inc | Geophysical surveying |
| CN104122588A (zh) * | 2014-07-30 | 2014-10-29 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 基于谱分解提高叠后地震资料分辨率的方法 |
| CN104181589A (zh) * | 2014-08-20 | 2014-12-03 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种非线性反褶积方法 |
| CN104237945A (zh) * | 2014-09-18 | 2014-12-24 | 中国石油集团东方地球物理勘探有限责任公司 | 一种地震资料自适应高分辨处理方法 |
| CN107436450A (zh) * | 2017-07-26 | 2017-12-05 | 西安交通大学 | 一种基于连续小波变换的地震信号带宽拓展方法 |
| CN109061764A (zh) * | 2018-09-07 | 2018-12-21 | 中国石油化工股份有限公司 | 一种分频融合波阻抗反演方法 |
| CN110174702A (zh) * | 2018-09-30 | 2019-08-27 | 中海油海南能源有限公司 | 一种海上地震数据低频弱信号恢复的方法和系统 |
| CN112213773A (zh) * | 2019-07-12 | 2021-01-12 | 中国石油化工股份有限公司 | 一种地震分辨率提高方法及电子设备 |
| CN112444867A (zh) * | 2019-08-30 | 2021-03-05 | 中国石油化工股份有限公司 | 一种基于压缩感知的地震资料频带拓展处理方法 |
| CN111025395A (zh) * | 2020-01-06 | 2020-04-17 | 中国石油化工股份有限公司 | 基于压缩感知的高斯频率域提高薄互层分辨率的方法 |
| CN111766624A (zh) * | 2020-06-22 | 2020-10-13 | 南京加宝囤信息科技有限公司 | 一种地震数据拓频处理方法、装置、储存介质及电子设备 |
| CN114428291A (zh) * | 2020-09-22 | 2022-05-03 | 中国石油化工股份有限公司 | 多道的反射系数脉冲确定方法、存储介质及电子设备 |
| CN112363222A (zh) * | 2020-10-28 | 2021-02-12 | 中国石油天然气集团有限公司 | 叠后自适应宽带约束波阻抗反演方法及装置 |
| CN112558161A (zh) * | 2020-12-24 | 2021-03-26 | 中国石油大学(华东) | 基于压缩感知的井约束地震低频恢复方法 |
| CN112925013A (zh) * | 2021-01-28 | 2021-06-08 | 中国石油化工股份有限公司 | 基于全频带延拓保真的地震数据高分辨率处理方法 |
| CN112882099A (zh) * | 2021-02-25 | 2021-06-01 | 中海石油深海开发有限公司 | 一种地震频带拓宽方法、装置、介质及电子设备 |
| CN114428280A (zh) * | 2021-12-20 | 2022-05-03 | 中国石油化工股份有限公司 | 一种地震数据低频信息补偿方法及其应用 |
| CN114966856A (zh) * | 2022-08-02 | 2022-08-30 | 中国科学院地质与地球物理研究所 | 基于多频带地震资料的碳封存场址优选方法、系统和设备 |
Non-Patent Citations (3)
| Title |
|---|
| 王笋等: "大容量气枪源广角地震数据的有效频带研究", 2018年中国地球科学联合学术年会论文集（二十二）——专题45：海洋地球物理、专题46：电磁地球物理学研究应用及其新进展, 21 October 2018 (2018-10-21) * |
| 王笋等: "大容量气枪震源广角地震数据的高频成分分析", 地球物理学进展, vol. 34, no. 1, 28 February 2019 (2019-02-28) * |
| 黄林军等: "基于对时变子波进行分频段处理的反褶积方法在薄层预测中的应用", 物探化探计算技术, vol. 41, no. 5, 30 September 2019 (2019-09-30) * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117434604B (zh) | 2024-08-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Gaci | The use of wavelet-based denoising techniques to enhance the first-arrival picking on seismic traces | |
| US20140019055A1 (en) | System and method for estimating and attenuating noise in seismic data | |
| CN110208862B (zh) | 一种基于混合高阶分数阶ATpV稀疏正则化的地震反演方法 | |
| CN113568058B (zh) | 一种基于多分辨率奇异值分解的大地电磁信噪分离方法及系统 | |
| Al-Dossary et al. | Lineament-preserving filtering | |
| US11662492B2 (en) | Seismic random noise attenuation | |
| Bagheri et al. | Denoising and improving the quality of seismic data using combination of DBM filter and FX deconvolution | |
| Nose-Filho et al. | Algorithms for sparse multichannel blind deconvolution | |
| CN117434604B (zh) | 地震数据处理方法、装置、存储介质及电子设备 | |
| CN113655522A (zh) | 频率域地震弱信号的增强方法 | |
| Qian et al. | Unsupervised intense VSP coupling noise suppression with iterative robust deep learning | |
| Jahanjooy et al. | Blind inversion of multidimensional seismic data using sequential Tikhonov and total variation regularizations | |
| Gómez et al. | Noise suppression in 2D and 3D seismic data with data-driven sifting algorithms | |
| CN116626760A (zh) | 基于自适应高阶最大熵wvd的地层不连续性检测方法及装置 | |
| Wang et al. | A Unidirectional Total Variation and Second‐Order Total Variation Model for Destriping of Remote Sensing Images | |
| CN117890980A (zh) | 一种测井数据标准化处理方法、装置及设备 | |
| CN117434592B (zh) | 地震数据处理方法、装置及电子设备 | |
| Nuha et al. | Bilateral filters with elliptical Gaussian kernels for seismic surveys denoising | |
| CN111766624B (zh) | 一种地震数据拓频处理方法、装置、储存介质及电子设备 | |
| Wu et al. | Attenuating coherent environmental noise in seismic data via the U-net method | |
| Guo et al. | Multi-trace acoustic impedance inversion with multiplicative regularization | |
| Wang et al. | An empirical mode decomposition based noise cancelation method for potential field data along with a new stopping criterion | |
| Zhang et al. | Ensemble-based data assimilation with curvelets regularization | |
| Yang et al. | Denoising Seismic Data via a Threshold Shrink Method in the Non‐Subsampled Contourlet Transform Domain | |
| Al-Shaikhi et al. | Enhanced Seismic Deconvolution by Side Lobe Suppression |
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 |