CN112666601A - 一种地震数据振幅拟合有效时窗的方法及其系统 - Google Patents
一种地震数据振幅拟合有效时窗的方法及其系统 Download PDFInfo
- Publication number
- CN112666601A CN112666601A CN201910976991.XA CN201910976991A CN112666601A CN 112666601 A CN112666601 A CN 112666601A CN 201910976991 A CN201910976991 A CN 201910976991A CN 112666601 A CN112666601 A CN 112666601A
- Authority
- CN
- China
- Prior art keywords
- data
- effective
- sample
- sampling points
- amplitude
- 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
- 238000000034 method Methods 0.000 title claims abstract description 23
- 238000005070 sampling Methods 0.000 claims abstract description 104
- 238000004364 calculation method Methods 0.000 claims abstract description 20
- 238000012216 screening Methods 0.000 claims abstract description 7
- 238000007781 pre-processing Methods 0.000 claims abstract description 6
- 238000004422 calculation algorithm Methods 0.000 abstract description 10
- 230000006870 function Effects 0.000 description 12
- 230000008707 rearrangement Effects 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了地震数据振幅拟合有效时窗的方法及其系统,该方法包括:S1、对道集数据按预设规则重新排列;S2、在排列后的道集数据中筛选出有效样点的排列;S3、利用所述有效样点的排列计算有效起始道位置和终止道位置。预处理模块,用于对道集数据按预设规则重新排列;该系统包括:第一计算模块,用于在排列后的道集数据中筛选出有效样点的排列;第二计算模块,用于利用所述有效样点的排列计算有效起始道位置和终止道位置。本发明大大提高了地震数据振幅拟合有效时窗的算法效率。
Description
技术领域
本发明属于勘探领域的地震数据的高效处理技术,具体涉及一种地震数据振幅拟合有效时窗的方法及其系统。
背景技术
围绕均值剔除拟合求AVO属性,即对每个时刻采样点各道振幅拟合Shuey三项式系数求属性的计算中,如图1所示,传统的方法步骤是:从整个道集数据的第一个样点开始,(1)判断在此样点的整个道集数据是否有效,(2)找出在此样点的整个道集数据中起始有效道位置,(3)找出在此样点的整个道集数据中终止有效道位置,以此类推,直到把所有的样点的起始有效道位置和终止有效道位置都找出来,然后使用算法围绕均值剔除误差大的点,并拟合Shuey三项式抛物线得抛物线系数。这种方法的算法复杂度O的值为N*M*2,其中N为道集中道数目,M为每道的样点数目,效率十分低下。
发明内容
本发明的特征和优点在下文的描述中部分地陈述,或者可从该描述显而易见,或者可通过实践本发明而学习。
为克服现有技术的问题,本发明提供一种地震数据振幅拟合有效时窗的方法,包括:
S1、对道集数据按预设规则重新排列;
S2、在排列后的道集数据中筛选出有效样点的排列;
S3、利用所述有效样点的排列计算有效起始道位置和终止道位置。
可选地,所述预设规则为对同一时刻采集的样点赋予同一个样点号。
可选地,所述步骤S2具体包括:
使用向量指令SIMD计算同一个样点号对应的所有样点的振幅之和,获取所述样点号对应的振幅和计算值;
去除所有振幅和计算值为零的样点号对应的样点,得到K个振幅和计算值不为零的样点号及其对应的样点,记为所述有效样点的排列。
可选地,所述步骤S2具体为:使用计算机的向量计算方法一次性计算出在道集数据中有效样点的排列。
可选地,所述步骤S3具体包括:
使用fabs(data[样点号][道号])>1.0e-10&&滤波函数(data[样点号][道号])==fabs(data[样点号][道号])的条件找到相应样点处有效起始道号;
使用fabs(data[样点号][道号])<1.0e-10&&滤波函数(data[样点号][道号])==当前样点振幅和计算值的条件找到相应样点处有效终止道号。
本发明提供一种地震数据振幅拟合有效时窗的系统,包括:
预处理模块,用于对道集数据按预设规则重新排列;
第一计算模块,用于在排列后的道集数据中筛选出有效样点的排列;
第二计算模块,用于利用所述有效样点的排列计算有效起始道位置和终止道位置。
可选地,所述预设规则为对同一时刻采集的样点赋予同一个样点号。
可选地,所述第一计算模块具体用于:
使用向量指令SIMD计算同一个样点号对应的所有样点的振幅之和,获取所述样点号对应的振幅和计算值;
去除所有振幅和计算值为零的样点号对应的样点,得到K个振幅和计算值不为零的样点号及其对应的样点,记为所述有效样点的排列。
可选地,所述第二计算模块具体用于:
使用fabs(data[样点号][道号])>1.0e-10&&滤波函数(data[样点号][道号])==fabs(data[样点号][道号])的条件找到相应样点处有效起始道号;
使用fabs(data[样点号][道号])<1.0e-10&&滤波函数(data[样点号][道号])==当前样点振幅和计算值的条件找到相应样点处有效终止道号。
本发明提供一种计算机可读存储介质,所述计算机可读存储介质存储有计算机可执行的至少一个程序,所述至少一个程序被所述计算机执行时使所述计算机执行本发明任一实施例所述的方法中的步骤。
本发明提供的地震数据振幅拟合有效时窗的方法及其系统,使用道集数据重排技术,对道集数据进行预处理,然后使用计算机的向量计算方法一次性计算出在整个道集数据中有效样点的排列,最后利用此有效样点,计算出有效起始道位置和终止道位置。从而设计出了高效的道集数据振幅拟合有效时窗计算方法。
附图说明
图1为现有技术的地震数据振幅拟合有效时窗的算法伪代码。
图2为本发明实施例提供的地震数据振幅拟合有效时窗的方法的流程示意图。
图3为图2对应的算法伪代码。
图4为本发明实施例提供的地震数据振幅拟合有效时窗的系统的结构示意图。
具体实施方式
下面结合附图对本发明作进一步详细描述:
如图2、图3所示,本发明提供一种地震数据振幅拟合有效时窗的方法,包括:
S1、对道集数据按预设规则重新排列;
更具体地,把按照道号顺序排列的道集数据,按照预设规则重新排列,从而使每个样点所对应的所有道的数据连续存放在内存中。
所述预设规则为对同一时刻采集的样点赋予同一个样点号,样点号可以为自然数,并随着时间的增长而增长。例如将0.1ms采集到的样点的样点号记为1;将0.2ms采集到的样点的样点号记为2,将0.3ms采集到的样点的样点号记为2,依此类推。
在具体实施时,同一时刻也可为同一时间间隔,例如,将0ms至1ms采集到的样点的样点号记为1;将1ms至2ms采集到的样点的样点号记为2,将2ms至3ms采集到的样点的样点号记为3,依此类推。
可见,每一个样点号对应多个样点,但每个样点都有对应的道号。
S2、在排列后的道集数据中筛选出有效样点的排列;
可以使用计算机的向量计算方法一次性计算出在道集数据中有效样点的排列。具体实施时,可以包括步骤:201、使用向量指令SIMD计算同一个样点号对应的所有样点的振幅之和,获取所述样点号对应的振幅和计算值;202、去除所有振幅和计算值为零的样点号对应的样点,得到K个振幅和计算值不为零的样点号及其对应的样点,记为所述有效样点的排列。
S3、利用所述有效样点的排列计算有效起始道位置和终止道位置。
当某一样点的绝对值大于0,且与滤波函数的值与其绝对值相等,则该样点对应道号为有效起始道号;当某一样点的绝对值小于0,且其滤波函数的值与当前样点振幅和计算值相等,则该样点对应的道号为有效终止道号。具体地,使用fabs(data[样点号][道号])>1.0e-10&&滤波函数(data[样点号][道号])==fabs(data[样点号][道号])的条件找到相应样点处有效起始道号;同理,使用fabs(data[样点号][道号])<1.0e-10&&滤波函数(data[样点号][道号])==当前样点振幅和计算值的条件找到相应样点处有效终止道号。
在步骤S3之后,可以进一步包括:把相应的有效起始和终止道号,带入道拟合系数计算算法中,计算相应的拟合系数。
本发明提一种计算机可读存储介质,所述计算机可读存储介质存储有计算机可执行的至少一个程序,所述至少一个程序被所述计算机执行时使所述计算机执行上述任一实施例提供的方法中的步骤。
如图4所示,本发明提供一种地震数据振幅拟合有效时窗的系统,包括:预处理模块10、第一计算模块20、第二计算模块30。其中:
预处理模块10用于对道集数据按预设规则重新排列;所述预设规则为对同一时刻采集的样点赋予同一个样点号,样点号可以为自然数,并随着时间的增长而增长。例如将0.1ms采集到的样点的样点号记为1;将0.2ms采集到的样点的样点号记为2,将0.3ms采集到的样点的样点号记为2,依此类推。
在具体实施时,同一时刻也可为同一时间间隔,例如,将0ms至1ms采集到的样点的样点号记为1;将1ms至2ms采集到的样点的样点号记为2,将2ms至3ms采集到的样点的样点号记为3,依此类推。
可见,每一个样点号对应多个样点,但每个样点都有对应的道号。
第一计算模块20与预处理模块10相连,第一计算模块20用于在排列后的道集数据中筛选出有效样点的排列;可以使用计算机的向量计算方法一次性计算出在道集数据中有效样点的排列。具体实施时,可以包括步骤:201、使用向量指令SIMD计算同一个样点号对应的所有样点的振幅之和,获取所述样点号对应的振幅和计算值;202、去除所有振幅和计算值为零的样点号对应的样点,得到K个振幅和计算值不为零的样点号及其对应的样点,记为所述有效样点的排列。
第二计算模块30与第一计算模块20相连,第二计算模块30用于利用所述有效样点的排列计算有效起始道位置和终止道位置。当某一样点的绝对值大于0,且与滤波函数的值与其绝对值相等,则该样点对应道号为有效起始道号;当某一样点的绝对值小于0,且其滤波函数的值与当前样点振幅和计算值相等,则该样点对应的道号为有效终止道号。具体地,使用fabs(data[样点号][道号])>1.0e-10&&滤波函数(data[样点号][道号])==fabs(data[样点号][道号])的条件找到相应样点处有效起始道号;同理,使用fabs(data[样点号][道号])<1.0e-10&&滤波函数(data[样点号][道号])==当前样点振幅和计算值的条件找到相应样点处有效终止道号。
所述第二计算模块30还可以用于把相应的有效起始和终止道号,带入道拟合系数计算算法中,计算相应的拟合系数。
本发明提供的地震数据振幅拟合有效时窗的方法及其系统,利用道集数据重排技术和计算机的向量一次性计算方法,筛选出在整个道集数据中有效样点的排列,利用此有效样点排列,计算出有效起始道位置和终止道位置。从而道集数据振幅拟合有效时窗计算方法的算法复杂度提高到了O’,O’=N+K*M,其中K为有效的样点排列。在考虑最坏算法复杂度的情况下,算法效率提高了近1倍。
上述技术方案只是本发明的一种实施方式,对于本领域内的技术人员而言,在本发明公开了应用方法和原理的基础上,很容易做出各种类型的改进或变形,而不仅限于本发明上述具体实施方式所描述的方法,因此前面描述的方式只是优选的,而并不具有限制性的意义。
Claims (10)
1.一种地震数据振幅拟合有效时窗的方法,其特征在于,包括:
S1、对道集数据按预设规则重新排列;
S2、在排列后的道集数据中筛选出有效样点的排列;
S3、利用所述有效样点的排列计算有效起始道位置和终止道位置。
2.根据权利要求1所述地震数据振幅拟合有效时窗的方法,其特征在于,所述预设规则为对同一时刻采集的样点赋予同一个样点号。
3.根据权利要求2所述地震数据振幅拟合有效时窗的方法,其特征在于,所述步骤S2具体包括:
使用向量指令SIMD计算同一个样点号对应的所有样点的振幅之和,获取所述样点号对应的振幅和计算值;
去除所有振幅和计算值为零的样点号对应的样点,得到K个振幅和计算值不为零的样点号及其对应的样点,记为所述有效样点的排列。
4.根据权利要求1所述地震数据振幅拟合有效时窗的方法,其特征在于,所述步骤S2具体为:使用计算机的向量计算方法一次性计算出在道集数据中有效样点的排列。
5.根据权利要求1所述地震数据振幅拟合有效时窗的方法,其特征在于,所述步骤S3具体包括:
使用fabs(data[样点号][道号])>1.0e-10&&滤波函数(data[样点号][道号])==fabs(data[样点号][道号])的条件找到相应样点处有效起始道号;
使用fabs(data[样点号][道号])<1.0e-10&&滤波函数(data[样点号][道号])==当前样点振幅和计算值的条件找到相应样点处有效终止道号。
6.一种地震数据振幅拟合有效时窗的系统,其特征在于,包括:
预处理模块,用于对道集数据按预设规则重新排列;
第一计算模块,用于在排列后的道集数据中筛选出有效样点的排列;
第二计算模块,用于利用所述有效样点的排列计算有效起始道位置和终止道位置。
7.根据权利要求6所述地震数据振幅拟合有效时窗的系统,其特征在于,所述预设规则为对同一时刻采集的样点赋予同一个样点号。
8.根据权利要求7所述地震数据振幅拟合有效时窗的系统,其特征在于,所述第一计算模块具体用于:
使用向量指令SIMD计算同一个样点号对应的所有样点的振幅之和,获取所述样点号对应的振幅和计算值;
去除所有振幅和计算值为零的样点号对应的样点,得到K个振幅和计算值不为零的样点号及其对应的样点,记为所述有效样点的排列。
9.根据权利要求6所述地震数据振幅拟合有效时窗的系统,其特征在于,所述第二计算模块具体用于:
使用fabs(data[样点号][道号])>1.0e-10&&滤波函数(data[样点号][道号])==fabs(data[样点号][道号])的条件找到相应样点处有效起始道号;
使用fabs(data[样点号][道号])<1.0e-10&&滤波函数(data[样点号][道号])==当前样点振幅和计算值的条件找到相应样点处有效终止道号。
10.一种计算机可读存储介质,所述计算机可读存储介质存储有计算机可执行的至少一个程序,其特征在于,所述至少一个程序被所述计算机执行时使所述计算机执行上述权利要求1~5任一项所述的方法中的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910976991.XA CN112666601B (zh) | 2019-10-15 | 2019-10-15 | 一种地震数据振幅拟合有效时窗的方法及其系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910976991.XA CN112666601B (zh) | 2019-10-15 | 2019-10-15 | 一种地震数据振幅拟合有效时窗的方法及其系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112666601A true CN112666601A (zh) | 2021-04-16 |
CN112666601B CN112666601B (zh) | 2024-03-22 |
Family
ID=75399737
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910976991.XA Active CN112666601B (zh) | 2019-10-15 | 2019-10-15 | 一种地震数据振幅拟合有效时窗的方法及其系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112666601B (zh) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6148264A (en) * | 1998-07-06 | 2000-11-14 | Exxonmobil Upstream Research Company | Method for removing seismic noise caused by external activity |
US20050033518A1 (en) * | 2003-08-07 | 2005-02-10 | Jenner Edward Louis | Method for wavelet-based seismic amplitude inversion |
CN106249292A (zh) * | 2015-06-09 | 2016-12-21 | 中国石油化工股份有限公司 | 一种共反射点道集优化处理方法及系统 |
CN106338763A (zh) * | 2016-08-23 | 2017-01-18 | 中国石油天然气股份有限公司 | 一种确定地震剖面显示数据的方法及装置 |
CN107065007A (zh) * | 2017-05-24 | 2017-08-18 | 中国石油天然气股份有限公司 | 一种地震资料振幅调整方法和装置 |
CN107843919A (zh) * | 2017-10-27 | 2018-03-27 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种微地震监测数据振幅均衡方法 |
CN107918145A (zh) * | 2016-10-10 | 2018-04-17 | 中国石油化工股份有限公司 | 地震炮能量的并行化处理方法及系统 |
US20180149764A1 (en) * | 2016-04-01 | 2018-05-31 | Petrochina Company Limited | Method and apparatus for processing seismic data |
US20180275302A1 (en) * | 2017-03-27 | 2018-09-27 | Ion Geophysical Corporation | Amplitude compensation of reverse time migration (rtm) gathers for avo/ava analysis |
-
2019
- 2019-10-15 CN CN201910976991.XA patent/CN112666601B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6148264A (en) * | 1998-07-06 | 2000-11-14 | Exxonmobil Upstream Research Company | Method for removing seismic noise caused by external activity |
US20050033518A1 (en) * | 2003-08-07 | 2005-02-10 | Jenner Edward Louis | Method for wavelet-based seismic amplitude inversion |
CN106249292A (zh) * | 2015-06-09 | 2016-12-21 | 中国石油化工股份有限公司 | 一种共反射点道集优化处理方法及系统 |
US20180149764A1 (en) * | 2016-04-01 | 2018-05-31 | Petrochina Company Limited | Method and apparatus for processing seismic data |
CN106338763A (zh) * | 2016-08-23 | 2017-01-18 | 中国石油天然气股份有限公司 | 一种确定地震剖面显示数据的方法及装置 |
CN107918145A (zh) * | 2016-10-10 | 2018-04-17 | 中国石油化工股份有限公司 | 地震炮能量的并行化处理方法及系统 |
US20180275302A1 (en) * | 2017-03-27 | 2018-09-27 | Ion Geophysical Corporation | Amplitude compensation of reverse time migration (rtm) gathers for avo/ava analysis |
CN107065007A (zh) * | 2017-05-24 | 2017-08-18 | 中国石油天然气股份有限公司 | 一种地震资料振幅调整方法和装置 |
CN107843919A (zh) * | 2017-10-27 | 2018-03-27 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种微地震监测数据振幅均衡方法 |
Non-Patent Citations (6)
Title |
---|
E.STUCCHI等: "Seismic preprocessing and amplitude cross‐calibration for a time‐lapse amplitude study on seismic data from the Oseberg reservoir", GEOPHYSICAL PROSPECTING, vol. 53, no. 02, pages 265 - 282, XP055355287, DOI: 10.1111/j.1365-2478.2004.00471.x * |
夏洪瑞等: "常规二次多项式拟合地震数据", 石油物探, no. 05, pages 492 - 496 * |
徐钰;段卫星;徐维秀;曾维辉;徐雷良;: "高精度初至自动拾取综合方法研究", 物探与化探, vol. 34, no. 05, pages 595 - 599 * |
李云典: "利用时窗拟合法求取AVO属性剖面", 石油地球物理勘探, vol. 31, no. 03, pages 448 - 453 * |
杨尚琴;许自龙;洪承煜: "基于多线程的地震相干体属性提取算法", 计算机系统应用, vol. 21, no. 11, pages 72 - 75 * |
许璐;刘志成;王小青;谢金娥;贾春梅;: "基于结构中值滤波的CRP道集优化处理技术", 地球物理学进展, vol. 30, no. 04, pages 1804 - 1810 * |
Also Published As
Publication number | Publication date |
---|---|
CN112666601B (zh) | 2024-03-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110865357B (zh) | 一种基于参数优化vmd的激光雷达回波信号降噪方法 | |
CN111723329B (zh) | 一种基于全卷积神经网络的震相特征识别波形反演方法 | |
CN108549078B (zh) | 一种雷达脉冲信号跨信道合并及检测方法 | |
CN105259410A (zh) | 一种强噪声干扰下的欠采样波形的频率估计方法及其装置 | |
CN104463057B (zh) | 一种磁卡解码方法及装置 | |
CN107729291A (zh) | 一种基于分解技术的双线性模型参数辨识方法 | |
CN115937071A (zh) | 一种图像检测方法、装置、设备及介质 | |
CN112666601A (zh) | 一种地震数据振幅拟合有效时窗的方法及其系统 | |
CN109186571B (zh) | 一种陀螺仪滤波去噪方法 | |
CN106199693A (zh) | 地震数据速度谱自动拾取方法和装置 | |
CN111596358B (zh) | 一种多次波的压制方法、设备以及系统 | |
Forouzan et al. | Iterative median filtering for restoration of images with impulsive noise | |
CN102170292B (zh) | 一种数据处理方法、数据处理系统以及相关设备 | |
CN114236610B (zh) | 基于深度卷积网络先验的迭代地震资料解混方法及系统 | |
CN111968620A (zh) | 算法的测试方法、装置、电子设备及存储介质 | |
CN115470504A (zh) | 结合人工智能的数据风险分析方法及服务器 | |
CN113820749A (zh) | 一种基于机器学习的地震数据速度场异常反演方法 | |
CN109633744B (zh) | 地震子波的提取方法、装置、设备及存储介质 | |
CN112540408A (zh) | 一种基于深度学习的地震数据静校正处理方法及系统 | |
CN105301655A (zh) | 一种共成像点道集线性噪声去除方法及装置 | |
CN111354372A (zh) | 一种基于前后端联合训练的音频场景分类方法及系统 | |
CN115034972B (zh) | 图像去噪方法、装置及设备 | |
CN111722274B (zh) | 拾取地震波时距曲线的方法及装置 | |
CN110398773B (zh) | 一种针对部分缺失的地震数据的恢复重构方法 | |
CN117574056B (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 |