CN107843926A - 一种变偏vsp炮记录加密重构方法及*** - Google Patents
一种变偏vsp炮记录加密重构方法及*** Download PDFInfo
- Publication number
- CN107843926A CN107843926A CN201610829564.5A CN201610829564A CN107843926A CN 107843926 A CN107843926 A CN 107843926A CN 201610829564 A CN201610829564 A CN 201610829564A CN 107843926 A CN107843926 A CN 107843926A
- Authority
- CN
- China
- Prior art keywords
- encrypted
- big guns
- vsp
- depth
- recording
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 58
- 230000006835 compression Effects 0.000 claims abstract description 28
- 238000007906 compression Methods 0.000 claims abstract description 28
- 238000012937 correction Methods 0.000 claims abstract description 26
- 238000005070 sampling Methods 0.000 claims description 16
- 238000013461 design Methods 0.000 claims description 7
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000005259 measurement Methods 0.000 claims description 6
- 238000001228 spectrum Methods 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 3
- 230000001788 irregular Effects 0.000 description 3
- 230000008447 perception Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000001149 cognitive effect Effects 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012360 testing method 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/32—Transforming one recording into another or one representation into another
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
一种变偏VSP炮记录加密重构方法及***,所述方法包括以下步骤:步骤1:对变偏VSP炮记录进行高程校正和能量均衡;步骤2:将经过高程校正和能量均衡的变偏VSP炮记录转换成共深度点道集;步骤3:基于压缩感知方法对所述共深度点道集进行加密重构,获得加密重构道集;步骤4:将所述加密重构道集重新排列,获得加密重构变偏VSP炮记录。本发明可以更简单地对实际采集的变偏VSP炮记录进行加密重构,增加炮记录数,使观测***排列内数据点更密集、更均匀。
Description
技术领域
本发明涉及VSP勘探采集和处理技术领域,特别涉及一种变偏VSP炮记录加密重构方法及***。
背景技术
在地震勘探中,变偏VSP资料存在炮点少、炮间距大、采集成本高、覆盖次数少、成像精度低等问题。目前通过室内VSP数据处理,重构加密炮记录,使炮间距缩小,可以节省采集成本,加快采集周期,提高数据成像精度,改善资料横向分辨率;通过加密重构变偏VSP炮记录可以增加固定排列内炮数量,增加叠加次数,改善剖面质量。但是由于加密重构炮记录并不能通过两相邻炮的简单插值实现,尤其对于炮分布不规则、炮间距不满足Shannon/Nyquist采样定理的情况,加密重构炮记录更难实现,也未查阅到VSP炮记录间加密重构的相关文献。
近年来出现的压缩感知(Compressed sensing,简称CS)理论是一种介于数学和信息科学之间的全新信号采集、编解码理论。该理论突破了传统信号采集必须满足Shannon/Nyquist采样定理的限制,可以用远少于传统采样定理所需的采样点数或观测点数高质量重构出满足采样定理的原信号或图像。运用到地震勘探中,可用于解决不完整地震数据规则化问题,也为变偏VSP炮记录加密重构提供一种新手段。
基于现有技术手段,需要一种更简单的方法及***来实现对实际采集的变偏VSP炮记录进行加密重构,以增加炮记录数,使观测***排列内数据点更密集、更规则。
发明内容
本发明的目的是提供一种变偏VSP炮记录加密重构方法及***,其针对变偏VSP资料存在炮点少,炮间距大,采集成本高,覆盖次数少,成像精度低问题,应用新近出现的压缩感知理论,提出VSP炮记录加密重构方法,对实际采集的变偏VSP炮记录进行加密重构,增加炮记录数,使观测***排列内数据点更密集、更均匀。
本发明的一方面提供一种变偏VSP炮记录加密重构方法,包括以下步骤:
步骤1:对变偏VSP炮记录进行高程校正和能量均衡;
步骤2:将经过高程校正和能量均衡的变偏VSP炮记录转换成共深度点道集;
步骤3:基于压缩感知方法对所述共深度点道集进行加密重构,获得加密重构道集;
步骤4:将所述加密重构道集重新排列,获得加密重构变偏VSP炮记录。
优选地,所述步骤3中,通过在所述共深度点道集的两道之间增加至少一个重构道进行所述加密重构。
优选地,所述步骤3中,通过曲波变换实现所述压缩感知方法中的信号稀疏表达。
优选地,所述步骤3中,通过规则采样方法实现所述压缩感知方法中的随机测量矩阵的设计。
优选地,所述步骤3中,通过谱投影梯度算法在所述共深度点道集的两道之间增加至少一个重构道。
本发明的另一方面提供一种变偏VSP炮记录加密重构***,包括:
校正模块,用于对各实际变偏VSP炮记录进行高程校正和能量均衡;
转换模块,用于将经过高程校正和能量均衡的变偏VSP炮记录转换成共深度点道集;
感知重构模块,用于基于压缩感知方法对所述共深度点道集进行加密重构,获得加密重构道集;
排列模块,用于将所述加密重构道集重新排列,获得加密重构变偏VSP炮记录。
优选地,所述感知重构模块通过在所述共深度点道集的两道之间增加至少一个重构道进行所述加密重构。
优选地,所述感知重构模块通过曲波变换实现所述压缩感知方法中的信号稀疏表达。
优选地,所述感知重构模块通过规则采样方法实现所述压缩感知方法中的随机测量矩阵的设计。
优选地,所述感知重构模块通过谱投影梯度算法在所述共深度点道集的两道之间增加至少一个重构道。
本发明的有益效果在于:(1)所采用的压缩感知方法突破了传统的Shannon-Nyqiust采样定理的限制,用采集的少量地震炮/道数据重构完整数据,可用于解决不完整地震数据规则化问题;(2)对变偏VSP炮记录进行加密重构,增加炮记录数,使观测***排列内数据点更密集、更均匀,极大地节省采集时间与成本,增加了叠加次数;(3)通过室内数据规则化处理重构加密炮记录,使勘探区域满足密集且规则的覆盖次数,可以提高变偏VSP资料的横向分辨率。
附图说明
通过结合附图对本发明示例性实施例进行更详细的描述,本发明的上述以及其它目的、特征和优势将变得更加明显。
图1示出了根据示例性实施例的变偏VSP炮记录加密重构方法的流程图;
图2(a)示出了根据示例性实施例的观测***示意图,图2(b)示出了根据示例性实施例的地面距井口120m处采集炮Z分量记录;
图3(a)示出了根据示例性实施例的采集炮在深度930m处点Z分量高程校正前道集,图3(b)为采集炮在深度930m处点Z分量高程校正后道集;
图4(a)示出了根据示例性实施例的采集炮在深度930m处点Z分量道集,图4(b)为采集炮在深度740m处点Z分量道集,图4(c)为采集炮在深度540m处点Z分量道集,图4(d)为采集炮在深度340m处点Z分量道集;
图5(a)示出了根据示例性实施例的采集炮在深度740m处点Z分量能量均衡前道集,图5(b)为采集炮在深度740m处点Z分量能量均衡后道集;
图6(a)示出了根据示例性实施例的采集炮在深度930m处点Z分量加密重构后的道集,图6(b)为采集炮在深度780m处点Z分量加密重构后的道集,图6(c)为采集炮在深度740m处点Z分量加密重构后的道集;
图7(a)示出了根据示例性实施例的地面距井口40m处采集炮Z分量记录,图7(b)为地面距井口60m处通过压缩感知方法重构炮Z分量记录,图7(c)为地面距井口80m处采集炮Z分量记录。
具体实施方式
下面将参照附图更详细地描述本发明的优选实施例。虽然附图中显示了本发明的优选实施例,然而应该理解,可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反,提供这些实施例是为了使本发明更加透彻和完整,并且能够将本发明的范围完整地传达给本领域的技术人员。
在地震勘探中,变偏VSP资料存在炮点少、炮间距大、采集成本高、覆盖次数少、成像精度低等问题。目前通过室内VSP数据处理,重构加密炮记录,使炮间距缩小,可以节省采集成本,加快采集周期,提高数据成像精度,改善资料横向分辨率;通过加密重构变偏VSP炮记录可以增加固定排列内炮数量,增加叠加次数,改善剖面质量。但是由于加密重构炮记录并不能通过两相邻炮的简单插值实现,尤其对于炮分布不规则、炮间距不满足Shannon/Nyquist采样定理的情况,加密重构炮记录更难实现,也未查阅到VSP炮记录间加密重构的相关文献。
针对上述问题,本发明一方面提出一种变偏VSP炮记录加密重构方法,应用新近出现的压缩感知理论,提出VSP炮记录加密重构方法,对实际采集的变偏VSP炮记录进行加密重构,增加炮记录数,使观测***排列内数据点更密集、更规则。
图1示出了根据示例性实施例的变偏VSP炮记录加密重构方法的流程图,其包括以下步骤:
步骤1:对变偏VSP炮记录进行高程校正和能量均衡。
可应用现有软件对变偏VSP炮记录进行高程校正和能量均衡。通过高程校正可以将炮点位置校正到同一海平面上,消除炮点受地表起伏影响产生的高程差,从而可以解决道集波形起跳点不同,波组特征不一致,相邻数据相关性差的问题,便于炮记录加密重构。通过能量均衡可以消除每一次放炮引起的能量差异,使相邻炮能量保持在一个合理的数量级上,便于利用相邻道重构中间插值道。进行能量均衡可以解决不同共深度道集上能量差异明显,影响重构效果的问题,使能量显示均衡,便于加密重构。
步骤2:将经过高程校正和能量均衡的变偏VSP炮记录转换成共深度点道集,以便于进行加密重构。
步骤3:基于压缩感知方法对所述共深度点道集进行加密重构。
具体来说,可以通过曲波变换实现压缩感知方法中的信号稀疏表达,通过规则采样方法实现压缩感知方法中的随机测量矩阵的设计,通过基于谱投影梯度算法实现压缩感知方法中的高效恢复算法,进而在共深度点道集的两道之间增加至少一个重构道,获取经过加密的重构道集。
步骤4:将加密重构道集重新排列,获得加密重构变偏VSP炮记录,可以得到更密集、更规则的炮点记录。
经过步骤1至3,已经对共深度点道集进行了加密重构,在本步骤中,对加密重构道集进行分选,重新抽取成各炮点记录,从而获得加密重构变偏VSP炮记录。
进一步地,在步骤3中,可以通过在共深度点道集的两道之间增加至少一个重构道进行加密重构,转换成炮记录就是两炮之间增加了至少一炮,即增加了炮密度。
进一步地,在步骤3中,可以通过曲波变换实现压缩感知方法中的信号稀疏表达,也可以通过如复氏变换等方式实现,可以根据实际情况选取合适的方法。
进一步地,在步骤3中,压缩感知方法可以采用规则采样方法实现随机测量矩阵的设计,也可以采用不规则随机采样等方法实现,可以根据实际情况选取合适的方法。
进一步地,在步骤3中,通过基于谱投影梯度算法的压缩感知方法对共深度点道集进行加密重构,获得加密重构道集,也可以根据实际情况选取其他快速重构算法。
根据上述方法,可以对实际变偏VSP炮记录进行加密重构处理,实现了变偏VSP炮记录加密重构,获取了加密重构炮记录,极大地节省采集时间与成本,增加了叠加次数,改善资料横向分辨率,为获取均匀、密集的变偏VSP炮记录提出了一种高效的室内处理方法。
本发明还提供一种变偏VSP炮记录加密重构***,根据示例性实施例的变偏VSP炮记录加密重构***包括:
校正模块,用于对各实际变偏VSP炮记录进行高程校正和能量均衡;
转换模块,用于将经过高程校正和能量均衡的变偏VSP炮记录转换成共深度点道集;
感知重构模块,用于基于压缩感知方法对所述共深度点道集进行加密重构,获得加密重构道集;
排列模块,用于将所述加密重构道集重新排列,获得加密重构变偏VSP炮记录。
实施例
图2(a)示出了根据示例性实施例的观测***示意图,图2(b)示出了根据示例性实施例的地面距井口120m处采集炮Z分量记录。
在该实施例中,实际变偏VSP资料地面放25炮,井中埋置64个检波器。炮间距40m,检波点距10m,最近的检波点距地面300m,炮记录长度2s,采样间距0.5ms,采样点数4001个。
由于VSP资料炮点位于地面,各炮记录高程不同且各炮能量不均衡,在做炮记录加密前需要进行高程校正和能量均衡,并将其转换成共深度点道集。
图3(a)为采集炮在深度930m处点Z分量高程校正前道集,图3(b)为采集炮在深度930m处点Z分量高程校正后道集。对比图3(a)与图3(b)可以看出,高程校正前道集波形起跳点不同,波组特征不一致,相邻数据相关性差,高程校正后道集波形起跳点一致,波组特征一致性、相关性得到明显改善,便于炮记录加密重构。
图4(a)~图4(d)分别显示应用示例中采集炮在深度930m、740m、540m和340m处点Z分量道集。通过图4(a)~图4(d)可以看出,不同共深度道集上能量差异明显,影响重构效果,因此对数据做均衡处理,使能量显示均衡,便于加密重构。
图5(a)和图5(b)分别显示实施例中的采集炮在深度740m处点Z分量能量均衡前道集和采集炮在深度740m处点Z分量能量均衡后道集,通过对比图形可以看出,均衡后的道集更均衡,能量差异已不明显,可以使重构效果更好。
在该实施例中,对各深度点道集应用压缩感知方法进行加密重构,其中,通过曲波变换实现压缩感知方法中的信号稀疏表达,通过规则采样方法实现压缩感知方法中的随机测量矩阵的设计,通过基于谱投影梯度算法实现压缩感知方法中的高效恢复算法,对共深度点道集进行加密重构。
在共深度点道集的2道之间增加1个重构道,转换成炮记录就是2炮之间的增加1炮,部分加密重构后的共深度点道集Z分量处理结果如图6(a)~图6(c)所示,其中,图6(a)为采集炮在深度930m处点Z分量加密重构后的道集,图6(b)为采集炮在深度780m处点Z分量加密重构后的道集,图6(c)为采集炮在深度740m处点Z分量加密重构后的道集。
将各深度点加密重构道集重新排列组成重构炮记录,并获得加密重构变偏VSP炮记录。图7(a)为地面距井口40m处采集炮Z分量记录,图7(b)为地面距井口60m处通过压缩感知方法重构炮Z分量记录,图7(c)为地面距井口80m处采集炮Z分量记录。
本发明的实际数据试验结果表明,根据示例性实施例的变偏VSP炮记录加密重构方法,有效解决了不同共深度道集上道集波形起跳点不同,波组特征不一致,相邻数据相关性差,能量差异明显等因素对加密重构的影响,实现了变偏VSP炮记录加密重构,极大地节省采集时间与成本,增加了叠加次数,改善资料横向分辨率,为获取均匀、密集的变偏VSP炮记录提出了一种高效的室内处理方法。
上述技术方案只是本发明的一种实施例,对于本领域内的技术人员而言,在本发明公开的原理的基础上,很容易做出各种类型的改进或变形,而不仅限于本发明上述具体实施例的描述,因此前面的描述只是优选的,而并不具有限制性的意义。
Claims (10)
1.一种变偏VSP炮记录加密重构方法,包括以下步骤:
步骤1:对变偏VSP炮记录进行高程校正和能量均衡;
步骤2:将经过高程校正和能量均衡的变偏VSP炮记录转换成共深度点道集;
步骤3:基于压缩感知方法对所述共深度点道集进行加密重构,获得加密重构道集;
步骤4:将所述加密重构道集重新排列,获得加密重构变偏VSP炮记录。
2.根据权利要求1所述的变偏VSP炮记录加密重构方法,其中,所述步骤3中,通过在所述共深度点道集的两道之间增加至少一个重构道进行所述加密重构。
3.根据权利要求1所述的变偏VSP炮记录加密重构方法,其中,所述步骤3中,通过曲波变换实现所述压缩感知方法中的信号稀疏表达。
4.根据权利要求1所述的变偏VSP炮记录加密重构方法,其中,所述步骤3中,通过规则采样方法实现所述压缩感知方法中的随机测量矩阵的设计。
5.根据权利要求2所述的变偏VSP炮记录加密重构方法,其中,所述步骤3中,通过谱投影梯度算法在所述共深度点道集的两道之间增加至少一个重构道。
6.一种变偏VSP炮记录加密重构***,包括:
校正模块,用于对各实际变偏VSP炮记录进行高程校正和能量均衡;
转换模块,用于将经过高程校正和能量均衡的变偏VSP炮记录转换成共深度点道集;
感知重构模块,用于基于压缩感知方法对所述共深度点道集进行加密重构,获得加密重构道集;
排列模块,用于将所述加密重构道集重新排列,获得加密重构变偏VSP炮记录。
7.根据权利要求6所述的变偏VSP炮记录加密重构***,其中,所述感知重构模块通过在所述共深度点道集的两道之间增加至少一个重构道进行所述加密重构。
8.根据权利要求6所述的变偏VSP炮记录加密重构***,其中,所述感知重构模块通过曲波变换实现所述压缩感知方法中的信号稀疏表达。
9.根据权利要求6所述的变偏VSP炮记录加密重构***,其中,所述感知重构模块通过规则采样方法实现所述压缩感知方法中的随机测量矩阵的设计。
10.根据权利要求6所述的变偏VSP炮记录加密重构***,其中,所述感知重构模块通过谱投影梯度算法在所述共深度点道集的两道之间增加至少一个重构道。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610829564.5A CN107843926A (zh) | 2016-09-18 | 2016-09-18 | 一种变偏vsp炮记录加密重构方法及*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610829564.5A CN107843926A (zh) | 2016-09-18 | 2016-09-18 | 一种变偏vsp炮记录加密重构方法及*** |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107843926A true CN107843926A (zh) | 2018-03-27 |
Family
ID=61656881
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610829564.5A Pending CN107843926A (zh) | 2016-09-18 | 2016-09-18 | 一种变偏vsp炮记录加密重构方法及*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107843926A (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140035582A1 (en) * | 2011-04-21 | 2014-02-06 | Koninklijke Philips N.V. | Magnetic resonance imaging of object in motion |
CN104345340A (zh) * | 2013-07-29 | 2015-02-11 | 中国石油化工股份有限公司 | 一种地震数据规则化中抖动稀疏采样方法 |
CN104422954A (zh) * | 2013-08-20 | 2015-03-18 | 中国石油化工股份有限公司 | 一种地震数据规则化中稀疏采样方法 |
US20150123659A1 (en) * | 2013-11-07 | 2015-05-07 | Sebastian Weingartner | System and method for free-breathing volumetric imaging of cardiac tissue |
-
2016
- 2016-09-18 CN CN201610829564.5A patent/CN107843926A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140035582A1 (en) * | 2011-04-21 | 2014-02-06 | Koninklijke Philips N.V. | Magnetic resonance imaging of object in motion |
CN104345340A (zh) * | 2013-07-29 | 2015-02-11 | 中国石油化工股份有限公司 | 一种地震数据规则化中抖动稀疏采样方法 |
CN104422954A (zh) * | 2013-08-20 | 2015-03-18 | 中国石油化工股份有限公司 | 一种地震数据规则化中稀疏采样方法 |
US20150123659A1 (en) * | 2013-11-07 | 2015-05-07 | Sebastian Weingartner | System and method for free-breathing volumetric imaging of cardiac tissue |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106291678B (zh) | 一种地震数据采集方法及*** | |
CN104345340B (zh) | 一种地震数据规则化中抖动稀疏采样方法 | |
CN106249291A (zh) | 一种基于二维非均匀曲波变换的高精度地震数据重建方法 | |
AU2015229106B2 (en) | Spatial sampling improvements and the formation of arrays using popcorn and simultaneous source acquisition | |
CN110361778A (zh) | 一种基于生成对抗网络的地震数据重建方法 | |
CN106970419A (zh) | 一种基于线性Bregman算法的非均匀曲波三维地震数据重建方法 | |
Koulakov et al. | Seismic tomography of volcanoes | |
CN108594302A (zh) | 一种地震子波的提取方法及处理终端 | |
CN104597499B (zh) | 可控震源独立同步激发地震数据邻炮干扰压制方法和装置 | |
CN102854532A (zh) | 三维叠前炮检域随机噪声压制方法 | |
Guo et al. | Compressed ghost edge imaging | |
CN107843926A (zh) | 一种变偏vsp炮记录加密重构方法及*** | |
CN106908836A (zh) | 采集脚印压制方法及*** | |
CN108376245B (zh) | 基于ud通道的时空序列图像震源识别方法 | |
Wright et al. | The interpretation of ionospheric radio drift measurements—II Kinesonde observations of microstructure and vertical motion in sporadic E | |
CN100453044C (zh) | 基于原始投影正弦图的ct射束硬化校正方法 | |
CN113484906B (zh) | 一种基于低频能谱数据驱动的高分辨率能量场重建方法 | |
CN109541687B (zh) | 一种熵约束数据驱动正规框架地震数据规则化方法 | |
CN104880732B (zh) | 一种十字子集的构建方法及装置 | |
Moldoveanu et al. | Quad coil towed streamer marine acquisition | |
El-Emam et al. | Maximizing the value of Single-Sensor measurements, Kuwait experience | |
CN106814395B (zh) | 确定背景炮的方法和装置 | |
CN104166160B (zh) | 一种地震数据组合方法、装置及*** | |
CN108844616A (zh) | 一种次声事件检测方法及*** | |
CN116736372B (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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180327 |
|
RJ01 | Rejection of invention patent application after publication |