[發(fā)明專利]光纖井中地震數(shù)據(jù)剖面的獲取方法及裝置

在審

• 申請?zhí)枺?/td>: 202011006244.2
• 申請日: 2020-09-23
• 公開（公告）號: CN112305616A
• 公開（公告）日: 2021-02-02
• 發(fā)明（設(shè)計(jì)）人: 蔡志東;李彥鵬;張建軍;徐剛;劉偉
• 申請（專利權(quán)）人: 中國石油天然氣集團(tuán)有限公司;中國石油集團(tuán)東方地球物理勘探有限責(zé)任公司
• 主分類號: G01V1/40
• 分類號: G01V1/40;G01V1/30;G01V1/22
• 代理公司: 北京三友知識產(chǎn)權(quán)代理有限公司 11127
• 代理人: 薛平;周曉飛

1.一種光纖井中地震數(shù)據(jù)剖面的獲取方法，其特征在于，包括：

獲取光纖井中地震數(shù)據(jù)；

根據(jù)光纖井中地震數(shù)據(jù)的下行波數(shù)據(jù)，計(jì)算反射系數(shù)；

根據(jù)光纖井中地震數(shù)據(jù)的上行波數(shù)據(jù)，獲得光纖井中地震上行波剖面；

根據(jù)所述反射系數(shù)和光纖井中地震上行波剖面，反演得到光纖井中地震數(shù)據(jù)剖面。

2.如權(quán)利要求1所述的方法，其特征在于，獲取光纖井中地震數(shù)據(jù)，包括：

獲取井中分布式光纖采集的零井源距井中地震數(shù)據(jù)。

3.如權(quán)利要求1所述的方法，其特征在于，獲取光纖井中地震數(shù)據(jù)之后，還包括：

對光纖井中地震數(shù)據(jù)進(jìn)行預(yù)處理，所述預(yù)處理保護(hù)光纖井中地震數(shù)據(jù)的原始振幅信息。

4.如權(quán)利要求3所述的方法，其特征在于，對光纖井中地震數(shù)據(jù)進(jìn)行預(yù)處理，包括：對光纖井中地震數(shù)據(jù)進(jìn)行光纜耦合噪聲壓制和/或隨機(jī)噪聲壓制處理。

5.如權(quán)利要求1所述的方法，其特征在于，根據(jù)光纖井中地震數(shù)據(jù)的下行波數(shù)據(jù)，計(jì)算反射系數(shù)，包括：

對光纖井中地震數(shù)據(jù)的下行波數(shù)據(jù)進(jìn)行初至拾取，得到隨深度變化的初至?xí)r間數(shù)列；

根據(jù)初至?xí)r間數(shù)列進(jìn)行振幅拾取，得到隨深度變化的振幅數(shù)列；

對振幅數(shù)列進(jìn)行冪函數(shù)擬合，得到深度域振幅衰減系數(shù)；

根據(jù)振幅數(shù)列和深度域振幅衰減系數(shù)，計(jì)算出標(biāo)準(zhǔn)振幅衰減數(shù)列；

根據(jù)振幅數(shù)列和標(biāo)準(zhǔn)振幅衰減數(shù)列，得到補(bǔ)償后振幅數(shù)據(jù)；

對補(bǔ)償后振幅數(shù)據(jù)中相鄰的振幅求取差值，得到深度振幅差關(guān)系數(shù)列；

利用深度振幅差關(guān)系數(shù)列，計(jì)算得到反射系數(shù)。

6.如權(quán)利要求5所述的方法，其特征在于，根據(jù)振幅數(shù)列和深度域振幅衰減系數(shù)，計(jì)算出標(biāo)準(zhǔn)振幅衰減數(shù)列，包括按如下公式計(jì)算出標(biāo)準(zhǔn)振幅衰減數(shù)列：

A′_h＝A₀h^k；

其中，A′_h為標(biāo)準(zhǔn)振幅衰減數(shù)列；A₀為振幅數(shù)列中零深度位置的振幅值；h為地層深度；k為深度域振幅衰減系數(shù)；

根據(jù)振幅數(shù)列和標(biāo)準(zhǔn)振幅衰減數(shù)列，得到補(bǔ)償后振幅數(shù)據(jù)，包括按如下公式計(jì)算得到補(bǔ)償后振幅數(shù)據(jù)：

A″_h＝A_h-A′_h＝A_h-A₀4π(h-h₀)²；

其中，A″_h為補(bǔ)償后振幅數(shù)據(jù)；A_h為振幅數(shù)列；h₀為初始地層深度；

對補(bǔ)償后振幅數(shù)據(jù)中相鄰的振幅求取差值，得到深度振幅差關(guān)系數(shù)列，包括按如下公式計(jì)算得到深度振幅差關(guān)系數(shù)列：

ΔA″_h＝A″_h-A″_h+Δh；

其中，ΔA″_h為深度振幅差關(guān)系數(shù)列；Δh為深度增量；A″_h+Δh為相鄰深度的補(bǔ)償后振幅數(shù)據(jù)；

利用深度振幅差關(guān)系數(shù)列，計(jì)算得到反射系數(shù)，包括按如下公式計(jì)算得到反射系數(shù)數(shù)列：

R_h＝(A″_h-ΔA″_h)/A″_h；

其中，R_h為反射系數(shù)數(shù)列。

7.如權(quán)利要求1所述的方法，其特征在于，根據(jù)光纖井中地震數(shù)據(jù)的上行波數(shù)據(jù)，獲得光纖井中地震上行波剖面，包括：

對光纖井中地震數(shù)據(jù)的上行波數(shù)據(jù)，依次進(jìn)行振幅補(bǔ)償處理、反褶積處理、波場分離處理和動(dòng)校正處理，獲得光纖井中地震上行波剖面。