基于迭代三参数小波变换的缝洞发育带检测方法
【专利摘要】本发明公开了一种基于迭代三参数小波变换的缝洞发育带检测方法，包括：采用三参数小波变换代替常规的小波变换进行计算，对2维地震切片的横向和纵向同时进行边缘检测处理，采用横向和纵向的边缘检测的结果的均方根作为最终的边缘检测结果，再进行2次小波变换，并对小波变换的结果数组进行保边去噪处理，得到最终的用于边缘检测的结果数组，将结果数组用于描述缝洞体发育带。本发明采用三参数小波变换代替常规的小波变换进行计算，克服现有技术对于地震资料匹配性较差的缺陷，提高边缘检测的精度。
【专利说明】基于迭代三参数小波变换的缝洞发育带检测方法
【技术领域】
[0001]本发明涉及一种油气地球物理勘探检测方法，尤其涉及一种基于迭代三参数小波变换的缝洞发育带检测方法。
【背景技术】
[0002]在油气地球物理勘探领域，边缘检测方法主要用于寻找缝洞发育带，其往往是有利的油气储存空间或运移通道。
[0003]自然界的岩石或岩层中普遍存在着孔隙、裂缝(隙)和溶蚀孔洞，其形态各异，尺度大小也相差悬殊。对于地震勘探而言，由于受分辨率的限制，无法识别出单个的孔、洞、缝，只能识别规模达到一定程度的裂缝发育带。事实上，单个的孔、洞、缝对油气的聚集所起的作用是微乎其微的，真正具有勘探、开发价值的是具有一定规模的缝洞发育带。
[0004]裂缝发育带是指相对围体岩层而言，其缝洞密度明显增大，并有一定延伸范围的岩体。因此，其在地震水平切片(或沿层切片)上，具有一定的分布范围和延伸方向。因此，可以采用各种边缘检测方法识别裂缝发育带。
[0005]在油气地球物理勘探领域，现有的基于小波变换的边缘检测方法主要有3类:(I)基于常规小波变换的方法，如黄捍东、胡光岷(2000)提出的多尺度边缘检测方法；(2)基于多尺度小波变换边缘检测，如贺振华、黄德济、文晓涛等(2007)提出的小波多尺度边缘检测；(3)基于噪音抑制的多尺度小波变换的方法，如文晓涛(2008)提出的小波域尺度积在裂缝检测中的应用。
[0006]上述3类方法的共同点:(I)都是基于常规的小波变换方法进行计算，其小波函数都是固定的小波函数，如Morlet小波等，即各个不同地区的地震资料都采用相同的小波函数；(2)都只针对2维切片的I个方向进行了处理(横向或纵向)。
[0007]上述3类方法的不同点:(I)第I类方法只采用单个尺度的小波变换结果进行边缘检测；(2)第2类方法在第I类方法的基础上进行了改进，采用两个相邻尺度的小波变换结果的乘积进行边缘检测；(3)第3类方法在第2类方法基础上进行了改进，引入了保边去噪方法进行预处理，即在边缘检测之前，先采用保边去噪方法对2维切片进行去噪处理，然后再采用第2类方法进行边缘检测。
[0008]但在实际应用中，现有技术还存在如下问题:
(I)现有技术都采用常规的小波变换方法(采用固定的母小波，如Morlet小波)进行计算，而不同地区的地震资料的地震子波往往具有一定的差异，其与地震资料的匹配性较差。
[0009](2)现有技术都只针对2维切片的I个方向进行了处理(横向或纵向)，边缘检测的精度较低；
(3)上述第I类方法只采用单个尺度的小波变换的结果进行边缘检测，其边缘检测的精度最低；
(4)上述第I类方法和第2类方法都没有考虑对噪声进行预处理，其抑制噪声的能力较差； (5)上述第3类方法虽然在边缘检测之前先进行了保边去噪的预处理，但是其抑制噪声的能力仍有待提闻。
【发明内容】
[0010]本发明的目的在于克服现有技术存在的上述问题，提供一种基于迭代三参数小波变换的缝洞发育带检测方法。本发明采用三参数小波变换代替常规的小波变换进行计算，克服现有技术对于地震资料匹配性较差的缺陷，提高边缘检测的精度。
[0011]为实现上述目的，本发明采用的技术方案如下:
一种基于迭代三参数小波变换的缝洞发育带检测方法，其特征在于:采用三参数小波变换代替常规的小波变换进行计算，对2维地震切片的横向和纵向同时进行边缘检测处理，采用横向和纵向的边缘检测的结果的均方根作为最终的边缘检测结果，再进行2次小波变换，并对小波变换的结果数组进行保边去噪处理，得到最终的用于边缘检测的结果数组，将结果数组用于描述缝洞体发育带。
[0012]所述方法具体包括如下步骤:
(O综合研究区地震资料的子波分析、时频分析资料，确定三参数小波
的三个参数€Γ、Τ和,，并得到用于计算的三参数小波函数伊'￠),三参数小波的计算公式为:
公式一:炉_(?) ={ΚΛ)[?κ(β?) — *(A)] +*?(A)sin(csf)}
公式二: *(Λ) = ^,/4" cas(pd) +...KA)
~2r Im ? Γ —a
公式三:ρ(Α) — — ^ 4 β e 如 CQS2(^fer).f 1-■■■■.ψ
L
1/4 ? Γ *?Λ
公式四:g(A)= — Η 4 ? χ -e 31 s?+ ^ '
L ΛΓ J
式中，/⑩是能量有限或平方可积信号H是尺度因子，b是位移因子，-m是以0.，r，於为参数的小波函数(即三参数小波)，?Γ是调制频率，T是能量衰减因子，於是能量延迟因子；
(2)将二维地震切片记为,其中hit为incline方向，ret的取值从I至Jfc, Jfc为切片的InLine总数呀为Crossline方向，V的取值从I至雄，雄为切片的CrossLine总数；针对目的层切片，分两个尺度对该切片进行边缘检测处理；
(3)设定离散小波变换的尺度■/=1，分两个方向(Incline方向和Crossline方向)对目的层切片进行处理；
(3-1)提取Inline号=1的I维数组，即SkceQ,=^ffy ;令该I维数组为
~......_y. (3-2)采用公式一至公式四中的三参数小波，计算步骤(3-1)中的I维数组SHcBl_X(ny)的二进小波变换，得到当前尺度下的细节分量数组(HF = L-=j%);并将I维数组沒kx—置—找卿}写入2为数组SUasl X EQ,哗)ιψ = I,.(3-3)更换Inline号,重复步骤(3_1)至步骤(3_2),直到所有的Inline号计算完成，得到当前等时切片的边缘检测的结果数组—=；
(3-4)提取Crossline号=1的I维数组，即;令该I维数组为 T(nx^TEc— I,~,A& ;
(3-5)采用公式一致公式四中的三参数小波，计算步骤(3)和(4)中的I维数组的二进小波变换，得到当前尺度下的细节分量数组_Γ_/Xhj?)(呼=.1.....为、;并将将I维数组5Κ(?_Γ _/】(呼)写入2为数组SSScsl Y_JS(Xτψ) τψ -1r^fNy.(3-6)更换C1ssline号，重复步骤(3-4广步骤(3-5)，直到所有的Crossline号计算完成，从而得到当前等时切片的边缘检测的结果数组SUasl — Y_E(ja:，mf) wn: = 1,.Zfo-Hy = I, ；
(3-7)计算当前尺度下的小波变换的结果，即上述2个方向(Inline方向和Crossline方向)的计算结果的均方根，即下式公式五公式五:
ss*?zici^ir)=?=1?施寒=1__
(4)设定离散小波变换的尺度/ = 2 ;重复步骤(3-1)至(3-7)，得到尺度_/ = 2时的小波变换的结果；
(5)计算相邻尺度J_和/+I的小波变换的结果的乘积，
公式六:SSfc* Xixx;呼)二 Stsas ElijsxLj ixp) x SSics- E2(ms^ wg/)- ■■■■■W —■■■■■■■.ffy
(6)针对步骤(5)中的2维数组S&K—ICmiiF),重复步骤(3)至步骤(5);得到去噪处理后的 2 维数组，记为Sfee—2(B^BjKXnxr = l…；
(7)对步骤(6)中的2维数组S0?_2(intiiF),根据噪音情况进行多次迭代的2维保边去噪处理，得到最终的用于边缘检测的结果数组，记为
_ Hr =呼=1?抱。
[0013]采用本发明的优点在于:
一、本发明是一种基于迭代三参数小波变换的缝洞发育带检测方法，与现有技术相比，本发明采用三参数小波变换代替常规的小波变换进行计算，克服现有技术对于地震资料匹配性较差的缺陷；对2维地震切片的横向和纵向同时进行边缘检测处理，采用横向和纵向的边缘检测的结果的均方根作为最终的边缘检测结果，提高现有技术的边缘检测的精度；进一步提高现有技术的边缘检测的精度，采用迭代处理，即进行2次小波变换计算；提高抑制噪声的能力，采用后去噪的计算流程，即对小波变换的结果数组进行保边去噪处理；其得到的结果精度更高，且抑制噪声的能力更强，具有推广使用的价值。
[0014]二、本发明能够预测地震资料边缘变化的精度，并且抑制噪声能力较强，与通常的小波变换边缘检测流程及效果不同。
[0015]三、通过本方法计算的边缘检测结果可以用于描述地质体的缝洞体发育参数，且本发明限制在地震勘探中应用，地震勘探资料分为叠前、叠后两种，本发明只限制在叠后资料中。
[0016]四、本发明采用三参数小波变换代替常规的小波变换进行计算，克服现有技术对于地震资料匹配性较差的缺陷，提高边缘检测的精度。
[0017]五、本发明对目的层切片的两个方向(Inline方向及Crossline方向)同时进行处理，进一步提闻边缘检测的精度。
[0018]六、本发明对缝洞发育带的边缘检测而言，采用单一的小尺度C/ =1)的离散小波变换的结果能较准确地指示边缘，但是其受噪声影响较大；采用单一的较大尺度〔J =Z)的离散小波变换的结果能较好地抑制噪声，但是指示边缘的能力较差。为此，本专利采用公式六所示的两种尺度小波变换的乘积进行缝洞发育带的边缘检测，既可以消除噪声的影响，又可以准确地指示边缘。
[0019]七、本发明结合多次迭代的保边去噪对边缘检测后的结果进行处理，进一步压制噪音的影响。
【具体实施方式】
[0020]本发明是要提供一种基于迭代三参数小波变换的缝洞发育带检测方法，它针对二维地震切片进行处理，用于缝洞发育带的识别，有效地指导油气地球物理勘探的储层预测。
[0021]本发明可以达到如下目的:
(1)采用三参数小波变换代替常规的小波变换进行计算，克服现有技术对于地震资料匹配性较差的缺陷；
(2)对2维地震切片的横向和纵向同时进行边缘检测处理，采用横向和纵向的边缘检测的结果的均方根作为最终的边缘检测结果，提高现有技术的边缘检测的精度；
(3)进一步提高现有技术的边缘检测的精度，采用迭代处理，即进行2次小波变换计算(即第I次的结果作为第2次的输入)；
(4)提高现有技术的第3类方法的抑制噪声的能力，采用后去噪的计算流程，即对小波变换的结果数组进行保边去噪处理。
[0022]与现有技术相比，本发明的创新点如下:
(1)采用三参数小波进行小波变换计算，可以有效地适应各个地区的地震资料，即小波函数可以灵活的修改或改变以匹配研究区的地震资料；
(2)对小波变换后的结果数组进行保边去噪处理，而不是在小波变换之前进行去噪处理，可以更有效地去除噪声，并提高边缘检测的精度；
(3)采用迭代的思想，进行2次计算，即将第一次小波变换的结果数组作为第二次小波变换的输入数组，可以进一步有效地提高了边缘检测的精度及抑制噪声。
[0023]概括说明本发明: 一种基于迭代三参数小波变换的缝洞发育带检测方法，采用三参数小波变换代替常规的小波变换进行计算，对2维地震切片的横向和纵向同时进行边缘检测处理，采用横向和纵向的边缘检测的结果的均方根作为最终的边缘检测结果，再进行2次小波变换，并对小波变换的结果数组进行保边去噪处理，得到最终的用于边缘检测的结果数组，将结果数组用于描述缝洞体发育带。将结果数组用于描述缝洞体发育带可采用现有技术。
[0024]下面对本发明作进一步说明:
本发明包括以下步骤:
(O综合研究区地震资料的子波分析、时频分析资料，确定三参数小波伊’ω的三个参数和於，并得到用于计算的三参数小波函数#'￠0 ,三参数小波的计算公式为:
公式:φ'{￡) = {KA)[cw<oi) - Jfc(A)] +^(A) Sb(OE)I
公式::Jc(A) = casim + *
p(.Aj
^ r I~ ?._ —i~1/2
2.T"w, ~
公式三.::/?(Λ) = — -1 4 e.21 —e & msP(/fer) +?
π
严?Γ上广
公式四:g(A) = H 4 β —Sr Sn^GScO + l-1T* ητ <■
κ _
式中，/?是能量有限或平方可积信号a是尺度因子，k是位移因子，Sm是以0.，r岁力参数的小波函数(即三参数小波)，是调制频率，Γ是能量衰减因子，於是能量延迟因子。
[0025]公式一所示的三参数小波通过三个参数可以灵活的构建各种类型的小波基函数(即任意频率或任意相位或复杂波形)，可以灵活的匹配各个地区的地震资料(注:不同的地区，由于地表激发条件或地下地质构造的区别，其地震子波往往存在一定的差异)。此外，当
CT= 5Λ, T = OJ和於= 0-0时,三参数小波炉'CO与常规的Morlet小波一致。
[0026](2)将二维地震切片记为数组，其中nr为incline方向，ret的取值从I至Al., At为切片的InLine总数为Crossline方向，》!F的取值从I至抑，抑为切片的CrossLine总数。针对目的层切片，分两个尺度对该切片进行边缘检测处理。
[0027](3)设定离散小波变换的尺度J=7I ,分两个方向(Incline方向和Crossline方向)对目的层切片进行处理；
(3-1)提取Inline号=1的I维数组，即版机呼)，呼:;令该I维数组为SSiicsl二X_;
(3-2)采用公式一至公式四中的三参数小波，计算步骤(3-1)中的I维数组Jr(I9T)的二进小波变换，得到当前尺度下的细节分量数组I—(?Τ = 1?Λ&τ);并将 I 维数组JTjl(IIF)写入 2 为数组Siicel X= 1-Jfy ；
(3-3)更换Inline号,重复步骤(3_1)至步骤(3_2),直到所有的Inline号计算完成，得到当前等时切片的边缘检测的结果数组呼3 ^-,Aftiy = L......^Ny.(3-4)提取Crossline号=1的I维数组，即滿εφα九抵= Hiic ;令该I维数组为SKc:￡l_Y(?r),?ir = I,;
(3-5)采用公式一致公式四中的三参数小波，计算步骤(3)和(4)中的I维数组Sfc?l_F0y)的二进小波变换，得到当前尺度下的细节分量数组驗1—Y—nc—y;并将将I维数组驗J_mw)写入2为数组
JT —iBPLiy) rgr — Ia,-■■■_ Jtiy.(3-6)更换Crossline号，重复步骤(3-4广步骤(3-5)，直到所有的Crossline号计算完成，从而得到当前等时切片的边缘检测的结果数组
Sficei r JE1(I^1Jw) ητ=】，—」_Α6? w = l._.(3-7)计算当前尺度下的小波变换的结果，即上述2个方向(Inline方向和Crossline方向)的计算结果的均方根，即下式公式五公式五:
SSw_J3Cie;>ipp) = ^Μιχ1_Χ_Ε31?ε,?ψ) + ?οβ1_Γ孤== 1?-,Λ5τ
(4)设定离散小波变换的尺度J_ = 2 ;重复步骤(3-1)至(3-7)，得到尺度7 = 2时的小波变换的结火5;
(5)计算相邻尺度■/和/+I的小波变换的结果的乘积，
公」弋六.SfiCif XiEtjHjOM ==
(6)针对步骤(5)中的2维数组重复步骤(3)至步骤(5);得到去噪处理后的2维数组，记为^36￡?_2(叹呼Xnr = I一取呼=1?抱
(7)对步骤(6)中的2维数纟|[52scer_2(H^njr),根据噪音情况进行多次迭代的2维保边去噪处理，得到最终的用于边缘检测的结果数组，记为Shceijix^TgrXrax = Ij JNx^ny =0
[0028]对缝洞发育带的边缘检测而言，采用单一的小尺度(J =1)的离散小波变换的结果能较准确地指示边缘，但是其受噪声影响较大；采用单一的较大尺度(■/ =2)的离散小波变换的结果能较好地抑制噪声，但是指示边缘的能力较差。为此，本专利采用公式六所示的两种尺度小波变换的乘积进行缝洞发育带的边缘检测，既可以消除噪声的影响，又可以准确地指示边缘。
[0029]本专利在塔钟某地区裂缝检测中的应用，与实钻资料对比，吻合率达94%。
[0030]对缝洞发育带的边缘检测而言，采用单一的小尺度C/ =1)的离散小波变换的结果能较准确地指示边缘，但是其受噪声影响较大；采用单一的较大尺度(■/ =2)的离散小波变换的结果能较好地抑制噪声，但是指示边缘的能力较差。为此，本专利采用公式六所示的两种尺度小波变换的乘积进行缝洞发育带的边缘检测，既可以消除噪声的影响，又可以准确地指示边缘。
[0031]应用实践表明，本发明采用三参数小波进行小波变换计算，能有效地适应各地区的地震资料；对小波变换的结果进行保边去噪处理，可更有效地去除噪音；对2维地震切片的横向和纵向的两个方向同时进行边缘检测处理，采用横向和纵向的边缘检测的结果的均方根作为最终的边缘检测结果，且采用迭代的思想，进行2次计算，即将第一次小波变换的结果数组作为第二次小波变换的输入数组，可以进一步有效地提高了边缘检测的精度及抑制噪声。
【权利要求】
1.一种基于迭代三参数小波变换的缝洞发育带检测方法，其特征在于:采用三参数小波变换代替常规的小波变换进行计算，对2维地震切片的横向和纵向同时进行边缘检测处理，采用横向和纵向的边缘检测的结果的均方根作为最终的边缘检测结果，再进行2次小波变换，并对小波变换的结果数组进行保边去噪处理，得到最终的用于边缘检测的结果数组，将结果数组用于描述缝洞体发育带。
2.根据权利要求1所述的基于迭代三参数小波变换的缝洞发育带检测方法，其特征在于:所述方法具体包括如下步骤: (O综合研究区地震资料的子波分析、时频分析资料，确定三参数小波炉'CO的三个参数U、T和於，并得到用于计算的三参数小波函数炉三参数小波的计算公式为: 公式一:φ\?) = {ΚΛ)[ο?<ο?) -fc(A)] +ig{A^aO] 公式二 MA) = ^ coscm + ?
ρ(.Α)
^ Γ I~?._ —ι~1/2
窜.—Sr 公式三:ρ(Λ) = — ■' 4 β.21 — e Cos5QSfcr) +1"
I χ" 公式四:5(Α) = η 4 e 21 —e ^5(jiScr)+ 1-Sr^r -
K 式中，/?是能量有限或平方可积信号，是尺度因子，b是位移因子，#*'?是以0.，r岁力参数的小波函数(即三参数小波)，0.是调制频率，Γ是能量衰减因子，岁是能量延迟因子； (2)将二维地震切片记为数组,其中nr为incline方向，rec的取值从I至A&, Jfc为切片的InLine总数W为Crossline方向’耶的取值从I至脚,场?为切片的CrossLine总数；针对目的层切片，分两个尺度对该切片进行边缘检测处理； (3)设定离散小波变换的尺度/=1,分两个方向(Incline方向和Crossline方向)对目的层切片进行处理； (3-1)提取Inline号=1的I维数组，即呼；令该I维数组为 _JT(nyXny ==.(3-2)采用公式一至公式四中的三参数小波，计算步骤(3-1)中的I维数组的二进小波变换，得到当前尺度下的细节分量数组S&xl_X_IKw) =;并将 I 维数组写入 2 为数组Sik^1._.hST_.ElQt ny)ny — I，■■■■■■―泛.jffy - (3-3)更换Inline号,重复步骤(3_1)至步骤(3_2),直到所有的Inline号计算完成，得到当前等时切片的边缘检测的结果数组呼)ret； (3-4)提取Crossline号=1的I维数组，即nr=Uifc ;令该I维数组为^Kcfii r(Kr)31 KET =X_^Nx ; (3-5)采用公式一致公式四中的三参数小波，计算步骤(3)和(4)中的I维数组的二进小波变换，得到当前尺度下的细节分量数组SUcel_T_Π(ιψ} =;并将将 I 维数组 Sice_￥ _η^ψ}写入 2 为数组 Jqy — I,」坤; (3-6)更换C1ssline号，重复步骤(3-4广步骤(3-5)，直到所有的Crossline号计算完成，从而得到当前等时切片的边缘检测的结果数组Γ—￡(?η^呼)Hr-17_^Νχ；7ψ — χ_,Jftr.(3-7)计算当前尺度下的小波变换的结果，即上述2个方向(Inline方向和Crossline方向)的计算结果的均方根，即下式公式五公式五: Mcm._EMmc,w} =Jf蚤織> ?r:l-^Μ?ιψ =(4)设定离散小波变换的尺度J_ = 2 ;重复步骤(3-1)至(3-7)，得到尺度i = 2时的小波变换的结果.....现取呼)； (5)计算相邻尺度/和的小波变换的结果的乘积， 公式六:5Sfc?ir_Xwti;#iy)= (6)针对步骤(5)中的2维数组重复步骤(3)至步骤(5);得到去噪处理后的2维数组，记为5Kc?___20n;jgO，Hr二I一； (7)对步骤(6)中的2维数组Skce_2(n^mr),根据噪音情况进行多次迭代的2维保边去噪处理，得到最终的用于边缘检测的结果数组，记为?3?τβ(ιπο mr) mo = I Jftriy = I _rNy o
【文档编号】G01V1/28GK104166162SQ201410413732
【公开日】2014年11月26日 申请日期:2014年8月21日 优先权日:2014年8月21日
【发明者】徐明华, 刘俊海, 欧阳诚, 吴大奎, 陈琴 申请人:中国石油集团川庆钻探工程有限公司