[發(fā)明專利]縱波和轉(zhuǎn)換波疊前聯(lián)合反演方法及系統(tǒng)有效
| 申請?zhí)枺?/td> | 201810857538.2 | 申請日: | 2018-07-31 |
| 公開(公告)號: | CN110780351B | 公開(公告)日: | 2021-11-23 |
| 發(fā)明(設(shè)計)人: | 陳天勝;張珺;季玉新 | 申請(專利權(quán))人: | 中國石油化工股份有限公司;中國石油化工股份有限公司石油勘探開發(fā)研究院 |
| 主分類號: | G01V1/50 | 分類號: | G01V1/50 |
| 代理公司: | 北京思創(chuàng)畢升專利事務(wù)所 11218 | 代理人: | 孫向民;廉莉莉 |
| 地址: | 100728 北*** | 國省代碼: | 北京;11 |
| 權(quán)利要求書: | 查看更多 | 說明書: | 查看更多 |
| 摘要: | |||
| 搜索關(guān)鍵詞: | 縱波 轉(zhuǎn)換 波疊前 聯(lián)合 反演 方法 系統(tǒng) | ||
1.一種縱波和轉(zhuǎn)換波疊前聯(lián)合反演方法,其特征在于,包括:
1)輸入縱波時間域的縱波角度道集和轉(zhuǎn)換波時間域的轉(zhuǎn)換波角度道集
2)提取各角度縱波子波Wpp和轉(zhuǎn)換波子波Wps;
3)建立縱波時間域的初始彈性參數(shù)模型其中,k為迭代次數(shù);
4)計算縱波時間域的縱波合成記錄誤差
5)計算時移算子矩陣Q和誤差反傳播算子矩陣QT;
6)計算轉(zhuǎn)換波時間域的轉(zhuǎn)換波合成記錄誤差
7)計算縱波時間域的縱波和轉(zhuǎn)換波的總誤差rk,如果總誤差小于或等于設(shè)定誤差門檻值,轉(zhuǎn)至步驟9),如果總誤差大于設(shè)定誤差門檻值,執(zhí)行步驟8);
8)更新縱波時間域的對數(shù)彈性參數(shù)模型,重復(fù)執(zhí)行步驟4)~步驟7);
9)輸出彈性參數(shù)反演結(jié)果
其中,為彈性參數(shù)模型,分別為在縱波時間域的對數(shù)縱波速度向量、對數(shù)橫波速度向量、對數(shù)密度向量;
在步驟8)中,采用誤差反傳播算子變權(quán)重共軛梯度迭代法更新縱波時間域的對數(shù)彈性參數(shù)模型,包括:
8.1)令其中,p(k)為彈性參數(shù)更新方向參數(shù),p(0)=r0;和是k次迭代的時移算子矩陣Q和誤差反傳播算子矩陣QT的函數(shù),即λ為縱波和轉(zhuǎn)換波誤差調(diào)節(jié)參數(shù);
8.2)計算與p(k)正交的更新方向參數(shù)p(k+1)=rk+1+αkp(k),
其中,
和是k+1次迭代的時移算子矩陣Q和誤反傳播算子矩陣QT的函數(shù),即
其中,為迭代次數(shù)為k時的彈性參數(shù)模型,為迭代次數(shù)為k+1時的彈性參數(shù)模型;
所述縱波時間域的縱波角度道集和所述轉(zhuǎn)換波時間域的轉(zhuǎn)換波角度道集分別通過公式(1)和公式(2)表示:
其中,Gpp=[WppAD WppBD WppCD],Gps=[0 WpsQED WpsQFD],
A(θi)、B(θi)、C(θi)為在縱波入射角為θi時縱波反射系數(shù)中的縱波速度變化率、橫波速度變化率和密度變化率的系數(shù),E(θi)、F(θi)為在縱波入射角為θi時轉(zhuǎn)換波反射系數(shù)中的橫波速度變化率和密度變化率的系數(shù),K為接收地震道數(shù)目;D為差分算子矩陣,
2.根據(jù)權(quán)利要求1所述的縱波和轉(zhuǎn)換波疊前聯(lián)合反演方法,其特征在于,進行反演的目標(biāo)函數(shù)為:
其中,λ為縱波和轉(zhuǎn)換波誤差調(diào)節(jié)參數(shù)。
3.根據(jù)權(quán)利要求2所述的縱波和轉(zhuǎn)換波疊前聯(lián)合反演方法,其特征在于,利用最小二乘法求解目標(biāo)函數(shù)(3),得到方程(4):
通過對所述方程(4)進行求解得到
4.根據(jù)權(quán)利要求1所述的縱波和轉(zhuǎn)換波疊前聯(lián)合反演方法,其特征在于,所述縱波時間域的縱波合成記錄誤差的表達式為:
其中,為迭代次數(shù)為k時的彈性參數(shù)模型。
5.根據(jù)權(quán)利要求1所述的縱波和轉(zhuǎn)換波疊前聯(lián)合反演方法,其特征在于,所述轉(zhuǎn)換波時間域的轉(zhuǎn)換波合成記錄誤差的表達式為:
其中,為迭代次數(shù)為k時的彈性參數(shù)模型。
6.根據(jù)權(quán)利要求1所述的縱波和轉(zhuǎn)換波疊前聯(lián)合反演方法,其特征在于,所述縱波時間域的縱波和轉(zhuǎn)換波的總誤差rk的表達式為:
其中,λ為縱波和轉(zhuǎn)換波誤差調(diào)節(jié)參數(shù)。
該專利技術(shù)資料僅供研究查看技術(shù)是否侵權(quán)等信息,商用須獲得專利權(quán)人授權(quán)。該專利全部權(quán)利屬于中國石油化工股份有限公司;中國石油化工股份有限公司石油勘探開發(fā)研究院,未經(jīng)中國石油化工股份有限公司;中國石油化工股份有限公司石油勘探開發(fā)研究院許可,擅自商用是侵權(quán)行為。如果您想購買此專利、獲得商業(yè)授權(quán)和技術(shù)合作,請聯(lián)系【客服】
本文鏈接:http://www.17sss.com.cn/pat/books/201810857538.2/1.html,轉(zhuǎn)載請聲明來源鉆瓜專利網(wǎng)。
- 圖像轉(zhuǎn)換設(shè)備、圖像轉(zhuǎn)換電路及圖像轉(zhuǎn)換方法
- 數(shù)模轉(zhuǎn)換電路及轉(zhuǎn)換方法
- 轉(zhuǎn)換設(shè)備和轉(zhuǎn)換方法
- 占空比轉(zhuǎn)換電路及轉(zhuǎn)換方法
- 通信轉(zhuǎn)換方法、轉(zhuǎn)換裝置及轉(zhuǎn)換系統(tǒng)
- 模數(shù)轉(zhuǎn)換和模數(shù)轉(zhuǎn)換方法
- 轉(zhuǎn)換模塊以及轉(zhuǎn)換電路
- 熱電轉(zhuǎn)換材料、熱電轉(zhuǎn)換元件和熱電轉(zhuǎn)換模塊
- 熱電轉(zhuǎn)換材料、熱電轉(zhuǎn)換元件及熱電轉(zhuǎn)換模塊
- 熱電轉(zhuǎn)換材料、熱電轉(zhuǎn)換元件及熱電轉(zhuǎn)換模塊
- 中繼網(wǎng)絡(luò)中的聯(lián)合方法、系統(tǒng)、中繼站及移動臺
- 中繼網(wǎng)絡(luò)中的聯(lián)合方法、系統(tǒng)、中繼站及移動臺
- 便攜信息終端、以及便攜信息終端的控制方法
- 電法與地震同步聯(lián)合反演方法及系統(tǒng)
- 銀行聯(lián)合放款的放款方法、裝置及計算機可讀存儲介質(zhì)
- 一種基于自主學(xué)習(xí)群搜索算法的聯(lián)合調(diào)度方法及系統(tǒng)
- 一種基于多目標(biāo)種群群搜索算法的聯(lián)合調(diào)度方法及系統(tǒng)
- 聯(lián)合貸款中的罰息與利息結(jié)算方法、裝置、介質(zhì)和設(shè)備
- 一種基坑支護用H型鋼構(gòu)件
- 用于云環(huán)境的語義元聯(lián)合代理





