CN112859168A - 一种增强厚煤层下赋煤层成像质量的反褶积应用方法 - Google Patents
一种增强厚煤层下赋煤层成像质量的反褶积应用方法 Download PDFInfo
- Publication number
- CN112859168A CN112859168A CN202110141821.7A CN202110141821A CN112859168A CN 112859168 A CN112859168 A CN 112859168A CN 202110141821 A CN202110141821 A CN 202110141821A CN 112859168 A CN112859168 A CN 112859168A
- Authority
- CN
- China
- Prior art keywords
- deconvolution
- coal seam
- frequency
- convolution
- application method
- 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
- 239000003245 coal Substances 0.000 title claims abstract description 70
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000003384 imaging method Methods 0.000 title claims abstract description 19
- 230000002708 enhancing effect Effects 0.000 title claims abstract description 7
- 238000012937 correction Methods 0.000 claims abstract description 19
- 238000012360 testing method Methods 0.000 claims description 21
- 230000000694 effects Effects 0.000 claims description 17
- 230000003068 static effect Effects 0.000 claims description 10
- 238000010183 spectrum analysis Methods 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000012544 monitoring process Methods 0.000 claims description 3
- 238000001228 spectrum Methods 0.000 claims description 3
- 230000002787 reinforcement Effects 0.000 claims description 2
- 238000012545 processing Methods 0.000 abstract description 9
- 238000010276 construction Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000006835 compression Effects 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 208000024891 symptom Diseases 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
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/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/40—Transforming data representation
- G01V2210/48—Other transforms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/53—Statics correction, e.g. weathering layer or transformation to a datum
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
本发明提供一种增强厚煤层下赋煤层成像质量的反褶积应用方法，其特征在于，包括以下步骤：反褶积之前将地震资料的低频信号全部保留；在应用第一种反褶积方法的基础上，对CRP数据动校正拉平后再次串联使用另一种反褶积方法。本发明以改善处理流程中反褶积应用方式来提高厚煤层下赋可采煤层成像质量，进而为矿方提供下赋可采煤层的赋存信息，省去了为得到下赋煤层信息而进行二次地震勘探施工的成本。
Description
技术领域
本发明涉及本发明涉及一种地震数据处理方法，特别涉及一种反褶积应用方法，适用于地质任务要求提高厚煤层下赋煤层成像质量的地震资料处理流程。
背景技术
目前资料处理技术对于厚度较大的煤层下面赋存的其它可采煤层的成像质量很难提高，距离越近成像质量越差。
问题和缺点：下部可采煤层受上赋厚煤层的屏蔽影响，反射波能量弱，成像差，很难进行资料解释，进而难以给矿方提供下赋可采煤层的可靠技术支持。
发明内容
为了有效解决厚煤层下赋可采煤层的成像质量差，难以解释该煤层赋存形态的问题，本发明提供一种增强厚煤层下赋煤层成像质量的反褶积应用方法。
为解决上述问题，本发明所采用的技术方案如下：
一种增强厚煤层下赋煤层成像质量的反褶积应用方法，其特征在于，包括以下步骤：
反褶积之前将地震资料的低频信号全部保留；
在应用一种反褶积方法的基础上，对CRP数据动校正拉平后再次串联使用另一种反褶积方法。
所述低频信号为频率为0-30Hz的低频信号。
以上方法具体为：
地表一致性反褶积：对褶积步长值进行测试，以等间隔测试多个步长值，以单炮整体及目的层频谱分析为主要监控方案，设定每个步长值对应的主频拓宽值，比较获得最佳步长值，以及对应的最佳主频拓宽值；
预测反褶积：进行单道预测与多道预测的充分对比，褶积时窗选取针对厚煤层下赋煤层为计算时窗；并且对褶积因子进行充分测试，褶积因子以等间隔选取多个进行测试，对比测试效果，最终确定效果最佳的褶积因子；
两次反褶积串联使用后，再做一次剩余静校正，校正时窗同第二次反褶积时窗。
与现有技术相比，本发明的有益效果是：
本发明以改善处理流程中反褶积应用方式来提高厚煤层下赋可采煤层成像质量，进而为矿方提供下赋可采煤层的赋存信息，省去了为得到下赋煤层信息而进行二次地震勘探施工的成本。
附图说明
图1-5为本发明实施例中山西某区15号煤层的预测反褶积中，对褶积因子进行测试时，褶积因子以20等间隔测试，测试了20、40、60、80、100，测试效果如图1-5所示；
图6是本发明实施例中山西某区15号煤层的两次反褶积后频谱分析频带逐步拓宽效果逐步优化示意图；
图7是本发明实施例中山西某区15号煤层的两次反褶积后子波压缩效果逐步优化的效果示意图；
图8是本发明实施例的山西某区提高厚煤层下赋15号煤层的成像质量方法流程图；
图9为本发明实施例中山西某区15号煤层的两次反褶积串联使用后，再做一次剩余静校正后的厚煤层下赋煤层成像改善顺序示意图；
图10为淮北某区针对厚煤层7号煤下部50m的8煤层做两次反褶积处理后的效果图。
具体实施方式
为了使本发明的技术手段、特征、成效易于了解，下面结合具体实施方式，进一步阐述本发明。
以山西某区为试验区，提供了如下的技术方案：
一种增强厚煤层下赋煤层成像质量的反褶积应用方法，包括：
反褶积之前将地震资料的30Hz以内的低频信号全部保留；
在应用一种反褶积方法的基础上，对CRP数据动校正拉平后再次串联使用另一种反褶积方法。
本方法要求保留低频信息，即在设置滤波档时不滤除低频信息，保留0-30Hz的低频信号。
具体实验为：
本区上赋可采厚煤层主要有3号煤和9号煤，厚度均达10m以上。本次地质任务主要目的层为15号煤，赋存于9号煤下部80m左右。由于上部3号及9号煤的厚度较大，对15号煤屏蔽作用强，15号煤无法形成连续性较好的同相轴，进而无法解释其赋存状态。以此为出发点，进行了本次方法试验。
试验方法的思路为：在原始数据充分保留低频信号的基础上，叠前噪音衰减后做地表一致性反褶积处理，然后对数据做CRP方向的动校正拉平，进而优化速度，做好剩余静校正与速度的迭代，在剩余静校正量充分优化的基础上，串联预测反褶积，此方法很大程度提高了厚煤层下赋煤层的成像质量(见图8方法应用流程图)。具体措施如下：
第一次反褶积即所述的地表一致性反褶积：主要对褶积步长值进行测试，步长值以等间隔4测试4、8、12、16、20，以单炮整体及目的层频谱分析为主要监控方案。测试结果为：步长4主频拓展拓宽到15-55Hz，步长值8主频拓宽至12至65Hz，步长值12主频拓宽至8至80Hz，步长值16及20拓宽效果同步长值12。经比较，本区测试步长值12最佳，对应的主频从25-45Hz拓宽至8-80Hz，见图6中a到b。需要说明的是，步长值与主频拓宽范围的变化因原始资料而异，只能通过对原始资料测试不同步长值，然后对单炮进行频谱分析来监控主频拓宽范围，主频拓宽范围最大者最佳。
由于原始数据受地表施工条件、环境、天气、激发因素等诸多随机因素影响，所产生的噪音也因资料而异，因此噪音衰减方法必须对症下药，不同噪音有不同的针对性衰减方法。煤田勘探中常见的噪音有以下几类：低频面波、声波、线性干扰、工业干扰、野值干扰、多次波干扰、地滚波干扰等等。现在市面上应用成熟的各类处理软件对以上噪音均有较好的处理方法。去噪的最优原则为在去除干扰噪音的同时，保留有效反射信号，即对有效反射信号的伤害最低。本实施方式主要针对反褶积的应用方法，与噪音衰减同属地震资料处理流程中的两大模块，使用顺序为反褶积在噪音衰减后。
第二次串联的反褶积即预测反褶积：进行单道预测与多道预测的充分对比，褶积时窗选取针对厚煤层下赋煤层为计算时窗，本区选取厚煤层9号煤层下10ms为起始计算时窗；并且对褶积因子进行充分测试，褶积因子以等间隔20测试，本区测试了20、40、60、80、100，测试效果如图1-5所示。以目的层频谱低频端与高频端振幅均得到较好的加强为比较标准来对比测试效果，最终确定多道预测反褶积60的褶积因子效果最佳。频谱分析第二次反褶积后目的层频谱低频端与高频端振幅均得到较好的加强，见图6中b到c。
经两次反褶积的优化组合使用，子波压缩效果逐步优化，见图7。
两次反褶积串联使用后，再做一次剩余静校正，校正时窗同第二次反褶积时窗(厚煤层下10ms起始)，经本方法处理后，最终使得厚煤层下赋煤层成像质量大大改善，见图9，为该煤层的赋存形态解释提供较好基础。所述剩余静校正为：计算时窗上起厚煤层下10ms，下至目的层下200ms左右；计算频率保低频，设置范围0Hz至100Hz；计算倾角选取参考处理区块煤层起伏情况。
本方法应用中关键点有：
单炮噪音衰减后进行第一次反褶机处理，要保留30Hz以内的全部低频成分；
第二次串联反褶积之前需要进行剩余静校正与速度分析的逐步迭代优化，迭代次数大于等于三次；组合使用反褶积与速度分析、剩余静校正迭代在处理流程如图8所示。
第二次串联反褶积后再做一次剩余静校正，校正时窗同第二次反褶积计算时窗，即厚煤层下10ms为起始，进一步提高资料成像质量。
本方法在其他地区具有一般适用性，在淮北某区针对厚煤层7号煤下部50m(距离更近)的8煤层做针对性处理，亦收到较好的应用效果，见图10。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明保护的范围之内。
Claims (5)
1.一种增强厚煤层下赋煤层成像质量的反褶积应用方法，其特征在于，包括以下步骤：
反褶积之前将地震资料的30Hz以内低频信号全部保留；
在应用第一种反褶积方法的基础上，对CRP数据动校正拉平后再次串联使用另一种反褶积方法。
2.根据权利要求1所述的反褶积应用方法，其特征在于，具体包括：
所述第一种反褶积方法为地表一致性反褶积：对褶积步长值进行测试，以等间隔测试多个步长值，以单炮整体及目的层频谱分析为主要监控方案，设定每个步长值对应的主频拓宽值，比较获得最佳步长值，以及对应的最佳主频拓宽值；
所述另一种反褶积方法为预测反褶积：进行单道预测与多道预测的充分对比，褶积时窗选取针对厚煤层下赋煤层为计算时窗；并且对褶积因子进行充分测试，褶积因子以等间隔选取多个进行测试，对比测试效果，最终确定效果最佳的褶积因子；
两次反褶积串联使用后，再做一次剩余静校正，校正时窗同第二次反褶积时窗。
3.根据权利要求1所述的反褶积应用方法，其特征在于，所述低频信号为频率为0-30Hz的低频信号。
4.根据权利要求2所述的反褶积应用方法，其特征在于，所述步长值与主频拓宽范围的变化因原始资料而异，通过对原始资料测试不同步长值，然后对单炮进行频谱分析来监控主频拓宽范围，主频拓宽范围最大者对应的步长值为最佳步长值。
5.根据权利要求2所述的反褶积应用方法，其特征在于，所述确定效果最佳的褶积因子，其标准是：以目的层频谱低频端与高频端振幅均得到较好的加强为比较标准来对比测试效果，最终确定多道预测反褶积60的褶积因子效果最佳。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110141821.7A CN112859168A (zh) | 2021-02-02 | 2021-02-02 | 一种增强厚煤层下赋煤层成像质量的反褶积应用方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110141821.7A CN112859168A (zh) | 2021-02-02 | 2021-02-02 | 一种增强厚煤层下赋煤层成像质量的反褶积应用方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112859168A true CN112859168A (zh) | 2021-05-28 |
Family
ID=75986163
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110141821.7A Pending CN112859168A (zh) | 2021-02-02 | 2021-02-02 | 一种增强厚煤层下赋煤层成像质量的反褶积应用方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112859168A (zh) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102323617A (zh) * | 2011-06-13 | 2012-01-18 | 中国石油化工股份有限公司 | 一种复杂地表的二维地震资料连片处理方法 |
| CN102854526A (zh) * | 2011-06-30 | 2013-01-02 | 中国石油化工股份有限公司 | 一种多分量地震资料的处理方法 |
| WO2016063125A1 (en) * | 2014-10-23 | 2016-04-28 | Cgg Services Sa | Imaging the near subsurface with surface consistent deconvolution operators |
| CN106371140A (zh) * | 2016-08-17 | 2017-02-01 | 中国石油化工股份有限公司 | 一种提高中深层地震资料分辨率的方法 |
| CN112198547A (zh) * | 2019-07-08 | 2021-01-08 | 中国石油天然气集团有限公司 | 深层或超深层地震资料处理方法及装置 |
-
2021
- 2021-02-02 CN CN202110141821.7A patent/CN112859168A/zh active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102323617A (zh) * | 2011-06-13 | 2012-01-18 | 中国石油化工股份有限公司 | 一种复杂地表的二维地震资料连片处理方法 |
| CN102854526A (zh) * | 2011-06-30 | 2013-01-02 | 中国石油化工股份有限公司 | 一种多分量地震资料的处理方法 |
| WO2016063125A1 (en) * | 2014-10-23 | 2016-04-28 | Cgg Services Sa | Imaging the near subsurface with surface consistent deconvolution operators |
| CN106371140A (zh) * | 2016-08-17 | 2017-02-01 | 中国石油化工股份有限公司 | 一种提高中深层地震资料分辨率的方法 |
| CN112198547A (zh) * | 2019-07-08 | 2021-01-08 | 中国石油天然气集团有限公司 | 深层或超深层地震资料处理方法及装置 |
Non-Patent Citations (3)
| Title |
|---|
| 张锐: "松辽盆地太30地区保幅高分辨率处理技术研究", 《中国优秀博硕士学位论文全文数据库(硕士)基础科学辑》 * |
| 王秀荣等: "三维地震资料二次精细处理技术在新元煤矿的应用", 《中国资源综合利用》 * |
| 闫艳丽: "提高分辨率处理技术在胜利油田东部探区的应用", 《石油石化物资采购》 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109425896B (zh) | 白云岩油气储层分布预测方法及装置 | |
| CN109061764B (zh) | 一种分频融合波阻抗反演方法 | |
| CN108614295B (zh) | 一种基于广义地震子波的地层q值计算方法 | |
| CN109143374B (zh) | 一种井周散射体成像方法及系统 | |
| CN107065013B (zh) | 一种地震尺度下的层速度确定方法及装置 | |
| CN109633752B (zh) | 基于三维快速Radon变换的海上拖缆资料自适应鬼波压制方法 | |
| CN106547020B (zh) | 一种地震数据的保幅处理方法 | |
| CN104793237A (zh) | 一种获得宽频可控震源扫描信号的方法和装置 | |
| CN103913770A (zh) | 基于vsp资料对地震数据进行处理的方法 | |
| CN107238864A (zh) | 基于优势频带约束的地震波能量时频补偿方法 | |
| CN116125530A (zh) | 一种黄土塬区地震资料处理方法、系统及电子设备 | |
| CN112859168A (zh) | 一种增强厚煤层下赋煤层成像质量的反褶积应用方法 | |
| CN114428343A (zh) | 基于归一化互相关的Marchenko成像方法及系统 | |
| CN104502977A (zh) | 一种井控保幅高分辨率地震资料处理方法 | |
| CN107831533A (zh) | 基于能量统计的采集方案检波器排列长度定量分析方法 | |
| CN112198547A (zh) | 深层或超深层地震资料处理方法及装置 | |
| CN109975874B (zh) | 基于阻尼雷克子波的可控震源扫描信号设计方法 | |
| CN113514889B (zh) | 一种提升海洋深反射地震数据中低频信号能量的处理方法 | |
| CN111929726B (zh) | 地震相干数据体处理方法及装置 | |
| CN109884705B (zh) | 双重约束时频域子波谱提高地震分辨率处理方法 | |
| CN113671565A (zh) | 一种针对巨厚储层气藏开发的地震多尺度储层预测方法 | |
| CN113156512B (zh) | 适合可控震源激发的三维地震观测方法、设备及系统 | |
| CN113917539B (zh) | 火山岩覆盖区地震数据叠前道集处理方法、系统和装置 | |
| CN113064205B (zh) | 菲涅尔带约束的浅水多次波衰减方法 | |
| CN112946750B (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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210528 |