用于同相轴自动拾取的倾斜叠加峰值振幅处边缘检测法
【技术领域】
[0001] 本发明涉及地震勘探技术领域，具体为一种用于同相轴自动拾取的倾斜叠加峰值 振幅处边缘检测法。
【背景技术】
[0002] 走时是衰减和速度层析成像的重要参数之一，我们需要拾取每个同相轴不同偏移 距处反射波的走时，包括大偏移距处的同相轴。可以看出，走时拾取的前提是同相轴拾取， 同相轴拾取也是地震解释的重要步骤之一。但由于噪声的影响，大偏移距处道集的信噪比 低，拾取连续同相轴比较困难。为了拾取正确的反射波走时，我们需要在大偏移距范围内拾 取连续同相轴，以便从记录上读取对应射线的走时。同相轴在地震数据解释中起重要作用。 现有的同相轴拾取方法有人工拾取方法及自动拾取方法等，人工拾取方法计算速度慢、精 度低，且受到主观因素的影响；自动拾取方法有互相关法、模式识别法、神经网络法等，这些 方法很难拾取连续同相轴，且抗噪性差。
【发明内容】
[0003] 针对现有技术中存在的问题，本发明提供一种信号的信噪比高，边缘清晰准确，完 整度及连续性好的，用于同相轴自动拾取的倾斜叠加峰值振幅处边缘检测法。
[0004] 本发明是通过以下技术方案来实现：
[0005] 本发明用于同相轴自动拾取的倾斜叠加峰值振幅处边缘检测法，包括如下步骤，
[0006] 步骤1，将叠前CSP数据用小波变换得到信号的解析部分；
[0007] 步骤2,根据步骤1得到的解析部分得到信号的瞬时振幅；
[0008] 步骤3,在瞬时振幅的剖面上进行局部线性Radon变换，拾取瞬时振幅的峰值振 幅，得到倾斜叠加峰值振幅剖面；
[0009] 步骤4,基于倾斜叠加峰值振幅剖面得到边缘检测剖面，用于在边缘检测剖面上拾 取同相轴及走时。
【主权项】
1. 用于同相轴自动拾取的倾斜叠加峰值振幅处边缘检测法，其特征在于，包括如下步 骤， 步骤1，将叠前CSP数据用小波变换得到信号的解析部分； 步骤2,根据步骤1得到的解析部分得到信号的瞬时振幅； 步骤3,在瞬时振幅的剖面上进行局部线性Radon变换，拾取瞬时振幅的峰值振幅，得 到倾斜叠加峰值振幅剖面； 步骤4,基于倾斜叠加峰值振幅剖面得到边缘检测剖面，用于在边缘检测剖面上拾取同 相轴及走时。
2. 根据权利要求1所述的用于同相轴自动拾取的倾斜叠加峰值振幅处边缘检测法，其 特征在于，步骤1中，将叠前CSP数据s(t)用小波变换得到信号的解析部分H[s(t)]如下，
其中，S(b，a)是地震道集数据的小波变换，心(叫为小波函数g(t)的Fourier变换的 实部，a为尺度因子，b为平移因子。
3. 根据权利要求2所述的用于同相轴自动拾取的倾斜叠加峰值振幅处边缘检测法，其 特征在于，s(t)关于g(t)的小波变换定义为：
式中，t,bGR，a> 0;g(t)GL1 (R,dt)nL2 (R,dt)，卵）是g⑴的复共轭。
4. 根据权利要求1所述的用于同相轴自动拾取的倾斜叠加峰值振幅处边缘检测法，其 特征在于，步骤2中，所述的瞬时振幅如下， A(t) =Is(t)+i?H[s(t)]I； 其中，H[s(t)]为用小波变换计算的信号的解析部分，A(t)为瞬时振幅的模，s(t)为叠 前CSP数据。
5. 根据权利要求1所述的用于同相轴自动拾取的倾斜叠加峰值振幅处边缘检测法，其 特征在于，步骤3中所述的局部线性Radon变换方法具体步骤如下， 步骤3. 1，在瞬时振幅剖面上选定参考道； 步骤3. 2,对于参考道上的某个时间截距T」，将参考道附近的几道沿np个具有不同斜 率Pj(j= 1，2,…，np)的直线进彳丁置加，斜率以Ap为间隔米样； 步骤3. 3,计算该时间截距处瞬时振幅沿不同方向叠加的和，将该叠加值记录在T -P坐标轴相应的位置（Tj，Pj)上，当选取的叠加斜率与同相轴的斜率接近或相等时，t-X域中 的记录沿该直线的叠加值最大； 步骤3. 4,将最大叠加值的平均值放置在t-x域中的对应位置（TjlXm)上，能够构造一 个超道集以增加信噪比，称该剖面为倾斜叠加峰值振幅剖面，对应于最大叠加值的斜率同 样放置在t-x域的对应位置上，构成射线梯度剖面；其中，t为时间，X为距离域。
6. 根据权利要求1所述的用于同相轴自动拾取的倾斜叠加峰值振幅处边缘检测法，其 特征在于，步骤3中，Radon变换频率域的离散形式为：
的傅里叶变换，Xm为参考道，m(T，p)是时间域的Radon变换。
7.根据权利要求1所述的用于同相轴自动拾取的倾斜叠加峰值振幅处边缘检测法，其 特征在于，步骤4中得到的边缘检测剖面由仿真软件Matlab中的edge函数利用canny微 分算子得到。
【专利摘要】本发明用于同相轴自动拾取的倾斜叠加峰值振幅处边缘检测法，包括如下步骤，1.将叠前CSP数据用小波变换得到信号的解析部分；2.根据步骤1得到的解析部分得到信号的瞬时振幅；3.在瞬时振幅的剖面上进行局部线性Radon变换，拾取瞬时振幅的峰值振幅，得到倾斜叠加峰值振幅剖面；4.基于倾斜叠加峰值振幅剖面得到边缘检测剖面，用于在边缘检测剖面上拾取同相轴及走时。在瞬时振幅剖面上进行倾斜叠加变换，因为瞬时振幅剖面相较原始地震记录更加圆滑，包含更多的低频分量，且可滤除部分高频毛刺噪声，而二者的同相轴位置相同。有效压制了噪声，克服了现有方法难以在大偏移距内拾取连续同相轴的困难，且操作方便，实现了同相轴自动拾取。
【IPC分类】G01V1-28
【公开号】CN104714251
【申请号】CN201510132019
【发明人】高静怀, 赵静
【申请人】西安交通大学
【公开日】2015年6月17日
【申请日】2015年3月24日