[發(fā)明專(zhuān)利]最大炮檢距確定方法

在審

• 申請(qǐng)?zhí)枺?/td>: 201310594960.0
• 申請(qǐng)日: 2013-11-21
• 公開(kāi)（公告）號(hào): CN103645501A
• 公開(kāi)（公告）日: 2014-03-19
• 發(fā)明（設(shè)計(jì)）人: 胡善政;李亞林;王勤耕;何光明;陳愛(ài)萍;敬龍江;陳燕雄;童濤;黎書(shū)琴;袁光銀;楊智超;王曉陽(yáng)
• 申請(qǐng)（專(zhuān)利權(quán)）人: 中國(guó)石油集團(tuán)川慶鉆探工程有限公司地球物理勘探公司
• 主分類(lèi)號(hào): G01V1/28
• 分類(lèi)號(hào): G01V1/28
• 代理公司: 北京銘碩知識(shí)產(chǎn)權(quán)代理有限公司 11286
• 代理人: 王兆賡;譚昌馳

1.一種最大炮檢距確定方法，其特征在于，包括以下步驟：

（a）從地質(zhì)模型中選取目標(biāo)反射層；

（b）在地質(zhì)模型的整個(gè)表面上均勻地布置多個(gè)勘探炮和多個(gè)檢波器；

（c）使所有勘探炮分別產(chǎn)生地震波，并且當(dāng)多個(gè)勘探炮中的一個(gè)勘探炮產(chǎn)生地震波時(shí)，多個(gè)檢波器分別接收該地震波被目標(biāo)反射層反射回來(lái)的反射波，從而獲得該地震波在各個(gè)檢波器中的炮檢對(duì)能量值；

（d）在所有勘探炮分別產(chǎn)生地震波之后，基于多個(gè)檢波器中的各個(gè)檢波器所獲得的炮檢對(duì)能量值，繪制炮檢對(duì)能量值圖；

（e）基于炮檢對(duì)能量值圖，確定針對(duì)當(dāng)前目標(biāo)反射層的最大炮檢距，

其中，炮檢對(duì)能量值是指，當(dāng)產(chǎn)生一次地震波時(shí)一個(gè)檢波器所接收到的反射波的能量的總和。

2.根據(jù)權(quán)利要求1所述的最大炮檢距確定方法，其特征在于，在步驟（c）中，利用有限差分算子，對(duì)當(dāng)前產(chǎn)生的地震波進(jìn)行正向延拓到目標(biāo)反射層，然后再?gòu)哪繕?biāo)反射層反向延拓到檢波器，來(lái)模擬地震波被目標(biāo)反射層反射回來(lái)的反射波，從而每個(gè)檢波器計(jì)算接收到的反射波的能量的總和。

3.根據(jù)權(quán)利要求1所述的最大炮檢距確定方法，其特征在于，在步驟（c）中，當(dāng)產(chǎn)生地震波時(shí)，地震波在頻率域中的傳播過(guò)程以下面的公式（1）來(lái)表示，并且從公式（1）中計(jì)算出地震波的波長(zhǎng)值U(r_s,r,ω)，

（k²+Δ）U(r_s,r,ω)=-f(ω)δ(r-r_s)????公式（1）

其中，k=ω/c₀(r)為波數(shù)，ω表示地震波的頻率，c₀(r)為位于目標(biāo)反射層的點(diǎn)r=(x,y,z)處的背景場(chǎng)速度值，Δ為拉普拉斯算子，f(ω)為勘探炮函數(shù)，δ(r-r_s)是單位脈沖函數(shù)，r_s表示勘探炮的三維坐標(biāo)。

4.根據(jù)權(quán)利要求3所述的最大炮檢距確定方法，其特征在于，在步驟（c）中，利用有限差分算子，對(duì)當(dāng)前產(chǎn)生的地震波進(jìn)行正向延拓到目標(biāo)反射層，然后再?gòu)哪繕?biāo)反射層反向延拓到檢波器，

其中，檢波器接收到的反向延拓的反射波的波場(chǎng)通過(guò)公式（2）來(lái)表示，

P(r_s,r_g,ω)=

2k²∫_V(r)G(r,r_s;ω)m(r)G(r,r_g;ω)d³r????公式（2）

其中，G(r,r_s;ω)和G(r,r_g;ω)分別是頻率域中勘探炮和檢波器位置的格林函數(shù)，r_g表示檢波器的三維坐標(biāo)，m(r)是位于目標(biāo)反射層的點(diǎn)r=(x,y,z)處速度模型擾動(dòng)，V(r)為目標(biāo)反射層的所有點(diǎn)的集合。

5.根據(jù)權(quán)利要求4所述的最大炮檢距確定方法，其特征在于，多個(gè)檢波器中的一個(gè)檢波器針對(duì)當(dāng)前產(chǎn)生地震波的勘探炮的炮檢對(duì)能量值通過(guò)公式（3）來(lái)表示，

P_sr(r_s,r_g,ω)=|P(r_s,r_g,ω)|????公式（3）。