CN107340536A - 地震反演剥离煤层强反射的方法和装置 - Google Patents
地震反演剥离煤层强反射的方法和装置 Download PDFInfo
- Publication number
- CN107340536A CN107340536A CN201610286684.5A CN201610286684A CN107340536A CN 107340536 A CN107340536 A CN 107340536A CN 201610286684 A CN201610286684 A CN 201610286684A CN 107340536 A CN107340536 A CN 107340536A
- Authority
- CN
- China
- Prior art keywords
- reflection
- coal seam
- mrow
- seismic
- strong
- 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
- 239000003245 coal Substances 0.000 title claims abstract description 78
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000013508 migration Methods 0.000 claims abstract description 20
- 230000005012 migration Effects 0.000 claims abstract description 20
- 238000001228 spectrum Methods 0.000 claims abstract description 16
- 230000004044 response Effects 0.000 claims description 15
- 238000012937 correction Methods 0.000 claims description 12
- 239000011159 matrix material Substances 0.000 claims description 8
- 238000013075 data extraction Methods 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 6
- 238000004422 calculation algorithm Methods 0.000 claims description 5
- 238000002922 simulated annealing Methods 0.000 claims description 4
- 238000013517 stratification Methods 0.000 claims description 4
- 230000002194 synthesizing effect Effects 0.000 claims description 3
- 230000002123 temporal effect Effects 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 15
- 230000006870 function Effects 0.000 description 14
- 238000012545 processing Methods 0.000 description 9
- 238000011161 development Methods 0.000 description 8
- 238000004590 computer program Methods 0.000 description 5
- 230000005540 biological transmission Effects 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 230000002159 abnormal effect Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 239000000835 fiber Substances 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000001902 propagating effect Effects 0.000 description 2
- 239000004576 sand Substances 0.000 description 2
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 1
- 230000005856 abnormality Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 229910052802 copper Inorganic materials 0.000 description 1
- 239000010949 copper Substances 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 230000003068 static effect 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
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
本发明公开了地震反演剥离煤层强反射的方法和装置。该方法包括：提取原始地震偏移成果数据体S(t)中每一道的统计子波Wn(fm)以及部分频谱针对每一道，基于使目标函数Ja极小反演得到顶界面等效反射系数rn,l,u、底界面等效反射系数rn,l,d以及rn,l,u和rn,l,d分别对应的时间位置和基于rn,l,u、rn,l,d、得到煤层强反射数据体S*(t)；以及基于下式得到储层弱反射数据体Sw(t)。应用本发明，能够精确估算煤层的反射能量和反射位置，实现准确剥离煤层强反射、突出储层弱反射信息的目的。
Description
技术领域
本发明涉及地震资料处理领域，更具体地，涉及一种地震反演剥离煤层强反射的方法和一种地震反演剥离煤层强反射的装置。
背景技术
煤层具有低速度低密度特性，与其他介质形成强反射界面，将储层反射淹没其中，所以煤系地层地震资料一般缺失内幕反射，为储层预测带了很大的困难。去除煤层地震响应，突显储层弱反射信息是目前煤系储层预测常用手段。
Mcglen(2005)探讨了煤线对砂岩储层地震反射响应特征的严重干涉，并采用测井资料正演出煤层地震响应，进而从地震剖面上消除煤层响应完成储层预测；Raguwanti(2009)利用地震、测井和测试资料，采用地震波阻抗反演方法进行煤系地层的砂体预测，并进行井位筛选和地震、地质动态跟踪研究；赵爽(2007)、佘刚(2013)等人利用多子波分解与重构方法把单道地震数分解为对应的多个不同频率雷克子波，并对地震数据进行分解、处理、解释，并根据研究需要分解频率域内任何单频或多频地震剖面重构，实现排除煤层干扰、提高分辨率处理的目的。李海山等(匹配追踪煤层强反射分离方法，石油地球物理勘探2014,49(5)：866-870)借鉴匹配追踪算法得到信号在过完备字典上的有效稀疏表示，并找到最优稀疏表示字典，将地震信号中的煤层强反射信息匹配出来，进而消除煤层强反射对目的层有效信号的屏蔽。
还有人从原始地震数据上估算出煤层响应特征，随后从原始地震数据上减除，该方法对于煤层发育位置及煤层厚度缺乏考虑，其地震响应位置及能量估算的精确度不够，将此从地震数据上野蛮性减除，会造成储层弱信息突显过多或欠缺的问题，进而影响储层预测精度。
因此，现有技术在凸显储层弱反射信息方面都不能令人满意。
发明内容
本发明提出了一种用于凸显储层弱反射信息的方法，本发明还提出了相应的装置。
根据本发明的一方面，提出了一种地震反演剥离煤层强反射的方法，该方法包括：提取原始地震偏移成果数据体S(t)中每一道的统计子波Wn(fm)以及部分频谱其中，原始地震偏移成果数据体S(t)包括N个道，下标n表示道号，n＝1,2,...,N，共M个频点，fm＝(m-1)×Δf+f1，[f1,fM]是预设的频谱区间；针对每一道，基于使目标函数Ja极小反演得到顶界面等效反射系数rn,l,u、底界面等效反射系数rn,l,d以及rn,l,u和rn,l,d分别对应的时间位置和
Ja＝||Gα-D||2，
其中，下标l表示煤层序号，l＝1,...,L，共L个煤层，G为M×2L的矩阵，G的第m行第2l-1列的元素第m行第2l列的元素α＝[rn,1,u,rn,1,d,...,rn,L,u,rn,L,d]T，基于rn，l，u、rn,l,d、得到煤层强反射数据体S*(t)；基于下式得到储层弱反射数据体Sw(t)：
Sw(t)＝S(t)-S*(t)。
根据本发明的另一方面，提出了一种地震反演剥离煤层强反射的装置，该装置包括：数据提取模块，用于提取原始地震偏移成果数据体S(t)中每一道的统计子波Wn(fm)以及部分频谱其中，原始地震偏移成果数据体S(t)包括N个道，下标n表示道号，n＝1,2,...,N，共M个频点，fm＝(m-1)×Δf+f1， [f1,fM]是预设的频谱区间；等效反射系数求解模块，针对每一道，用于基于使目标函数Ja极小反演得到顶界面等效反射系数rn,l,u、底界面等效反射系数rn,l,d以及rn,l,u和rn,l,d分别对应的时间位置和
Ja＝||Gα-D||2，
其中，下标l表示煤层序号，l＝1,...,L，共L个煤层，G为M×2L的矩阵，G的第m行第2l-1列的元素第m行第2l列的元素α＝[rn,1,u,rn,1,d,...,rn,L,u,rn,L,d]T，强反射估算模块，用于基于rn,l,u、rn,l,d、得到煤层强反射数据体S*(t)；弱反射获取模块，用于基于下式得到储层弱反射数据体Sw(t)：
Sw(t)＝S(t)-S*(t)。
本发明通过各方面通过精确估算煤层的反射能量和反射位置，从而实现准确剥离煤层强反射、突出储层弱反射信息的目的。
附图说明
通过结合附图对本发明示例性实施方式进行更详细的描述，本发明的上述以及其它目的、特征和优势将变得更加明显，其中，在本发明示例性实施方式中，相同的参考标号通常代表相同部件。
图1示出了根据本发明的一个实施例的地震反演剥离煤层强反射的方法的流程图。
图2(a)～(c)示出了应用本发明的模型试算示意图。
图3(a)～(b)示出了应用本发明获取某实际研究区的储层弱反射信息的示意图。
具体实施方式
下面将参照附图更详细地描述本发明的优选实施方式。虽然附图中显示了本发明的优选实施方式，然而应该理解，可以以各种形式实现本发明而不应被这里阐述的实施方式所限制。相反，提供这些实施方式是为了使本发明更加透彻和完整，并且能够将本发明的范围完整地传达给本领域的技术人员。
此处，先对本发明的基本原理进行简单介绍。
我们将煤层顶界面(简称为顶界面)时间记作tn,l,u，底界面(简称为底界面)时间记为tn,l,d，此处下标n表示道号，用于区分不同的道(即地震道)，也可以采用线号+道号的方式来区分不同的道；下标l表示煤层序号，设共L个煤层，则l＝1,...,L。考虑到煤层较薄且在地震调谐厚度以下，可以将顶界面等效反射系数的记作rn,l,u，底界面的等效反射系数大小记作rn,l,d。可以根据褶积模型和随机非线性反演方法，在校正顶、底界面的时间位置的同时，反演得到顶、底界面等效反射系数rn,l,u、rn,l,d。再根据褶积公式合成煤层强反射数据体(也称为煤层强地震响应)S*(t)，将其从原始地震偏移成果数据体S(t)中减除，突显出储层若反射数据体(也称为储层弱反射信息)，从而为下一步储层预测提供更能反映储层岩特征性质的地震数据。
考虑地震数据的噪声影响，高频和低频数据不能作为真实的地震响应，因此，可考虑将地震偏移数据转换到频率域，提取瞬时频率数据体，得到其部分频谱观测数据D的表达式。为描述简要，以下描述以仅一个煤层为例进行原理说明。设地震偏移成果数据体S(t)包括N个道，针对每一道，有：
其中，下标n表示道号，n＝1,2,...,N；共M个频点，fm＝(m-1)×Δf+f1， [f1,fM]是预设的频谱区间；表示单道信号的部分频谱；Sn(fm)，m＝1,2,...,M表示单道信号；Wn(fm)，m＝1,2,...,M表示单道信号的统计子波；Rn(fm)，m＝1,2,...,M表示单道信号的频谱；Nn(fm)，m＝1,2,...,M表示单道信号的噪音。
公式(1)可以写成矩阵形式：
其中，G为M×2的矩阵，G的第m行第1列的元素第m行第2列的元素α＝[rn,u,rn,d]T，
公式(2)中，等效反射系数的时间位置和大小被分裂开，公式(2)是求取4个未知数rn,u、rn,d、的非线性方程组。
可以根据公式(2)建立数据匹配平方最小的目标函数：
Ja＝||Gα-D||2， (3)
当求取的rn,u、rn,d、使得Ja取得极小时，即为公式(3)的解。
需要注意的是，上述结果是在假设L＝1的情况下得到的。实际应用中，可能只有一个煤层(即L＝1)或多个煤层(L>1)。
实施例1
图1示出了根据本发明的一个实施例的地震反演剥离煤层强反射的方法的流程图，该方法包括：
步骤101，提取原始地震偏移成果数据体S(t)中每一道的统计子波Wn(fm)以及部分频谱其中，原始地震偏移成果数据体S(t)包括N个道，下标n表示道号，n＝1,2,...,N，共M个频点，fm＝(m-1)×Δf+f1，[f1,fM]是预设的频谱区间；
步骤102，针对每一道，基于使目标函数Ja极小反演得到顶界面等效反射系数rn,l,u、底界面等效反射系数rn,l,d以及rn,l,u和rn,l,d分别对应的时间位置和
Ja＝||Gα-D||2，
其中，下标l表示煤层序号，l＝1,...,L，共L个煤层，G为M×2L的矩阵，G的第m行第2l-1列的元素第m行第2l列的元素α＝[rn,1,u,rn,1,d,...,rn,L,u,rn,L,d]T，
步骤103，基于rn,l,u、rn,l,d、得到煤层强反射数据体S*(t)；以及
步骤104，基于下式得到储层弱反射数据体Sw(t)：
Sw(t)＝S(t)-S*(t)。
本实施例通过精确估算煤层的反射能量和反射位置，从而实现准确剥离煤层强反射、突出储层弱反射信息的目的。
在步骤101中，提取每一道的Wn(fm)和可以包括：可以结合地质的分层、测井曲线以及原始地震偏移成果数据体S(t)，进行井震标定，以提取每一道的Wn(fm)和
在步骤102中，可以采用非常快速模拟退火(VFSA)算法求解目标函数Ja以得到rn,l,u、rn,l,d、在求解过程中，可以设置初始的时间位置分别是校正前相应的顶、底界面时间；并且可以设置第1道的初始的顶、底界面等效反射系数均为0，以及因煤层发育横向连续性较好，可以设置除第1道外每一道的煤层l的初始的顶、底界面等效反射系数为上一道求解得到的rn-1,l,u和rn-1,l,d，以加快搜索速度。
进一步地，发明人经大量实验发现，可以令时间位置与校正前相应的顶界面时间的差值小于5ms、令时间位置与校正前相应的底界面时间的差值小于5ms，有利于进一步提高估算的准确度。
在步骤103中，得到煤层强反射数据体S*(t)可以包括：
针对每一道，可以基于下式得到该道的强地震响应
其中，符号“*”表示褶积运算，Wn(t)是Wn(fm)，m＝1,2,...,M对应的时域信号，Rn(t)是基于rn,l,u、rn,l,d、得到的该道的反射系数；
可以通过合成所述N个道的强地震响应得到煤层强反射数据体S*(t)。
实施例2
本发明还公开了一种地震反演剥离煤层强反射的装置。该装置包括数据提取模块、等效反射系数求解模块、强反射估算模块和弱反射获取模块。数据提取模块用于提取原始地震偏移成果数据体S(t)中每一道的统计子波Wn(fm)以及部分频谱其中，原始地震偏移成果数据体S(t)包括N个道，下标n表示道号，n＝1,2,...,N，共M个频点，fm＝(m-1)×Δf+f1，[f1,fM]是预设的频谱区间。针对每一道，等效反射系数求解模块用于基于使目标函数Ja极小反演得到顶界面等效反射系数rn,l,u、底界面等效反射系数rn,l,d以及rn,l,u和rn,l,d分别对应的时间位置和
Ja＝||Gα-D||2，
其中，下标l表示煤层序号，l＝1,...,L，共L个煤层，G为M×2L的矩阵，G的第m行第2l-1列的元素第m行第2l列的元素α＝[rn,1,u,rn,1,d,...,rn,L,u,rn,L,d]T，强反射估算模块用于基于rn,l,u、rn,l,d、得到煤层强反射数据体S*(t)。弱反射获取模块用于基于下式得到储层弱反射数据体Sw(t)：
Sw(t)＝S(t)-S*(t)。
在一种可能的实施方式中，在数据提取模块中，提取每一道的Wn(fm)和可以包括：可以结合地质的分层、测井曲线以及原始地震偏移成果数据体S(t)，进行井震标定，以提取每一道的Wn(fm)和
在一种可能的实施方式中，在等效反射系数求解模块中，可以采用非常快速模拟退火算法求解目标函数Ja以得到rn,l,u、rn,l,d、其中，在求解过程中，可以设置初始的时间位置分别是校正前相应的顶、底界面时间，并且可以设置第1道的初始的顶、底界面等效反射系数均为0，以及可以设置除第1道外每一道的煤层l的初始的顶、底界面等效反射系数为上一道求解得到的rn-1,l,u和rn-1,l,d。
进一步地，针对每一道，时间位置与校正前的顶界面时间的差值可以小于5ms，时间位置与校正前的底界面时间的差值可以小于5ms。
在强反射估算模块中，得到煤层强反射数据体S*(t)可以包括：
针对每一道，可以基于下式得到该道的强地震响应
其中，符号“*”表示褶积运算，Wn(t)是Wn(fm)，m＝1,2,...,M对应的时域信号，Rn(t)是基于rn,l,u、rn,l,d、得到的该道的反射系数；可以通过合成所述N个道的强地震响应得到煤层强反射数据体S*(t)。
应用示例
为便于理解本发明实施例的方案及其效果，以下给出一个具体应用示例。本领域技术人员应理解，该示例仅为了便于理解本发明，其任何具体细节并非意在以任何方式限制本发明。
图2(a)～(c)为本发明的二维模型试算图件。图2(a)是二维模型——砂泥岩地层地震记录，波谷下倾严重的地方为发育的砂岩层；图2(b)为在原来二维砂泥岩模型中，设计砂岩上覆一层煤层，其地震剖面上能量显著增强，在砂岩发育的地方的波谷位置为强能量，看不出异常；图2(c)为应用本发明去除煤层强反射数据体后的地震数据，可以看出砂岩发育的波谷位置异常明显，波谷异常发育位置与图2(a)一致，且能量相当，表明本发明在理论上具有较好的精度。
图3为本发明在实际数据运算图件。图3(a)为实际研究区内，煤系储层中沿目的层段提取的波峰振幅属性，受煤层影响，储层发育位置为白色显示的弱振峰能量(椭圆内)，而异常的黑色区域实际钻井揭示是煤层发育区域；图3(b)是应用本发明去除煤层强反射后的目的层段波峰振幅属性平面图，储层发育异常得以突见(椭圆内)，厚煤层发育区域能量得以消除。实际数据运算结果再次证明本发明在煤系储层预测中有很好的应用效果。
本发明可以是系统、方法和/或计算机程序产品。计算机程序产品可以包括计算机可读存储介质，其上载有用于使处理器实现本发明的各个方面的计算机可读程序指令。
计算机可读存储介质可以是可以保持和存储由指令执行设备使用的指令的有形设备。计算机可读存储介质例如可以是――但不限于――电存储设备、磁存储设备、光存储设备、电磁存储设备、半导体存储设备或者上述的任意合适的组合。计算机可读存储介质的更具体的例子(非穷举的列表)包括：便携式计算机盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、静态随机存取存储器(SRAM)、便携式压缩盘只读存储器(CD-ROM)、数字多功能盘(DVD)、记忆棒、软盘、机械编码设备、例如其上存储有指令的打孔卡或凹槽内凸起结构、以及上述的任意合适的组合。这里所使用的计算机可读存储介质不被解释为瞬时信号本身，诸如无线电波或者其他自由传播的电磁波、通过波导或其他传输媒介传播的电磁波(例如，通过光纤电缆的光脉冲)、或者通过电线传输的电信号。
这里所描述的计算机可读程序指令可以从计算机可读存储介质下载到各个计算/处理设备，或者通过网络、例如因特网、局域网、广域网和/或无线网下载到外部计算机或外部存储设备。网络可以包括铜传输电缆、光纤传输、无线传输、路由器、防火墙、交换机、网关计算机和/或边缘服务器。每个计算/处理设备中的网络适配卡或者网络接口从网络接收计算机可读程序指令，并转发该计算机可读程序指令，以供存储在各个计算/处理设备中的计算机可读存储介质中。
用于执行本发明操作的计算机程序指令可以是汇编指令、指令集架构(ISA)指令、机器指令、机器相关指令、微代码、固件指令、状态设置数据、或者以一种或多种编程语言的任意组合编写的源代码或目标代码，所述编程语言包括面向对象的编程语言—诸如Smalltalk、C++等，以及常规的过程式编程语言—诸如“C”语言或类似的编程语言。计算机可读程序指令可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行。在涉及远程计算机的情形中，远程计算机可以通过任意种类的网络—包括局域网(LAN)或广域网(WAN)—连接到用户计算机，或者，可以连接到外部计算机(例如利用因特网服务提供商来通过因特网连接)。在一些实施例中，通过利用计算机可读程序指令的状态信息来个性化定制电子电路，例如可编程逻辑电路、现场可编程门阵列(FPGA)或可编程逻辑阵列(PLA)，该电子电路可以执行计算机可读程序指令，从而实现本发明的各个方面。
这里参照根据本发明实施例的方法、装置(系统)和计算机程序产品的流程图和/或框图描述了本发明的各个方面。应当理解，流程图和/或框图的每个方框以及流程图和/或框图中各方框的组合，都可以由计算机可读程序指令实现。
这些计算机可读程序指令可以提供给通用计算机、专用计算机或其它可编程数据处理装置的处理器，从而生产出一种机器，使得这些指令在通过计算机或其它可编程数据处理装置的处理器执行时，产生了实现流程图和/或框图中的一个或多个方框中规定的功能/动作的装置。也可以把这些计算机可读程序指令存储在计算机可读存储介质中，这些指令使得计算机、可编程数据处理装置和/或其他设备以特定方式工作，从而，存储有指令的计算机可读介质则包括一个制造品，其包括实现流程图和/或框图中的一个或多个方框中规定的功能/动作的各个方面的指令。
也可以把计算机可读程序指令加载到计算机、其它可编程数据处理装置、或其它设备上，使得在计算机、其它可编程数据处理装置或其它设备上执行一系列操作步骤，以产生计算机实现的过程，从而使得在计算机、其它可编程数据处理装置、或其它设备上执行的指令实现流程图和/或框图中的一个或多个方框中规定的功能/动作。
附图中的流程图和框图显示了根据本发明的多个实施例的系统、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上，流程图或框图中的每个方框可以代表一个模块、程序段或指令的一部分，所述模块、程序段或指令的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。在有些作为替换的实现中，方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如，两个连续的方框实际上可以基本并行地执行，它们有时也可以按相反的顺序执行，这依所涉及的功能而定。也要注意的是，框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合，可以用执行规定的功能或动作的专用的基于硬件的系统来实现，或者可以用专用硬件与计算机指令的组合来实现。
以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择，旨在最好地解释各实施例的原理、实际应用或对市场中的技术的改进，或者使本技术领域的其它普通技术人员能理解本文披露的各实施例。
Claims (10)
1.一种地震反演剥离煤层强反射的方法，该方法包括：
提取原始地震偏移成果数据体S(t)中每一道的统计子波Wn(fm)以及部分频谱其中，原始地震偏移成果数据体S(t)包括N个道，下标n表示道号，n＝1,2,...,N，共M个频点，fm＝(m-1)×Δf+f1，m＝1,2,...,M，[f1,fM]是预设的频谱区间；
针对每一道，基于使目标函数Ja极小反演得到顶界面等效反射系数rn,l,u、底界面等效反射系数rn,l,d以及rn,l,u和rn,l,d分别对应的时间位置和
Ja＝||Gα-D||2，
其中，下标l表示煤层序号，l＝1,...,L，共L个煤层，G为M×2L的矩阵，G的第m行第2l-1列的元素第m行第2l列的元素α＝[rn,1,u,rn,1,d,...,rn,L,u,rn,L,d]T，
基于rn,l,u、rn,l,d、得到煤层强反射数据体S*(t)；
基于下式得到储层弱反射数据体Sw(t)：
Sw(t)＝S(t)-S*(t)。
2.根据权利要求1所述的方法，其中，提取每一道的Wn(fm)和包括：
结合地质的分层、测井曲线以及原始地震偏移成果数据体S(t)，进行井震标定，以提取每一道的Wn(fm)和
3.根据权利要求1所述的方法，其中，采用非常快速模拟退火算法求解目标函数Ja以得到rn,l,u、rn,l,d、其中，在求解过程中，设置初始的时间位置分别是校正前相应的顶、底界面时间，并且设置第1道的初始的顶、底界面等效反射系数均为0，以及除第1道外每一道的煤层l的初始的顶、底界面等效反射系数为上一道求解得到的rn-1,l,u和rn-1,l,d。
4.根据权利要求3所述的方法，其中，针对每一道，时间位置与校正前相应的顶界面时间的差值小于5ms，时间位置与校正前的相应的底界面时间的差值小于5ms。
5.根据权利要求1所述的方法，其中，得到煤层强反射数据体S*(t)包括：
针对每一道，基于下式得到该道的强地震响应
<mrow>
<msubsup>
<mi>S</mi>
<mi>n</mi>
<mo>*</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>W</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>*</mo>
<msub>
<mi>R</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
其中，符号“*”表示褶积运算，Wn(t)是Wn(fm)，m＝1,2,...,M对应的时域信号，Rn(t)是基于rn,l,u、rn,l,d、l＝1,...,L得到的该道的反射系；
通过合成所述N个道的强地震响应得到煤层强反射数据体S*(t)。
6.一种地震反演剥离煤层强反射的装置，该装置包括：
数据提取模块，用于提取原始地震偏移成果数据体S(t)中每一道的统计子波Wn(fm)以及部分频谱其中，原始地震偏移成果数据体S(t)包括N个道，下标n表示道号，n＝1,2,...,N，共M个频点，fm＝(m-1)×Δf+f1，m＝1,2,...,M，[f1,fM]是预设的频谱区间；
等效反射系数求解模块，针对每一道，用于基于使目标函数Ja极小反演得到顶界面等效反射系数rn,l,u、底界面等效反射系数rn,l,d以及rn,l,u和rn,l,d分别对应的时间位置和
Ja＝||Gα-D||2，
其中，下标l表示煤层序号，l＝1,...,L，共L个煤层，G为M×2L的矩阵，G的第m行第2l-1列的元素第m行第2l列的元素α＝[rn,1,u,rn,1,d,...,rn,L,u,rn,L,d]T，
强反射估算模块，用于基于rn,l,u、rn,l,d、得到煤层强反射数据体S*(t)；
弱反射获取模块，用于基于下式得到储层弱反射数据体Sw(t)：
Sw(t)＝S(t)-S*(t)。
7.根据权利要求6所述的装置，其中，在数据提取模块中，提取每一道的Wn(fm)和包括：
结合地质的分层、测井曲线以及原始地震偏移成果数据体S(t)，进行井震标定，以提取每一道的Wn(fm)和
8.根据权利要求6所述的装置，其中，在等效反射系数求解模块中，采用非常快速模拟退火算法求解目标函数Ja以得到rn,l,u、rn,l,d、其中，在求解过程中，设置初始的时间位置分别是校正前相应的顶、底界面时间，并且设置第1道的初始的顶、底界面等效反射系数均为0，以及除第1道外每一道的煤层l的初始的顶、底界面等效反射系数为上一道求解得到的rn-1,l,u和rn-1,l,d。
9.根据权利要求8所述的装置，其中，针对每一道，时间位置与校正前相应的顶界面时间的差值小于5ms，时间位置与校正前相应的底界面时间的差值小于5ms。
10.根据权利要求6所述的装置，其中，在强反射估算模块中，得到煤层强反射数据体S*(t)包括：
针对每一道，基于下式得到该道的强地震响应
<mrow>
<msubsup>
<mi>S</mi>
<mi>n</mi>
<mo>*</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>W</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>*</mo>
<msub>
<mi>R</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
其中，符号“*”表示褶积运算，Wn(t)是Wn(fm)，m＝1,2,...,M对应的时域信号，Rn(t)是基于rn,l,u、rn,l,d、l＝1,...,L得到的该道的反射系数；
通过合成所述N个道的强地震响应得到煤层强反射数据体S*(t)。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610286684.5A CN107340536B (zh) | 2016-05-03 | 2016-05-03 | 地震反演剥离煤层强反射的方法和装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610286684.5A CN107340536B (zh) | 2016-05-03 | 2016-05-03 | 地震反演剥离煤层强反射的方法和装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107340536A true CN107340536A (zh) | 2017-11-10 |
| CN107340536B CN107340536B (zh) | 2019-02-01 |
Family
ID=60222927
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610286684.5A Active CN107340536B (zh) | 2016-05-03 | 2016-05-03 | 地震反演剥离煤层强反射的方法和装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107340536B (zh) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109801543A (zh) * | 2018-12-27 | 2019-05-24 | 中国石油大学(北京) | 用于铺设滑脱层的组件和方法 |
| CN110794461A (zh) * | 2018-08-01 | 2020-02-14 | 中国石油化工股份有限公司 | 一种强地震反射界面下的小尺度缝洞信息凸显方法 |
| CN111025393A (zh) * | 2019-12-28 | 2020-04-17 | 中海石油(中国)有限公司上海分公司 | 针对含薄煤层地层的储层预测方法、装置、设备及介质 |
| CN111505708A (zh) * | 2020-04-28 | 2020-08-07 | 西安交通大学 | 一种基于深度学习的强反射层剥离方法 |
| CN111610562A (zh) * | 2019-02-22 | 2020-09-01 | 中国石油天然气集团有限公司 | 地震数据中强反射信号分离方法及装置 |
| CN112305584A (zh) * | 2019-07-24 | 2021-02-02 | 中国石油化工股份有限公司 | 一种基于相移子波的薄煤层强反射分离方法 |
| CN113126155A (zh) * | 2021-04-01 | 2021-07-16 | 中国石油化工股份有限公司 | 一种针对分布于煤岩间受强反射影响的砂岩储层预测方法 |
| CN114428361A (zh) * | 2020-09-21 | 2022-05-03 | 中国石油化工股份有限公司 | 基于低频模型优化的表层岩溶缝洞体预测方法及电子设备 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE102010007655A1 (de) * | 2010-02-11 | 2011-08-11 | Marmot Passive Monitoring Technologies SA, Vaud | Verfahren und Vorrichtung zum Überwachen von natürlichem CO2-Lagerstätten |
| CN105093315A (zh) * | 2014-04-25 | 2015-11-25 | 中国石油化工股份有限公司 | 一种去除煤层强反射信号的方法 |
| CN105223608A (zh) * | 2015-08-07 | 2016-01-06 | 中国石油大学(华东) | 一种含煤强屏蔽缝洞型储层的地震预测与描述方法 |
| US9268047B2 (en) * | 2012-10-12 | 2016-02-23 | Rock Solid Images, Inc | Geophysical surveying |
| CN105353411A (zh) * | 2015-10-27 | 2016-02-24 | 中国石油大学(华东) | 超深层煤上储层砂体刻画与描述方法 |
-
2016
- 2016-05-03 CN CN201610286684.5A patent/CN107340536B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE102010007655A1 (de) * | 2010-02-11 | 2011-08-11 | Marmot Passive Monitoring Technologies SA, Vaud | Verfahren und Vorrichtung zum Überwachen von natürlichem CO2-Lagerstätten |
| US9268047B2 (en) * | 2012-10-12 | 2016-02-23 | Rock Solid Images, Inc | Geophysical surveying |
| CN105093315A (zh) * | 2014-04-25 | 2015-11-25 | 中国石油化工股份有限公司 | 一种去除煤层强反射信号的方法 |
| CN105223608A (zh) * | 2015-08-07 | 2016-01-06 | 中国石油大学(华东) | 一种含煤强屏蔽缝洞型储层的地震预测与描述方法 |
| CN105353411A (zh) * | 2015-10-27 | 2016-02-24 | 中国石油大学(华东) | 超深层煤上储层砂体刻画与描述方法 |
Non-Patent Citations (5)
| Title |
|---|
| YANGHUA WANG: "Multichannel matching pursuit for seismic trace decomposition", 《GEOPHYSICS》 * |
| YANGHUA WANG: "Seismic time-frequency spectral decomposition by matching pursuit", 《GEOPHYSICS》 * |
| 张在金等: "煤系地层地震强反射剥离方法研究及低频伴影分析", 《石油地球物理勘探》 * |
| 朱博华等: "匹配追踪强反射层分离方法及应用", 《石油物探》 * |
| 赵爽等: "多子波分解技术检测含煤砂岩储层", 《天然气工业》 * |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110794461B (zh) * | 2018-08-01 | 2021-09-24 | 中国石油化工股份有限公司 | 一种强地震反射界面下的小尺度缝洞信息凸显方法 |
| CN110794461A (zh) * | 2018-08-01 | 2020-02-14 | 中国石油化工股份有限公司 | 一种强地震反射界面下的小尺度缝洞信息凸显方法 |
| CN109801543A (zh) * | 2018-12-27 | 2019-05-24 | 中国石油大学(北京) | 用于铺设滑脱层的组件和方法 |
| CN111610562B (zh) * | 2019-02-22 | 2023-10-27 | 中国石油天然气集团有限公司 | 地震数据中强反射信号分离方法及装置 |
| CN111610562A (zh) * | 2019-02-22 | 2020-09-01 | 中国石油天然气集团有限公司 | 地震数据中强反射信号分离方法及装置 |
| CN112305584A (zh) * | 2019-07-24 | 2021-02-02 | 中国石油化工股份有限公司 | 一种基于相移子波的薄煤层强反射分离方法 |
| CN111025393A (zh) * | 2019-12-28 | 2020-04-17 | 中海石油(中国)有限公司上海分公司 | 针对含薄煤层地层的储层预测方法、装置、设备及介质 |
| CN111505708B (zh) * | 2020-04-28 | 2021-04-20 | 西安交通大学 | 一种基于深度学习的强反射层剥离方法 |
| US11243320B2 (en) | 2020-04-28 | 2022-02-08 | Xi'an Jiaotong University | Method of stripping strong reflection layer based on deep learning |
| CN111505708A (zh) * | 2020-04-28 | 2020-08-07 | 西安交通大学 | 一种基于深度学习的强反射层剥离方法 |
| CN114428361A (zh) * | 2020-09-21 | 2022-05-03 | 中国石油化工股份有限公司 | 基于低频模型优化的表层岩溶缝洞体预测方法及电子设备 |
| CN113126155A (zh) * | 2021-04-01 | 2021-07-16 | 中国石油化工股份有限公司 | 一种针对分布于煤岩间受强反射影响的砂岩储层预测方法 |
| CN113126155B (zh) * | 2021-04-01 | 2024-03-01 | 中国石油化工股份有限公司 | 一种针对分布于煤岩间受强反射影响的砂岩储层预测方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107340536B (zh) | 2019-02-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107340536B (zh) | 地震反演剥离煤层强反射的方法和装置 | |
| Gómez et al. | A simple method inspired by empirical mode decomposition for denoising seismic data | |
| Höcker et al. | Fast structural interpretation with structure-oriented filtering | |
| Li et al. | A method for low-frequency noise suppression based on mathematical morphology in microseismic monitoring | |
| Puryear et al. | Layer-thickness determination and stratigraphic interpretation using spectral inversion: Theory and application | |
| Bekara et al. | Random and coherent noise attenuation by empirical mode decomposition | |
| Alaei et al. | Seismic resolution enhancement using scale transform in the time-frequency domain | |
| Ghaderpour et al. | Antileakage least-squares spectral analysis for seismic data regularization and random noise attenuation | |
| EP3259620B1 (en) | Multistage full wavefield inversion process that generates a multiple free data set | |
| Wang et al. | Enhancing resolution of nonstationary seismic data by molecular-Gabor transform | |
| Chen | Automatic velocity analysis using high-resolution hyperbolic Radon transform | |
| Li et al. | Adaptive multiple subtraction based on 3D blind separation of convolved mixtures | |
| CN111060961B (zh) | 基于多信息约束反演的品质因子确定方法、装置及系统 | |
| Yue et al. | Suppression of periodic interference during tunnel seismic predictions via the Hankel-SVD-ICA method | |
| EP3217354B1 (en) | Seismic modeling system providing seismic survey data frequency domain inpainting and related methods | |
| DK3211594T3 (en) | SEISMIC MODELING SYSTEM, PROVIDING INPAINTING SEISMIC SURVEY DATA BASED ON COMPARISON BORDERS OF THE REGION UNDER MISCELLANEOUS AND RELATED PROCEDURES | |
| Naghizadeh et al. | Seismic data reconstruction using multidimensional prediction filters | |
| CN108508481B (zh) | 一种纵波转换波地震数据时间匹配的方法、装置及系统 | |
| US20170248713A1 (en) | Seismic modeling system providing seismic survey data spatial domain exemplar inpainting and related methods | |
| Li et al. | Removing abnormal environmental noise in nodal land seismic data using deep learning | |
| CN112213782B (zh) | 分相位地震数据的处理方法、装置和服务器 | |
| US20190369277A1 (en) | System for Generating a Model of an Underground Formation from Seismic Data | |
| CN106547021B (zh) | 基于单井褶积算法建立初始模型的方法和装置 | |
| CN114779332A (zh) | 地震数据沉积背景去除方法、装置及电子设备 | |
| Nguyen et al. | Seismic noise attenuation using curvelet transform and dip map data structure |
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 |