用于同相軸自動拾取的傾斜疊加峰值振幅處邊緣檢測法
【技術(shù)領(lǐng)域】
[0001] 本發(fā)明涉及地震勘探技術(shù)領(lǐng)域，具體為一種用于同相軸自動拾取的傾斜疊加峰值 振幅處邊緣檢測法。
【背景技術(shù)】
[0002] 走時是衰減和速度層析成像的重要參數(shù)之一，我們需要拾取每個同相軸不同偏移 距處反射波的走時，包括大偏移距處的同相軸。可以看出，走時拾取的前提是同相軸拾取， 同相軸拾取也是地震解釋的重要步驟之一。但由于噪聲的影響，大偏移距處道集的信噪比 低，拾取連續(xù)同相軸比較困難。為了拾取正確的反射波走時，我們需要在大偏移距范圍內(nèi)拾 取連續(xù)同相軸，以便從記錄上讀取對應(yīng)射線的走時。同相軸在地震數(shù)據(jù)解釋中起重要作用。 現(xiàn)有的同相軸拾取方法有人工拾取方法及自動拾取方法等，人工拾取方法計算速度慢、精 度低，且受到主觀因素的影響；自動拾取方法有互相關(guān)法、模式識別法、神經(jīng)網(wǎng)絡(luò)法等，這些 方法很難拾取連續(xù)同相軸，且抗噪性差。

【發(fā)明內(nèi)容】

[0003] 針對現(xiàn)有技術(shù)中存在的問題，本發(fā)明提供一種信號的信噪比高，邊緣清晰準(zhǔn)確，完 整度及連續(xù)性好的，用于同相軸自動拾取的傾斜疊加峰值振幅處邊緣檢測法。
[0004] 本發(fā)明是通過以下技術(shù)方案來實現(xiàn)：
[0005] 本發(fā)明用于同相軸自動拾取的傾斜疊加峰值振幅處邊緣檢測法，包括如下步驟，
[0006] 步驟1，將疊前CSP數(shù)據(jù)用小波變換得到信號的解析部分；
[0007] 步驟2,根據(jù)步驟1得到的解析部分得到信號的瞬時振幅；
[0008] 步驟3,在瞬時振幅的剖面上進行局部線性Radon變換，拾取瞬時振幅的峰值振 幅，得到傾斜疊加峰值振幅剖面；
[0009] 步驟4,基于傾斜疊加峰值振幅剖面得到邊緣檢測剖面，用于在邊緣檢測剖面上拾 取同相軸及走時。
【主權(quán)項】
1. 用于同相軸自動拾取的傾斜疊加峰值振幅處邊緣檢測法，其特征在于，包括如下步 驟， 步驟1，將疊前CSP數(shù)據(jù)用小波變換得到信號的解析部分； 步驟2,根據(jù)步驟1得到的解析部分得到信號的瞬時振幅； 步驟3,在瞬時振幅的剖面上進行局部線性Radon變換，拾取瞬時振幅的峰值振幅，得 到傾斜疊加峰值振幅剖面； 步驟4,基于傾斜疊加峰值振幅剖面得到邊緣檢測剖面，用于在邊緣檢測剖面上拾取同 相軸及走時。
2. 根據(jù)權(quán)利要求1所述的用于同相軸自動拾取的傾斜疊加峰值振幅處邊緣檢測法，其 特征在于，步驟1中，將疊前CSP數(shù)據(jù)s(t)用小波變換得到信號的解析部分H[s(t)]如下，
其中，S(b，a)是地震道集數(shù)據(jù)的小波變換，心(叫為小波函數(shù)g(t)的Fourier變換的 實部，a為尺度因子，b為平移因子。
3. 根據(jù)權(quán)利要求2所述的用于同相軸自動拾取的傾斜疊加峰值振幅處邊緣檢測法，其 特征在于，s(t)關(guān)于g(t)的小波變換定義為：
式中，t,bGR，a> 0;g(t)GL1 (R,dt)nL2 (R,dt)，卵）是g⑴的復(fù)共軛。
4. 根據(jù)權(quán)利要求1所述的用于同相軸自動拾取的傾斜疊加峰值振幅處邊緣檢測法，其 特征在于，步驟2中，所述的瞬時振幅如下， A(t) =Is(t)+i?H[s(t)]I； 其中，H[s(t)]為用小波變換計算的信號的解析部分，A(t)為瞬時振幅的模，s(t)為疊 前CSP數(shù)據(jù)。
5. 根據(jù)權(quán)利要求1所述的用于同相軸自動拾取的傾斜疊加峰值振幅處邊緣檢測法，其 特征在于，步驟3中所述的局部線性Radon變換方法具體步驟如下， 步驟3. 1，在瞬時振幅剖面上選定參考道； 步驟3. 2,對于參考道上的某個時間截距T」，將參考道附近的幾道沿np個具有不同斜 率Pj(j= 1，2,…，np)的直線進彳丁置加，斜率以Ap為間隔米樣； 步驟3. 3,計算該時間截距處瞬時振幅沿不同方向疊加的和，將該疊加值記錄在T -P坐標(biāo)軸相應(yīng)的位置（Tj，Pj)上，當(dāng)選取的疊加斜率與同相軸的斜率接近或相等時，t-X域中 的記錄沿該直線的疊加值最大； 步驟3. 4,將最大疊加值的平均值放置在t-x域中的對應(yīng)位置（TjlXm)上，能夠構(gòu)造一 個超道集以增加信噪比，稱該剖面為傾斜疊加峰值振幅剖面，對應(yīng)于最大疊加值的斜率同 樣放置在t-x域的對應(yīng)位置上，構(gòu)成射線梯度剖面；其中，t為時間，X為距離域。
6. 根據(jù)權(quán)利要求1所述的用于同相軸自動拾取的傾斜疊加峰值振幅處邊緣檢測法，其 特征在于，步驟3中，Radon變換頻率域的離散形式為：
的傅里葉變換，Xm為參考道，m(T，p)是時間域的Radon變換。
7.根據(jù)權(quán)利要求1所述的用于同相軸自動拾取的傾斜疊加峰值振幅處邊緣檢測法，其 特征在于，步驟4中得到的邊緣檢測剖面由仿真軟件Matlab中的edge函數(shù)利用canny微 分算子得到。
【專利摘要】本發(fā)明用于同相軸自動拾取的傾斜疊加峰值振幅處邊緣檢測法，包括如下步驟，1.將疊前CSP數(shù)據(jù)用小波變換得到信號的解析部分；2.根據(jù)步驟1得到的解析部分得到信號的瞬時振幅；3.在瞬時振幅的剖面上進行局部線性Radon變換，拾取瞬時振幅的峰值振幅，得到傾斜疊加峰值振幅剖面；4.基于傾斜疊加峰值振幅剖面得到邊緣檢測剖面，用于在邊緣檢測剖面上拾取同相軸及走時。在瞬時振幅剖面上進行傾斜疊加變換，因為瞬時振幅剖面相較原始地震記錄更加圓滑，包含更多的低頻分量，且可濾除部分高頻毛刺噪聲，而二者的同相軸位置相同。有效壓制了噪聲，克服了現(xiàn)有方法難以在大偏移距內(nèi)拾取連續(xù)同相軸的困難，且操作方便，實現(xiàn)了同相軸自動拾取。
【IPC分類】G01V1-28
【公開號】CN104714251
【申請?zhí)枴緾N201510132019
【發(fā)明人】高靜懷, 趙靜
【申請人】西安交通大學(xué)
【公開日】2015年6月17日
【申請日】2015年3月24日