CN112649864A - 地震剖面变密度分段色标显示方法及系统 - Google Patents

地震剖面变密度分段色标显示方法及系统 Download PDF

Info

Publication number
CN112649864A
CN112649864A CN201910968681.3A CN201910968681A CN112649864A CN 112649864 A CN112649864 A CN 112649864A CN 201910968681 A CN201910968681 A CN 201910968681A CN 112649864 A CN112649864 A CN 112649864A
Authority
CN
China
Prior art keywords
amplitude
seismic
seismic data
color scale
normalized
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.)
Pending
Application number
CN201910968681.3A
Other languages
English (en)
Inventor
刘旭跃
张兵
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Chemical Corp
Sinopec Geophysical Research Institute
Original Assignee
China Petroleum and Chemical Corp
Sinopec Geophysical Research Institute
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Petroleum and Chemical Corp, Sinopec Geophysical Research Institute filed Critical China Petroleum and Chemical Corp
Priority to CN201910968681.3A priority Critical patent/CN112649864A/zh
Publication of CN112649864A publication Critical patent/CN112649864A/zh
Pending legal-status Critical Current

Links

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/30—Analysis
    • G01V1/301—Analysis for determining seismic cross-sections or geostructures

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)计算归一化地震数据样点的均方根振幅：
Figure BDA0002231349320000021
其中，AMPrms为均方根振幅，t1为时窗初始时间，R是时窗半径，AMPi是样点的振幅值。
优选地，通过公式(2)计算所述剩余补偿因子：
Figure BDA0002231349320000022
其中，F为剩余补偿因子，ns＝采样点总数-当前采样点下标+1，Δdt为采样间隔。
优选地，根据所述补偿后振幅确定色标包括：确定振幅值与所述色标的对应关系；将每一个采样点的补偿振幅值代入所述对应关系，标定所述色标。
根据本发明的另一方面，提出了一种地震剖面变密度分段色标显示系统，其特征在于，该系统包括：存储器，存储有计算机可执行指令；处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：读取地震数据，针对所述地震数据进行归一化处理，获得归一化地震数据；针对所述归一化地震数据的振幅进行能量补偿，获得补偿后振幅；根据所述补偿后振幅确定色标；根据所述归一化地震数据与所述色标，绘制地震剖面图。
优选地，通过以下步骤进行所述能量补偿：确定时窗，设定时窗的起始时间和时窗半径，进而计算归一化地震数据样点的均方根振幅；根据所述均方根振幅，计算剩余补偿因子；根据所述剩余补偿因子在道集内进行叠加，获得所述补偿后振幅。
优选地，通过公式(1)计算归一化地震数据样点的均方根振幅：
Figure BDA0002231349320000031
其中，AMPrms为均方根振幅，t1为时窗初始时间，R是时窗半径，AMPi是样点的振幅值。
优选地，通过公式(2)计算所述剩余补偿因子：
Figure BDA0002231349320000032
其中，F为剩余补偿因子，ns＝采样点总数-当前采样点下标+1，Δdt为采样间隔。
优选地，根据所述补偿后振幅确定色标包括：确定振幅值与所述色标的对应关系；将每一个采样点的补偿振幅值代入所述对应关系，标定所述色标。
本发明的方法和装置具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的，或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述，这些附图和具体实施方式共同用于解释本发明的特定原理。
附图说明
通过结合附图对本发明示例性实施例进行更详细的描述，本发明的上述以及其它目的、特征和优势将变得更加明显，其中，在本发明示例性实施例中，相同的参考标号通常代表相同部件。
图1示出了根据本发明的地震剖面变密度分段色标显示方法的步骤的流程图。
图2示出了根据本发明的一个实施例的色标的示意图。
图3示出了根据本发明的一个实施例的地震剖面图。
图4示出了根据本发明的一个实施例的能量补偿后获得的地震剖面图。
具体实施方式
下面将参照附图更详细地描述本发明。虽然附图中显示了本发明的优选实施例，然而应该理解，可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了使本发明更加透彻和完整，并且能够将本发明的范围完整地传达给本领域的技术人员。
图1示出了根据本发明的地震剖面变密度分段色标显示方法的步骤的流程图。
在该实施例中，根据本发明的地震剖面变密度分段色标显示方法可以包括：步骤101，读取地震数据，针对地震数据进行归一化处理，获得归一化地震数据；步骤102，针对归一化地震数据的振幅进行能量补偿，获得补偿后振幅；步骤103，根据补偿后振幅确定色标；步骤104，根据归一化地震数据与色标，绘制地震剖面图。
在一个示例中，通过以下步骤进行能量补偿：确定时窗，设定时窗的起始时间和时窗半径，进而计算归一化地震数据样点的均方根振幅；根据均方根振幅，计算剩余补偿因子；根据剩余补偿因子在道集内进行叠加，获得补偿后振幅。
在一个示例中，通过公式(1)计算归一化地震数据样点的均方根振幅：
Figure BDA0002231349320000041
其中，AMPrms为均方根振幅，t1为时窗初始时间，R是时窗半径，AMPi是样点的振幅值。
在一个示例中，通过公式(2)计算剩余补偿因子：
Figure BDA0002231349320000051
其中，F为剩余补偿因子，ns＝采样点总数-当前采样点下标+1，Δdt为采样间隔。
在一个示例中，根据补偿后振幅确定色标包括：确定振幅值与色标的对应关系；将每一个采样点的补偿振幅值代入对应关系，标定色标。
具体地，根据本发明的地震剖面变密度分段色标显示方法可以包括：
读取地震数据，当前窗口显示的道数是tracenum，地震采样点数是sample，将地震数据归一化，从-1到1，存储在数组seis_data[tracenum*sample]中，获得归一化地震数据。
确定时窗，设定时窗的起始时间和时窗半径，进而通过公式(1)计算归一化地震数据样点的均方根振幅；由于地震信号在地层传播中会随着深度增加而衰减，所以剩余补偿因子F采用线性值，根据均方根振幅，通过公式(2)计算剩余补偿因子；根据剩余补偿因子在CMP道集内进行叠加，获得补偿后振幅。
色标colorTab有n种颜色，确定振幅值与色标的对应关系，Min(n)对应补偿后振幅最小值，Max(n)对应补偿后振幅最大值；将每一个采样点的补偿振幅值代入对应关系，标定色标；将归一化地震数据内的所有补偿振幅值根据色标确定颜色，绘制地震剖面图。
本方法通过分段色标显示，形象高效地绘制地震剖面图，通过地震能量补偿，增强图像对比性，使地震剖面图显示更加清晰明了，有利于地震处理解释工作的顺利开展，提高工作效率。
应用示例
为便于理解本发明实施例的方案及其效果，以下给出一个具体应用示例。本领域技术人员应理解，该示例仅为了便于理解本发明，其任何具体细节并非意在以任何方式限制本发明。
根据本发明的地震剖面变密度分段色标显示方法包括：
读取地震数据，当前窗口显示的道数是tracenum，地震采样点数是sample，将地震数据归一化，从-1到1，存储在数组seis_data[tracenum*sample]中，获得归一化地震数据。
确定时窗，设定时窗的起始时间和时窗半径，进而通过公式(1)计算归一化地震数据样点的均方根振幅；由于地震信号在地层传播中会随着深度增加而衰减，所以剩余补偿因子F采用线性值，根据均方根振幅，通过公式(2)计算剩余补偿因子；根据剩余补偿因子在CMP道集内进行叠加，获得补偿后振幅。
图2示出了根据本发明的一个实施例的色标的示意图。
色标colorTab有n种颜色，屏幕像素数组为SubBit[nWidth*nHeight*4]，用RGB色彩组合表示颜色，包括颜色值R、颜色值G、颜色值B，则色标如图2所示。确定振幅值与色标的对应关系，Min(n)对应补偿后振幅最小值，Max(n)对应补偿后振幅最大值；将每一个采样点的补偿振幅值代入对应关系，标定色标；将归一化地震数据内的所有补偿振幅值根据色标确定颜色，绘制地震剖面图。
图3示出了根据本发明的一个实施例的未经能量补偿获得的地震剖面图。
图4示出了根据本发明的一个实施例的能量补偿后获得的地震剖面图。
二者对比可知，经过能量补偿后获得的地震剖面图更加清晰，显示的数据更多。
综上所述，本发明通过分段色标显示，形象高效地绘制地震剖面图，通过地震能量补偿，增强图像对比性，使地震剖面图显示更加清晰明了，有利于地震处理解释工作的顺利开展，提高工作效率。
本领域技术人员应理解，上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果，并不意在将本发明的实施例限制于所给出的任何示例。
根据本发明的实施例，提供了一种地震剖面变密度分段色标显示系统，其特征在于，该系统包括：存储器，存储有计算机可执行指令；处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：读取地震数据，针对地震数据进行归一化处理，获得归一化地震数据；针对归一化地震数据的振幅进行能量补偿，获得补偿后振幅；根据补偿后振幅确定色标；根据归一化地震数据与色标，绘制地震剖面图。
在一个示例中，通过以下步骤进行能量补偿：确定时窗，设定时窗的起始时间和时窗半径，进而计算归一化地震数据样点的均方根振幅；根据均方根振幅，计算剩余补偿因子；根据剩余补偿因子在道集内进行叠加，获得补偿后振幅。
在一个示例中，通过公式(1)计算归一化地震数据样点的均方根振幅：
Figure BDA0002231349320000071
其中，AMPrms为均方根振幅，t1为时窗初始时间，R是时窗半径，AMPi是样点的振幅值。
在一个示例中，通过公式(2)计算剩余补偿因子：
Figure BDA0002231349320000072
其中，F为剩余补偿因子，ns＝采样点总数-当前采样点下标+1，Δdt为采样间隔。
在一个示例中，根据补偿后振幅确定色标包括：确定振幅值与色标的对应关系；将每一个采样点的补偿振幅值代入对应关系，标定色标。
本系统通过分段色标显示，形象高效地绘制地震剖面图，通过地震能量补偿，增强图像对比性，使地震剖面图显示更加清晰明了，有利于地震处理解释工作的顺利开展，提高工作效率。
本领域技术人员应理解，上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果，并不意在将本发明的实施例限制于所给出的任何示例。
以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。

Claims (10)

1.一种地震剖面变密度分段色标显示方法，其特征在于，包括：
读取地震数据，针对所述地震数据进行归一化处理，获得归一化地震数据；
针对所述归一化地震数据的振幅进行能量补偿，获得补偿后振幅；
根据所述补偿后振幅确定色标；
根据所述归一化地震数据与所述色标，绘制地震剖面图。
2.根据权利要求1所述的地震剖面变密度分段色标显示方法，其中，通过以下步骤进行所述能量补偿：
确定时窗，设定时窗的起始时间和时窗半径，进而计算归一化地震数据样点的均方根振幅；
根据所述均方根振幅，计算剩余补偿因子；
根据所述剩余补偿因子在道集内进行叠加，获得所述补偿后振幅。
3.根据权利要求2所述的地震剖面变密度分段色标显示方法，其中，通过公式(1)计算归一化地震数据样点的均方根振幅：
Figure FDA0002231349310000011
其中，AMPrms为均方根振幅，t1为时窗初始时间，R是时窗半径，AMPi是样点的振幅值。
4.根据权利要求2所述的地震剖面变密度分段色标显示方法，其中，通过公式(2)计算所述剩余补偿因子：
Figure FDA0002231349310000012
其中，F为剩余补偿因子，ns＝采样点总数-当前采样点下标+1，Δdt为采样间隔。
5.根据权利要求1所述的地震剖面变密度分段色标显示方法，其中，根据所述补偿后振幅确定色标包括：
确定振幅值与所述色标的对应关系；
将每一个采样点的补偿振幅值代入所述对应关系，标定所述色标。
6.一种地震剖面变密度分段色标显示系统，其特征在于，该系统包括：
存储器，存储有计算机可执行指令；
处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：
读取地震数据，针对所述地震数据进行归一化处理，获得归一化地震数据；
针对所述归一化地震数据的振幅进行能量补偿，获得补偿后振幅；
根据所述补偿后振幅确定色标；
根据所述归一化地震数据与所述色标，绘制地震剖面图。
7.根据权利要求6所述的地震剖面变密度分段色标显示系统，其中，通过以下步骤进行所述能量补偿：
确定时窗，设定时窗的起始时间和时窗半径，进而计算归一化地震数据样点的均方根振幅；
根据所述均方根振幅，计算剩余补偿因子；
根据所述剩余补偿因子在道集内进行叠加，获得所述补偿后振幅。
8.根据权利要求7所述的地震剖面变密度分段色标显示系统，其中，通过公式(1)计算归一化地震数据样点的均方根振幅：
Figure FDA0002231349310000031
其中，AMPrms为均方根振幅，t1为时窗初始时间，R是时窗半径，AMPi是样点的振幅值。
9.根据权利要求7所述的地震剖面变密度分段色标显示系统，其中，通过公式(2)计算所述剩余补偿因子：
Figure FDA0002231349310000032
其中，F为剩余补偿因子，ns＝采样点总数-当前采样点下标+1，Δdt为采样间隔。
10.根据权利要求6所述的地震剖面变密度分段色标显示系统，其中，根据所述补偿后振幅确定色标包括：
确定振幅值与所述色标的对应关系；
将每一个采样点的补偿振幅值代入所述对应关系，标定所述色标。
CN201910968681.3A 2019-10-12 2019-10-12 地震剖面变密度分段色标显示方法及系统 Pending CN112649864A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910968681.3A CN112649864A (zh) 2019-10-12 2019-10-12 地震剖面变密度分段色标显示方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910968681.3A CN112649864A (zh) 2019-10-12 2019-10-12 地震剖面变密度分段色标显示方法及系统

Publications (1)

Publication Number Publication Date
CN112649864A true CN112649864A (zh) 2021-04-13

Family

ID=75343733

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910968681.3A Pending CN112649864A (zh) 2019-10-12 2019-10-12 地震剖面变密度分段色标显示方法及系统

Country Status (1)

Country Link
CN (1) CN112649864A (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104502965A (zh) * 2014-12-22 2015-04-08 中国石油天然气集团公司 一种振幅补偿因子的反演方法
WO2016141598A1 (zh) * 2015-03-12 2016-09-15 中国矿业大学(北京) 共姿态道集地表一致性振幅补偿方法及装置
CN108181649A (zh) * 2017-11-23 2018-06-19 中国石油天然气集团公司 一种介质分界面振幅补偿方法及装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104502965A (zh) * 2014-12-22 2015-04-08 中国石油天然气集团公司 一种振幅补偿因子的反演方法
WO2016141598A1 (zh) * 2015-03-12 2016-09-15 中国矿业大学(北京) 共姿态道集地表一致性振幅补偿方法及装置
CN108181649A (zh) * 2017-11-23 2018-06-19 中国石油天然气集团公司 一种介质分界面振幅补偿方法及装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
周文辉等: "SEGY 格式地震数据的三维可视化", 计算机应用与软件, vol. 36, no. 2, pages 79 - 81 *
王德安: "剩余能量补偿及面元均化技术在吉林探区的应用", 化工管理, pages 76 *

Similar Documents

Publication Publication Date Title
US20060165311A1 (en) Spatial standard observer
US11010883B2 (en) Automated analysis of petrographic thin section images using advanced machine learning techniques
CN103278864B (zh) 洞缝型储层的地质特征参数及分布的测定方法及装置
US20080112614A1 (en) Histogram tile map for gpu based histogram computation
US20240169505A1 (en) Optical fiber image calibration method and apparatus, and computer- readable storage medium
CN105447879A (zh) 乳房图像中检测胸肌的方法及装置
US20100027848A1 (en) Structure-independent analysis of 3-D seismic random noise
CN103955906A (zh) 基于蝙蝠算法的Criminisi图像修复方法
CN116152209A (zh) 一种耳机罩缺陷检测方法、装置、设备及存储介质
CN114841927A (zh) 一种页岩储层裂缝识别方法、装置、设备及存储介质
CN113255607A (zh) 基于多模型融合的遥感房屋检测方法及装置
CN115683431A (zh) 基于直线跟踪算法的拉索索力确定方法、装置及设备
CN115343685A (zh) 应用于病害识别的多维度探地雷达检测方法、装置、设备
EP2536123A1 (en) Image processing method and image processing apparatus
CN112649864A (zh) 地震剖面变密度分段色标显示方法及系统
JP7413049B2 (ja) 地中レーダーのデータ処理方法、データ処理プログラム及び地中レーダー装置
CN104463898A (zh) 一种物流平台用图像识别计数系统
CN111142161B (zh) 基于地震数据的复数域地质影像方法及电子设备
EP2799893B1 (en) Adaptive Frequency-Domain Windowing
GB2599762A (en) Fault skeletonization for fault identification in a subterranean environment
CN107589456B (zh) 一种获取地震数据的方法和装置及计算机可读存储介质
CN105241886A (zh) 一种基于hvs与引导滤波器的土遗址裂缝检测方法
CN110737024B (zh) 地震波信号的处理方法、装置、设备及存储介质
CN112417177B (zh) 一种基于调研数据的知识图谱构建方法
Titarenko et al. Study of the ability of neural networks to extract and use semantic information when they are trained to reconstruct noisy images

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