CN107367758B - 一种提高信噪比的可控震源加权相关的方法 - Google Patents
一种提高信噪比的可控震源加权相关的方法 Download PDFInfo
- Publication number
- CN107367758B CN107367758B CN201610309760.XA CN201610309760A CN107367758B CN 107367758 B CN107367758 B CN 107367758B CN 201610309760 A CN201610309760 A CN 201610309760A CN 107367758 B CN107367758 B CN 107367758B
- Authority
- CN
- China
- Prior art keywords
- dominant frequency
- scanning signal
- wavelet
- slope
- signal
- 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
- 238000000034 method Methods 0.000 title claims abstract description 44
- 238000001228 spectrum Methods 0.000 claims abstract description 28
- 238000013461 design Methods 0.000 claims description 5
- 238000010183 spectrum analysis Methods 0.000 claims description 5
- 238000007796 conventional method Methods 0.000 abstract description 4
- 230000008569 process Effects 0.000 description 8
- 230000000694 effects Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 2
- 230000002596 correlated effect Effects 0.000 description 2
- 230000001737 promoting effect Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 239000002360 explosive Substances 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
- 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/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
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)
- Complex Calculations (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种提高信噪比的可控震源加权相关的方法，该方法包括以下步骤：获得可控震源单炮的主频和主频部分的斜率；改造参考信号的频谱，提升主频部分的能量，得到加权后的频谱曲线；将该加权后的频谱曲线转化为加权扫描信号；采用该加权扫描信号与原始记录进行相关，得到加权相关后的可控震源单炮记录。该方法通过在相关前提高参考信号的主频部分能量，转化为加权扫描信号，再将加权扫描信号与原始母记录进行相关得到单炮记录，使用该方法得到的单炮记录较常规方法得到的单炮记录信噪比更高。
Description
技术领域
本发明涉及油气地球物理勘探，更具体地，涉及一种提高信噪比的可控震源加权相关的方法。
背景技术
可控震源在目前的地震采集中应用越来越广泛，尤其是目前西部地区对炸药使用的敏感性，更加需要扩大可控震源的使用。可控震源不同于炸药震源激发，需要在地表激发，这就使得可控震源的地震波传播能量两次受到近地表的影响，因此可控震源单炮信噪比是重点提高的方面。
提高可控震源信噪比的方法主要是通过可控震源台次、扫描频率、出力大小等试验，从中优选出最佳的施工因素，另外还有通过改变扫描信号，在激发阶段改变某一频率的输出，理论上增加某个频率的能量有2种方法，第1种是通过增加某个频率的扫描时间，通过时间积累来获得该频率的高能量，这就需要设计非线性扫描信号，而这种设计在增加某个频率的扫描时间的同时，又会降低其他频率的扫描时间。第2种是在扫描信号中直接提高某一频率的出力大小，扫描时间不变，但是基于可控震源机械工作的限制，可控震源在扫描过程中不可能改变出力大小，只能设置一个出力大小，例如70％，80％，那么这种信号无法应用于可控震源的工作中，但可以应用于相关过程。
因此，本领域亟需一种提高可控震源单炮信噪比的方法。
发明内容
本发明的目的是发明一种新的可控震源加权相关方法，提高可控震源单炮信噪比和记录面貌，推动可控震源技术进步。本发明将相关过程中的参考扫描信号进行改造，提升主频部分的能量大小，实践证明提升主频是非常有效的，通过将主频能量提升后获得加权扫描信号，与母记录进行相关获得单炮记录。
根据本发明的一个方面，提供一种提高信噪比的可控震源加权相关的方法，该方法包括以下步骤：
获得可控震源单炮的主频和主频部分的斜率；
改造参考信号的频谱，提升主频部分的能量，得到加权后的频谱曲线；
将该加权后的频谱曲线转化为加权扫描信号；
采用该加权扫描信号与原始记录进行相关，得到加权相关后的可控震源单炮记录。
进一步地，在可控震源施工得到的资料中，选择目的层部分开时窗进行频谱分析，获得可控震源单炮的主频和主频部分的斜率。
进一步地，以分贝的形式分析可控震源单炮的频谱。
进一步地，根据获得的主频部分的斜率进行计算主频部分的能量的提升量。
进一步地，提升主频部分的能量后，主频的斜率和单炮主频的斜率相当。
进一步地，将该加权后的频谱曲线转化为加权扫描信号的步骤包括：首先将扫描信号自相关得到扫描信号的子波，然后设计出一个雷克子波，采用前面分析得到的主频作为该雷克子波的频率，将该雷克子波乘以不同的系数，用以调节雷克子波振幅的大小；将不同振幅大小的雷克子波与扫描信号子波相加，相加后资料的频谱中要求主频突出，斜率和实际单炮主频的斜率相当。
进一步地，如果相加后资料频率不满足要求，继续调节雷克子波的系数。
进一步地，形式为一个子波，将该子波与扫描信号褶积，获得加权扫描信号。
本发明的方法通过在相关前提高参考信号的主频部分能量，转化为加权扫描信号，再将加权扫描信号与原始母记录进行相关得到单炮记录，使用该方法得到的单炮记录较常规方法得到的单炮记录信噪比更高。
附图说明
通过结合附图对本公开示例性实施方式进行更详细的描述，本公开的上述以及其它目的、特征和优势将变得更加明显，其中，在本公开示例性实施方式中，相同的参考标号通常代表相同部件。
图1为根据本发明实施例的方法流程图。
图2为根据本发明实施例的加权扫描信号图。
图3为应用本发明方法前后的效果图，3a为原始单炮效果图，3b为方法应用后单炮效果图。
图4为根据本发明实施例的确定主频率斜率示意图，4a为单炮目的层段开时窗图，4b为主频率斜率选取图。
图5为根据本发明应用实施例的扫描信号波形图。
图6为根据本发明应用实施例的相加子波的波形图。
附图仅用于示例性说明，不能理解为对本专利的限制；为了更好说明本实施例，附图某些部件会有省略、放大或缩小，并不代表实际产品的尺寸；对于本领域技术人员来说，附图中某些公知结构及其说明可能省略是可以理解的；相同或相似的标号对应相同或相似的部件；附图中描述位置关系的用语仅用于示例性说明，不能理解为对本专利的限制。
具体实施方式
下面将参照附图更详细地描述本公开的优选实施方式。虽然附图中显示了本公开的优选实施方式，然而应该理解，可以以各种形式实现本公开而不应被这里阐述的实施方式所限制。相反，提供这些实施方式是为了使本公开更加透彻和完整，并且能够将本公开的范围完整地传达给本领域的技术人员。
本发明是一种针对提升单炮信噪比的可控震源加权相关技术，该方法通过在相关前提高参考信号的主频部分能量，转化为加权扫描信号，再将加权扫描信号与原始母记录进行相关得到单炮记录，使用该方法得到的单炮记录较常规方法得到的单炮记录信噪比更高。
本发明是提高信噪比的可控震源加权相关方法，其方法包括首先需要选择需要加权相关的频率，通过频谱分析得到可控震源单炮的主频以及主频部分的斜率。例如，参照图4，可以根据实际单炮的频谱进行计算，在单炮目的层段开时窗(如图4a所示)，根据主频的幅度计算出其斜率(如图4b所示)，图中箭头就是斜率。
一般以分贝的形式分析单炮的频谱，主频就是需要加权相关的频率，将参考信号(线性扫描)的频谱进行改造，提高主频部分的能量，提高量的大小按单炮中主频部分的斜率进行计算。提高主频能量后，其主频的斜率和前面单炮主频的斜率相当。
得到加权后的扫描信号的频谱，将该加权后的频谱转化为加权扫描信号。作为一种实施方式，首先将扫描信号自相关得到扫描信号的子波，然后通过处理软件中的模块设计出一个雷克子波，采用前面分析得到的主频作为该雷克子波的频率，然后将该雷克子波乘以不同的系数(例如1.2、1.4、2.3….等)，用以调节雷克子波振幅的大小，将不同振幅大小的雷克子波与扫描信号子波相加，相加后资料的频谱中要求主频如前所述，主频突出，斜率和实际单炮主频的斜率相当，如果相加后资料频率不满足要求，可继续调节雷克子波的系数，直到满意为止。例如，获得的相加后资料，其形式仍然是一个子波，将该子波与扫描信号褶积，获得加权扫描信号，如图2所示。
采用该加权扫描信号与母记录进行相关得到单炮记录，该单炮记录比常规方法得到的单炮信噪比更高，资料面貌更好，可提升低信噪比地区的可控震源单炮效果。例如，可以在处理软件中，将该加权扫描信号与母记录进行互相关运算，处理软件中通常提供互相关运算模块。
为使本发明的目的、特征和优点能更明显易懂，下文配合所附图式，作详细说明如下：
如图1所示，图1为本发明的提高信噪比的可控震源加权相关方法的流程图。
在步骤101，在可控震源施工得到的资料中，选择目的层部分开时窗进行频谱分析，流程进入到步骤102。
在步骤102，通过频谱分析获得目的层段的主频和主频部分的斜率，流程进入到步骤103。
在步骤103，将参考信号(一般为线性扫描信号)的频谱进行改造，提升主频部分的能量，其提升量的大小根据前面分析的主频部分的斜率进行计算，得到加权后的频谱曲线，流程进入到步骤104。
在步骤104，将频谱曲线分段拟合成公式，并通过公式将该加权后的频谱转化为加权扫描信号(如图2)，流程进入到步骤105。
在步骤105，采用该加权扫描信号与原始记录进行相关。流程进入到步骤106。
在步骤106，得到加权相关后的可控震源单炮，该方法得到的单炮信噪比更高，单炮记录面貌更好(如图3)，流程结束。
为便于理解本发明实施例的方案及其效果，以下给出一个具体应用示例。本领域技术人员应理解，该示例仅为了便于理解本发明，其任何具体细节并非意在以任何方式限制本发明。
在西部某山前带工区，采用的是线性扫描信号6-70Hz，16s扫描，频谱及扫描信号波形如图5所示。根据实际单炮进行频谱分析，在主要目的层段开时窗，并计算得到主频及其斜率，如图4所示。
将扫描信号进行自相关得到扫描信号子波，并设计出主频为单炮主频的雷克子波，改变该雷克子波的振幅大小与扫描信号子波相加，得到的相加子波的频谱需要满足主频突出、主频斜率与单炮主频斜率相当，如图6所示，如不满足继续调整雷克子波的振幅。将最终获得的相加子波与扫描信号进行互相关，得到加权后的扫描信号，如图2所示，将加权后的扫描信号与母记录进行互相关运算，得到单炮记录。
显然，本发明的上述实施例仅仅是为清楚地说明本发明所作的举例，而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明权利要求的保护范围之内。
Claims (8)
1.一种提高信噪比的可控震源加权相关的方法，其特征在于，该方法包括以下步骤：
获得可控震源单炮的主频和主频部分的斜率；
改造参考信号的频谱，提升主频部分的能量，得到加权后的频谱曲线；
将该加权后的频谱曲线转化为加权扫描信号；
采用该加权扫描信号与原始记录进行相关，得到加权相关后的可控震源单炮记录；
其中，将该加权后的频谱曲线转化为加权扫描信号的步骤包括：首先将扫描信号自相关得到扫描信号的子波，然后设计出一个雷克子波，采用前面分析得到的主频作为该雷克子波的频率，将该雷克子波乘以不同的系数，用以调节雷克子波振幅的大小；
将不同振幅大小的雷克子波与扫描信号子波相加，相加后资料的频谱中要求主频突出，斜率和实际单炮主频的斜率相当。
2.根据权利要求1所述的方法，其特征在于，在可控震源施工得到的资料中，选择目的层部分开时窗进行频谱分析，获得可控震源单炮的主频和主频部分的斜率。
3.根据权利要求2所述的方法，其特征在于，以分贝的形式分析可控震源单炮的频谱。
4.根据权利要求1所述的方法，其特征在于，根据获得的主频部分的斜率进行计算主频部分的能量的提升量。
5.根据权利要求1所述的方法，其特征在于，参考信号为线性扫描信号。
6.根据权利要求1所述的方法，其特征在于，提升主频部分的能量后，主频的斜率和单炮主频的斜率相当。
7.根据权利要求1所述的方法，其特征在于，如果相加后资料频率不满足要求，继续调节雷克子波的系数。
8.根据权利要求1所述的方法，其特征在于，获得的相加后资料，形式为一个子波，将该子波与扫描信号褶积，获得加权扫描信号。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610309760.XA CN107367758B (zh) | 2016-05-11 | 2016-05-11 | 一种提高信噪比的可控震源加权相关的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610309760.XA CN107367758B (zh) | 2016-05-11 | 2016-05-11 | 一种提高信噪比的可控震源加权相关的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107367758A CN107367758A (zh) | 2017-11-21 |
| CN107367758B true CN107367758B (zh) | 2019-06-04 |
Family
ID=60304555
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610309760.XA Active CN107367758B (zh) | 2016-05-11 | 2016-05-11 | 一种提高信噪比的可控震源加权相关的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107367758B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109975874B (zh) * | 2017-12-28 | 2020-05-19 | 中石化石油工程技术服务有限公司 | 基于阻尼雷克子波的可控震源扫描信号设计方法 |
| CN110888165B (zh) * | 2018-09-10 | 2021-08-24 | 中国石油化工股份有限公司 | 一种可控震源次谐波噪声的压制方法及系统 |
| CN110954956B (zh) * | 2018-09-26 | 2021-10-22 | 中国石油化工股份有限公司 | 观测系统的采集痕迹评价方法及计算机可读存储介质 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104199093A (zh) * | 2014-09-01 | 2014-12-10 | 中国海洋石油总公司 | 基于时频域能量自适应加权的地震信号分辨率增强方法 |
| CN105301634A (zh) * | 2014-06-26 | 2016-02-03 | 中石化石油工程地球物理有限公司胜利分公司 | 利用近震源检波器记录提高可控震源单炮相关质量的方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8000168B2 (en) * | 2006-12-08 | 2011-08-16 | Conocophillips Company | Dynamic source parameter selection for seismic vibrator data acquisition |
-
2016
- 2016-05-11 CN CN201610309760.XA patent/CN107367758B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105301634A (zh) * | 2014-06-26 | 2016-02-03 | 中石化石油工程地球物理有限公司胜利分公司 | 利用近震源检波器记录提高可控震源单炮相关质量的方法 |
| CN104199093A (zh) * | 2014-09-01 | 2014-12-10 | 中国海洋石油总公司 | 基于时频域能量自适应加权的地震信号分辨率增强方法 |
Non-Patent Citations (1)
| Title |
|---|
| "可控震源高保真采集数值模拟及炮集分离";张洁 等;《石油物探》;20120331;第51卷(第2期);第178-183页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107367758A (zh) | 2017-11-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Bagaini | Low-frequency vibroseis data with maximum displacement sweeps | |
| Odaka et al. | A new method of quickly estimating epicentral distance and magnitude from a single seismic record | |
| CA1216353A (en) | Vibratory signal sweep seismic prospecting method and apparatus | |
| US7330401B2 (en) | Seismic vibratory acquisition method and apparatus | |
| CA2277119C (en) | Seismic data acquisition and processing using non-linear distortion in a groundforce signal | |
| O’Connell et al. | Interferometric multichannel analysis of surface waves (IMASW) | |
| CN109782335B (zh) | 低畸变宽频扫描信号设计方法 | |
| JP3617036B2 (ja) | 地盤探査装置及びそれに使用される解析プログラム | |
| CN107367758B (zh) | 一种提高信噪比的可控震源加权相关的方法 | |
| WO2013138956A1 (zh) | 一种确定可控震源最佳低频扫描信号的方法 | |
| US6665619B2 (en) | Noise estimator for seismic exploration | |
| CN105572723B (zh) | 基于自相关子波的可控震源扫描信号的设计方法 | |
| MX2013004281A (es) | Metodos y aparato para que por prospeccion geofisica detecte cuerpos de fluidos en formaciones subterraneas. | |
| KR100412097B1 (ko) | 탄성파 자료 시변 최적 옵셋 처리 방법, 시스템 및 장치 | |
| Atkinson et al. | Attenuation and source characteristics of the 23 June 2010 M 5.0 Val-des-Bois, Quebec, earthquake | |
| CN100501450C (zh) | 一种对地震信号进行波形拉伸补偿的方法 | |
| Satoh | Empirical frequency-dependent radiation pattern of the 1998 Miyagiken-Nanbu earthquake in Japan | |
| Pitarka et al. | Validation of a 3D velocity model of the Puget Sound region based on modeling ground motion from the 28 February 2001 Nisqually earthquake | |
| CN104793237A (zh) | 一种获得宽频可控震源扫描信号的方法和装置 | |
| Taipodia et al. | Impact of strike energy on the resolution of dispersion image in active MASW survey | |
| WO2003032008A2 (en) | System and method for assigning exploration risk to seismic attributes | |
| Riga et al. | Assessing the applicability of deconvolution of borehole records for determining near‐surface shear‐wave attenuation | |
| CN100412569C (zh) | 利用地震微测井进行地震信号高频补偿方法 | |
| CN109143345A (zh) | 基于模拟退火的品质因子q非线性反演方法及系统 | |
| CN110082820B (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 |