[發(fā)明專利]一種井震聯(lián)合初始波阻抗建模方法

有效

• 申請(qǐng)?zhí)枺?/td>: 201310354864.9
• 申請(qǐng)日: 2013-08-15
• 公開(kāi)（公告）號(hào): CN103399346A
• 公開(kāi)（公告）日: 2013-11-20
• 發(fā)明（設(shè)計(jì)）人: 彭真明;范弘毅;周晗;王圣川;孔德輝;張爽
• 申請(qǐng)（專利權(quán)）人: 電子科技大學(xué)
• 主分類號(hào): G01V1/30
• 分類號(hào): G01V1/30;G06F19/00
• 代理公司: 成都華典專利事務(wù)所(普通合伙) 51223
• 代理人: 徐豐;楊保剛

1.一種井震聯(lián)合初始波阻抗建模方法，其特征在于，該建模方法包括如下步驟：?

步驟一、從文件讀入地震數(shù)據(jù)體、測(cè)井曲線以及層位數(shù)據(jù)；?

步驟二、對(duì)步驟一讀入的地震數(shù)據(jù)體進(jìn)行抽道操作，間隔整體道數(shù)的5%-10%抽取一道，得到一組新的地震道，在新的地震道上建立偽井道集；?

步驟三、對(duì)步驟二抽取出的偽井道集逐道進(jìn)行重采樣和閾值處理，選擇線性插值方法來(lái)進(jìn)行測(cè)井?dāng)?shù)據(jù)的重采樣，所述重采樣的測(cè)井阻抗通過(guò)下式求取反射系數(shù)：?

其中，r為反射系數(shù)點(diǎn)，I為測(cè)井曲線上的阻抗點(diǎn)，n為采樣位置；并將反射系數(shù)進(jìn)行稀疏化或平滑得到反射系數(shù)序列，稀疏公式為：?

步驟四：將步驟三中得到的反射系數(shù)序列進(jìn)行反射界面與地震道中振幅點(diǎn)的匹配；?

步驟五：在得到了匹配完成的偽井之后，通過(guò)加權(quán)來(lái)將每一道偽井外推進(jìn)行逐道的建模；所述加權(quán)為距離反比的加權(quán)，其公式為：?

其中，Z_p為目標(biāo)點(diǎn)的模型阻抗，Z_i為各個(gè)井上相應(yīng)點(diǎn)的阻抗，d_i為目標(biāo)點(diǎn)位置和各個(gè)井之間的距離，定義為：xline為數(shù)據(jù)的CDP?號(hào)。?

步驟六：將步驟五所獲得的模型進(jìn)行平滑整理，所用的平滑方式為五點(diǎn)平滑的空域?yàn)V波器，即把N個(gè)數(shù)據(jù)看成一個(gè)隊(duì)列，隊(duì)列的長(zhǎng)度固定為N，每進(jìn)行一次新的采樣，把新采樣結(jié)果放入隊(duì)尾，去掉排隊(duì)首的一個(gè)數(shù)據(jù)，對(duì)于不同的數(shù)據(jù)加以相同的權(quán)，這樣每進(jìn)行一次滑動(dòng)就可算得新的濾波值，則第n次采樣經(jīng)濾波后的輸出為：?

其中，N為滑動(dòng)平均項(xiàng)數(shù)，這里設(shè)為5點(diǎn)滑動(dòng)平均，X_n-i是未經(jīng)濾波的第n-i次采樣值；?

步驟七：將獲得的每一道模型進(jìn)行排列和顯示，之后交由后續(xù)的反演流程進(jìn)行反演處理。?

2.根據(jù)權(quán)利要求1所述的井震練習(xí)初始波阻抗建模方法，其特征在于，所述步驟二中偽井的數(shù)量是一個(gè)待選擇的參數(shù)。

3.根據(jù)權(quán)利要求1所述的井震練習(xí)初始波阻抗建模方法，其特征在于，所述步驟四中的匹配，其匹配方法如下：?

①先對(duì)應(yīng)搜索絕對(duì)值最大的峰谷值點(diǎn)，搜索到之后，進(jìn)行反射界面的搬移，搬移過(guò)程中地震道和反射系數(shù)道被分成三段，對(duì)這三段繼續(xù)進(jìn)行搜索絕對(duì)值最大的峰谷值點(diǎn)的操作；當(dāng)全部大的反射系數(shù)進(jìn)行匹配之后，整個(gè)地震記錄被分為若干段；?

②對(duì)分成段的地震記錄進(jìn)行對(duì)應(yīng)反射系數(shù)段上的插值處理得到分辨率較高的偽井模型；?

③將進(jìn)行了峰值匹配的偽井通過(guò)遞推公式遞推得到偽井的阻抗序列。?