CN102928872B - 一种提高滩坝砂储层地震识别性能与描述精度的方法 - Google Patents
一种提高滩坝砂储层地震识别性能与描述精度的方法 Download PDFInfo
- Publication number
- CN102928872B CN102928872B CN201210405357.9A CN201210405357A CN102928872B CN 102928872 B CN102928872 B CN 102928872B CN 201210405357 A CN201210405357 A CN 201210405357A CN 102928872 B CN102928872 B CN 102928872B
- Authority
- CN
- China
- Prior art keywords
- wavelet
- sigma
- tau
- lambda
- seismic
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 239000004576 sand Substances 0.000 title claims abstract description 20
- 238000000605 extraction Methods 0.000 claims description 7
- 238000012545 processing Methods 0.000 claims description 7
- 230000004927 fusion Effects 0.000 claims description 5
- 230000001427 coherent effect Effects 0.000 claims description 3
- 239000003086 colorant Substances 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 abstract description 8
- 238000004458 analytical method Methods 0.000 abstract description 6
- 230000000694 effects Effects 0.000 description 7
- 238000006243 chemical reaction Methods 0.000 description 5
- 239000000284 extract Substances 0.000 description 4
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种提高滩坝砂储层地震识别性能与描述精度的方法,属于地震勘探领域。其是先采用三参数小波来提高地震分辨薄互层的能力,再通过多子波技术对地震数据进行变换,得到成果数据,后期利用沿层切片或者属性分析等技术,以提高滩坝砂储层的地震识别性能与描述精度。本发明能有效提高薄层识别能力,提高主频,展宽频带范围。经多子波变化后,能更加突出滩坝砂储层信息,经属性融合显示后,能大大提高后期储层描述的精度。
Description
技术领域
本发明涉及一种提高滩坝砂储层地震识别性能与描述精度的方法,属于地震勘探领域。
背景技术
地震勘探已经从构造油气藏向岩性油气藏过渡,大部分滩坝砂油气藏属于岩性油气藏,由于该类油藏埋藏深、单层薄、横向变化大、地震分辨率低等难点,制约了滩坝砂油藏的勘探开发,油田增储上产难度加大。
发明内容
本发明的任务在于提供一种提高滩坝砂储层地震识别性能与描述精度的方法。其是先采用三参数小波来提高地震分辨薄互层的能力,再通过多子波技术对地震数据进行变换,得到成果数据,后期利用沿层切片或者属性分析等技术,以提高滩坝砂储层的地震识别性能与描述精度。
其技术解决方案是:
一种提高滩坝砂储层地震识别性能与描述精度的方法,包括以下步骤:
a、首先利用三参数小波对地震资料进行提高分辨率处理;三参数小波公式为
其中,t表示时间,i表示虚数,σ为小波函数的调制频率,τ为能量衰减因子,β为能量延迟因子,(σ,τ,β∈R,σ,τ≥0),Λ=(σ,τ,β)记参数σ,τ,β集合;上式中的p(Λ),q(Λ)和k(Λ)分别为
b、根据选取的三个小波参数Λ=(σ,τ,β),得到三参数小波函数φ(t,Λ);设地震数据为X(t,j),j=1,2,...,M,其中M为道数、j表示第j道,按道进行褶积处理得到
s(t,j)=X(t,j)*φ(t,Λ),j=1,2,...,M
得到三参数小波处理后的数据;
c、对处理后的地震数据信号,首先估算滩坝砂储层位置的子波的频带范围[fmin,fmax],则可求出中心频率fc和频带宽度fw,设子波在该频带内的能量占子波总能量的λ2,则
将该式转化为求解以下方程
Aw=λw
其中,A={am,n}N×N,则得到
其中,N为子波样点数,Δt为采样间隔;取λ接近于1所对应的特征向量w,记为{w},{w}的频域能量限制在频带[fmin,fmax]以内,{w}为所要求取的多子波;
d、求取了标准子波以后,利用所求得的子波与地震道做相关处理
其中s为地震道,w为所求子波,C(t,j)便是相关信号,即为最后多子波变换数据;
e、对以上数据进行属性提取,对提取的属性进行多属性融合。
上述步骤a中,Λ=(σ,τ,β)可以根据实际情况进行选择,我们这里选择Λ=(4,0.5,0)。
上述步骤c中,求取的多子波一般是2~3个,后期处理时,应根据实际情况,选取合适的多子波进行分析。
上述步骤e中,所述多属性融合为二属性融合;即把二维颜色面板横坐标方向看作是一个属性,纵坐标方向是另外一个属性,这样通过两个属性的值的大小来选择二维面板里的颜色,多属性RGB显示技术本质上是一种图像处理技术,通过映射函数
IRGB(X)=S[IR(X),IG(X),IB(X)]
X=(x,y,z)代表在数据体中的位置,IR(X),IG(X),IB(X)分别代表红、绿、蓝三种颜色上归一化的属性分量,通过三维颜色体,将三个属性分量映射成一个输出值,得到最终结果。
进一步的,在上述步骤e中,属性间融合顺序的不同会产生不同的效果。
本发明具有以下有益技术效果:
本发明针对滩坝砂油藏埋藏深、单层薄、横向变化大、地震分辨率低等难点,提出了利用三参数小波提高薄层识别能力的方法,并使用多子波检测技术对薄互层进行进一步识别。本发明能有效提高薄层识别能力,提高主频,展宽频带范围。经多子波变化后,能更加突出滩坝砂储层信息,经属性融合显示后,能大大提高后期储层描述的精度。
附图说明
下面结合附图与具体实施方式对本发明作更进一步说明:
图1为本发明的一种实施方式的流程框图。
图2示出了本发明中使用的三参数小波函数,通过参数调节,可以得到不同的小波函数。图2a选用参数为Λ=(0.5,2,0),图2b选用参数为Λ=(4,0.5,0)。
图3示出了本发明中通过三参数小波处理后得到的剖面对比,图3a为原始剖面,图3b是三参数小波处理剖面,选用参数Λ=(4,0.5,0),可见处理剖面的分辨率得到明显的提高。
图4示出了本发明中求取的多子波的波形特征,通过设定主频和频宽,可得到多个子波,图4展示了计算得到的三种子波。
图5示出了本发明中多子波变换后地震剖面,使用的地震数据为三参数小波处理后的数据。黑色波形显示为原始地震数据,背景图像为多子波处理变换剖面,两者进行了叠合显示,这样更突出处理后的效果。
图6示出了本发明应用在某一滩坝砂工区的应用效果。利用处理完的数据,进行属性分析等后期处理。提取了某层位0-20ms时窗内的均方根振幅属性,图6a是原始数据提取结果,图6b是三参数小波处理后提取结果,图6c是三参数小波处理后经多子波变换后的提取结果,图6d是将属性提取结果与层位等t0图进行融合显示的结果。
具体实施方式
结合附图,一种提高滩坝砂储层地震识别性能与描述精度的方法,包括以下步骤:
a、首先利用三参数小波对地震资料进行提高分辨率处理。三参数小波公式为
其中,t表示时间,i表示虚数,σ为小波函数的调制频率,τ为能量衰减因子,β为能量延迟因子,(σ,τ,β∈R,σ,τ≥0),Λ=(σ,τ,β)记参数σ,τ,β集合。上式中的p(Λ),q(Λ)和k(Λ)分别为
b、选取合适的三个小波参数Λ=(σ,τ,β),得到三参数小波函数φ(t,Λ)。设地震数据为X(t,j),j=1,2,...,M,其中M为道数、j表示第j道,按道进行褶积处理得到
s(t,j)=X(t,j)*φ(t,Λ),j=1,2,...,M
得到三参数小波处理后的数据。
c、对处理后的地震数据信号,首先估算滩坝砂储层位置的子波的频带范围[fmin,fmax],则可求出中心频率fc和频带宽度fw,设子波在该频带内的能量占子波总能量的λ2,则
将该式转化为求解以下方程
Aw=λw
其中,A={am,n}N×N,则得到
其中,N为子波样点数,Δt为采样间隔。取λ接近于1所对应的特征向量w,记为{w},一般为多个,{w}的频域能量限制在频带[fmin,fmax]以内,{w}为所要求取的多子波。
d、求取了标准子波以后,利用所求得的子波与地震道做相关处理
其中s为地震道,w为所求子波,C(t,j)便是相关信号,即为最后多子波变换数据。
e、对以上数据进行属性提取,对提取的属性进行多属性融合。二属性融合时,我们把二维颜色面板横坐标方向看作是一个属性,纵坐标方向看做是另外一个属性,这样通过两个属性的值的大小来选择二维面板里的颜色。多属性RGB显示技术本质上是一种图像处理技术,通过映射函数
IRGB(X)=S[IR(X),IG(X),IB(X)]
X=(x,y,z)代表在数据体中的位置,IR(X),IG(X),IB(X)分别代表Red、Green、Blue三种颜色上归一化的属性分量,通过三维颜色体,将三个属性分量映射成一个输出值,得到最终结果。
在上述步骤a中,Λ=(σ,τ,β)可以根据实际情况进行选择,我们这里选择Λ=(4,0.5,0)。
在上述步骤c中,求取的多子波一般是2~3个,后期处理时,应根据实际情况,选取合适的多子波进行分析。
在上述步骤e中,属性间融合顺序的不同会产生不同的效果。
本发明获取的最终处理剖面可参见图6。经对比研究可以发现,本发明能够有效圈定滩坝砂体发育范围,提高了地震识别的性能和后期描述的精度。
更进一步地,发明人还针对本发明中某一些要素发生变化时,这些变化对处理结果的影响作了更深入的研究。诸如:
一、本发明中三参数小波参数改变时的波形特征,不同的参数选择对应不同的波形特征,也会对后期的薄层识别能力产生影响。
二、本发明中主频和频宽有关要素改变时的多子波变化状态,根据设定不同的多子波主频和频宽,可得到多个子波,如图4。图5展示的是子波3的处理剖面,从剖面上看,多子波变换剖面比原始剖面有所改进,通过后期属性分析及其多属性融合技术,展示出了该技术在实际处理中良好的应用效果。
下面是本发明的具体应用例:
将本发明应用于某一滩坝砂工区,利用处理完的数据,进行属性分析等后期处理。提取了某层位0-20ms时窗内的均方根振幅属性,图6a是原始数据提取结果,图6b是三参数小波处理后提取结果,图6c是三参数小波处理后经多子波变换后的提取结果,图6d是将图6c的属性提取结果与层位等t0图进行融合显示的结果。从分析结果可知,图6c的计算结果能比较好的圈定滩坝砂体发育范围,通过属性融合显示,应用效果比较明显。
上述方式中未述及的有关技术内容采取或借鉴已有技术即可实现。
需要说明的是,在本说明书的教导下本领域技术人员还可以做出这样或那样的容易变化方式,诸如等同方式,或明显变形方式。上述的变化方式均应在本发明的保护范围之内。
Claims (3)
1.一种提高滩坝砂储层地震识别性能与描述精度的方法,其特征在于包括以下步骤:
a、首先利用三参数小波对地震资料进行提高分辨率处理;三参数小波公式为
其中,t表示时间,i表示虚数,σ为小波函数的调制频率,τ为能量衰减因子,β为能量延迟因子,(σ,τ,β∈R,σ,τ≥0),Λ=(σ,τ,β)记参数σ,τ,β集合;上式中的p(Λ),q(Λ)和k(Λ)分别为
b、根据选取的三个小波参数Λ=(σ,τ,β),得到三参数小波函数φ(t,Λ);设地震数据为X(t,j),j=1,2,...,M,其中M为道数、j表示第j道,按道进行褶积处理得到
s(t,j)=X(t,j)*φ(t,Λ),j=1,2,...,M
得到三参数小波处理后的数据;
c、对处理后的地震数据信号,首先估算滩坝砂储层位置的子波的频带范围[fmin,fmax],则可求出中心频率fc和频带宽度fw,设子波在该频带内的能量占子波总能量的λ2,则
将该式转化为求解以下方程
Aw=λw
其中,A={am,n}N×N,则得到
其中,N为子波样点数,Δt为采样间隔;取λ接近于1所对应的特征向量w,w即为求取的子波,其集合记为{w},{w}的频域能量限制在频带[fmin,fmax]以内,{w}即为所要求取的多子波;
d、求取了标准子波以后,利用所求得的子波与地震道做相关处理
其中s为地震道,w为所求子波,C(t,j)便是相关信号,即为最后多子波变换数据;
e、对以上数据进行属性提取,对提取的属性进行多属性融合;
上述步骤e中,所述多属性融合为二属性融合;即把二维颜色面板横坐标方向看作是一个属性,纵坐标方向是另外一个属性,这样通过两个属性的值的大小来选择二维面板里的颜色,多属性RGB显示技术本质上是一种图像处理技术,通过映射函数
IRGB(X)=S[IR(X),IG(X),IB(X)]
X=(x,y,z)代表在数据体中的位置,IR(X),IG(X),IB(X)分别代表红、绿、蓝三种颜色上归一化的属性分量,通过三维颜色体,将三个属性分量映射成一个输出值,得到最终结果。
2.根据权利要求1所述的一种提高滩坝砂储层地震识别性能与描述精度的方法,其特征在于:上述步骤a中,选择Λ=(4,0.5,0)。
3.根据权利要求1所述的一种提高滩坝砂储层地震识别性能与描述精度的方法,其特征在于:上述步骤c中,求取的多子波是2~3个。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210405357.9A CN102928872B (zh) | 2012-10-23 | 2012-10-23 | 一种提高滩坝砂储层地震识别性能与描述精度的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210405357.9A CN102928872B (zh) | 2012-10-23 | 2012-10-23 | 一种提高滩坝砂储层地震识别性能与描述精度的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102928872A CN102928872A (zh) | 2013-02-13 |
CN102928872B true CN102928872B (zh) | 2016-04-06 |
Family
ID=47643712
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210405357.9A Expired - Fee Related CN102928872B (zh) | 2012-10-23 | 2012-10-23 | 一种提高滩坝砂储层地震识别性能与描述精度的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102928872B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103513286B (zh) * | 2013-10-15 | 2016-05-18 | 中国石油大学(华东) | 一种地质模型约束下的滩坝结构单元判别方法 |
CN103969685B (zh) * | 2014-04-23 | 2016-07-13 | 长江大学 | 一种薄互层地震信号的处理方法 |
CN104698497B (zh) * | 2015-04-09 | 2015-10-14 | 中国石油大学(华东) | 一种砂泥岩薄互层层位的精细标定方法 |
CN111562615B (zh) * | 2020-06-02 | 2022-02-11 | 中国石油化工股份有限公司 | 退积型滩坝砂体发育区确定方法、油藏识别方法及装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080066392A1 (en) * | 2006-09-15 | 2008-03-20 | Bradford Tyler Sorensen | Variable floorplan shelters for previously unbuildable types of land |
-
2012
- 2012-10-23 CN CN201210405357.9A patent/CN102928872B/zh not_active Expired - Fee Related
Non-Patent Citations (3)
Title |
---|
小波变换在提高地震资料分辨率中的应用综述;王斌等;《内蒙古石油化工》;20080228(第04期);全文 * |
小波变换域地震资料瞬时频率分析方法;高静怀等;《石油物探》;20071115(第06期);第1-3部分 * |
滩坝砂油藏储层预测方法研究;张瑞芳;《中国优秀硕士学位论文全文数据库》;20120531;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN102928872A (zh) | 2013-02-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9964654B2 (en) | Seismic attribute color model transform | |
Du et al. | Seismic facies analysis based on self-organizing map and empirical mode decomposition | |
CN104020492A (zh) | 一种三维地震资料的保边滤波方法 | |
CN104537678B (zh) | 一种从单幅遥感图像中去除云雾的方法 | |
CN108985304B (zh) | 一种基于浅剖数据的沉积层结构自动提取方法 | |
CN102928872B (zh) | 一种提高滩坝砂储层地震识别性能与描述精度的方法 | |
CN104766090B (zh) | 一种基于bemd和sofm的探地雷达数据可视化方法 | |
CN107255831A (zh) | 一种叠前频散属性的提取方法 | |
CN107609573A (zh) | 基于低秩分解和空谱约束的高光谱图像时变特征提取方法 | |
CN103116881A (zh) | 基于PCA与Shearlet变换的遥感图像融合方法 | |
CN107179550B (zh) | 一种数据驱动的地震信号零相位反褶积方法 | |
CN103364832A (zh) | 一种基于自适应最优核时频分布的地震衰减定性估计方法 | |
CN107356967A (zh) | 一种压制地震资料强屏蔽干扰的稀疏优化方法 | |
CN104252549A (zh) | 一种基于克里金插值的分析布井方法 | |
CN107193046A (zh) | 一种基于谱反演的砂体厚度预测方法及系统 | |
CN103136733B (zh) | 基于多尺度图像分割和颜色传递的遥感图像彩色增强方法 | |
CN106291682A (zh) | 一种基于基追踪方法的叠后声波阻抗反演方法 | |
CN112324422B (zh) | 一种电成像测井缝洞识别方法、系统及孔隙结构表征方法 | |
CN105513060A (zh) | 一种视觉感知启发的高分辨率遥感图像分割方法 | |
CN105259571A (zh) | 一种地层倾角检测方法 | |
WO2019169932A1 (zh) | 基于粗糙集理论与色相法的多波油气地震响应表征方法 | |
CN105319593A (zh) | 基于曲波变换和奇异值分解的联合去噪方法 | |
CN102866425A (zh) | 一种基于盲源地震信号平稳叠加模型的盲分离方法 | |
CN110058298B (zh) | 一种三维地质体空间插值方法及系统 | |
CN105093312B (zh) | 基于频率域多阶微分的地震相对波阻抗预测方法和装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160406 |