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

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

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