[發(fā)明專利]一種三維無射線追蹤回折波層析成像方法及裝置

有效

• 申請?zhí)枺?/td>: 201410431526.5
• 申請日: 2014-08-28
• 公開（公告）號: CN104181593B
• 公開（公告）日: 2017-01-11
• 發(fā)明（設計）人: 熊登;錢忠平;王成祥;程春華;王仕儉
• 申請（專利權）人: 中國石油天然氣集團公司;中國石油集團東方地球物理勘探有限責任公司
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 北京三友知識產權代理有限公司11127
• 代理人: 湯在彥

1.一種三維無射線追蹤回折波層析成像方法，其特征在于，該方法包括：

拾取單炮初至，對拾取的每個單炮初至按照方位角度間隔選定多個方位角的初至，并對每一個方位角初至按照偏移距遞增方式排序，形成0～x范圍內的初至曲線；其中，x為最大偏移距；

以偏移距為橫坐標、初至時間為縱坐標，對0～x范圍內的初至曲線進行低階多項式最小二乘法整體擬合，得到每一單炮初至的方位角對應地多項式系數；

根據所述多項式系數獲得對應地每一單炮初至的方位角慢度曲線；

根據所述每一單炮初至的方位角慢度曲線獲得拾取的每個單炮初至所選定的方位角對應地回折波回折點的速度和深度；

根據拾取的每個單炮初至所選定的方位角對應地回折波回折點的速度和深度插值成三維速度場。

2.如權利要求1所述的方法，其特征在于，所述拾取的每個單炮初至所選定的方位角對應地回折波回折點的速度獲取方法為：

對拾取的每一單炮初至的方位角慢度曲線在0～x范圍內任取M點為回折波回折點；

求出點i對應位置的慢度曲線值的倒數，得到點i的速度；其中，i＝1,2，…，M。

3.如權利要求2所述的方法，其特征在于，所述拾取的每個單炮初至所選定的方位角對應地回折波回折點的深度獲取方法為：

對點i的深度采用積分公式獲得；其中，所述積分公式為：

zi=1π∫0xicosh-1(dtdx(dtdx)i)dx]]>

其中，z_i為回折點i對應的深度，x_i為點i對應的偏移距。

4.如權利要求1所述的方法，其特征在于，所述每一單炮初至的方位角慢度曲線獲取的方法為：

根據所述多項式系數獲得多項式；

對所述多項式進行求導，獲得每一單炮初至的方位角慢度曲線。

5.如權利要求1所述的方法，其特征在于，所述方位角度間隔為相當或不相等的方位角度間隔。

6.一種三維無射線追蹤回折波層析成像裝置，其特征在于，該裝置包括：

初至曲線獲取單元，用于拾取單炮初至，對拾取的每個單炮初至按照方位角度間隔選定至少兩個方位角的初至，并對每一個方位角初至按照偏移距遞增方式排序，形成0～x范圍內的初至曲線；其中，x為最大偏移距；

擬合單元，用于以偏移距為橫坐標、初至時間為縱坐標，對0～x范圍內的初至曲線進行低階多項式最小二乘法整體擬合，得到每一單炮初至的方位角對應地多項式系數；

方位角慢度曲線獲取單元，用于根據所述多項式系數獲得對應地每一單炮初至的方位角慢度曲線；

回折波回折點的速度和深度獲取單元，用于根據所述每一單炮初至的方位角慢度曲線獲得拾取的每個單炮初至所選定的方位角對應地回折波回折點的速度和深度；

插值單元，用于根據拾取的每個單炮初至所選定的方位角對應地回折波回折點的速度和深度插值成三維速度場。

7.如權利要求6所述的裝置，其特征在于，所述回折波回折點的速度和深度獲取單元包括：回折波回折點獲取模塊和第一計算模塊；其中，

所述回折波回折點獲取模塊，用于對拾取的每一單炮初至的方位角慢度曲線在0～x范圍內任取M點為回折波回折點；

所述第一計算模塊，用于求出點i對應位置的慢度曲線值的倒數，得到點i的速度；其中，i＝1,2，…，M。