[發明專利]各向同性彈性參數保幅反演方法及裝置

在審

• 申請號: 201811188186.2
• 申請日: 2018-10-12
• 公開（公告）號: CN109541682A
• 公開（公告）日: 2019-03-29
• 發明（設計）人: 毛偉建;歐陽威;詹毅;李合群
• 申請（專利權）人: 中國石油天然氣集團有限公司;中國石油集團東方地球物理勘探有限責任公司;中國科學院測量與地球物理研究所
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 北京三友知識產權代理有限公司 11127
• 代理人: 李輝

1.一種各向同性彈性參數保幅反演方法，其特征在于，包括：

獲取各向同性彈性介質模型中每個成像點的成像反演結果；所述成像反演結果包括P波散射場上的逆廣義拉東變換反投影算子；

基于所述逆廣義拉東變換反投影算子構建非線性反演系統；

確定所述各向同性彈性介質模型的照明矩陣和二次項透射系數；

根據所述逆廣義拉東變換反投影算子、所述照明矩陣和所述二次項透射系數，求解所述非線性反演系統，獲得所述各向同性彈性介質模型中每個成像點處的二次非線性保幅反演值。

2.如權利要求1所述的各向同性彈性參數保幅反演方法，其特征在于，所述各向同性彈性介質模型為三維模型；

相應的，三維空間下，所述非線性反演系統，包括：

其中，分別為三維模型中不同權重的逆廣義Radon變換反投影算子；U^P(z)為P波散射位移場；為照明矩陣；a₁₁(z)、a₁₂(z)、a₁₃(z)、a₂₁(z)、a₂₂(z)、a₂₃(z)、a₃₁(z)、a₃₂(z)、a₃₃(z)分別為照明矩陣中的元素；f₁(z)、f₂(z)、f₃(z)分別為三維模型中各向同性彈性參數的非線性二次組合函數，z為成像點。

3.如權利要求2所述的各向同性彈性參數保幅反演方法，其特征在于，所述二次項透射系數，包括：

其中，分別為三維模型中的二次項透射系數，x為散射點；ω₀為參考頻率；r為接收點；l(x)為與散射區域和散射點x有關的標量；s為震源；為符號函數；為關于x的梯度；ρ⁰(x)為三維背景模型中密度參數；e₃、e₂分別為三維中的單位向量，且e₃＝(0,0,1)，e₂＝(0,1,0)；為三維背景模型中P波速度；e為自然常數；i為虛數單位；ξ為三維空間中單位球面上的積分變量；λ⁰(x)和μ⁰(x)為三維背景模型中拉梅常數。

4.如權利要求1所述的各向同性彈性參數保幅反演方法，其特征在于，所述各向同性彈性介質模型為二維模型；

相應的，二維空間下，所述非線性反演系統，包括：

其中，分別為二維模型中不同權重的逆廣義Radon變換反投影算子；U^P(z)為P波散射位移場；為照明矩陣；分別為照明矩陣的元素，分別為二維模型中各向同性彈性參數的非線性二次組合函數，z為成像點。

5.如權利要求4所述的各向同性彈性參數保幅反演方法，其特征在于，所述二次項透射系數，包括：

其中，分別為二維模型中的二次項透射系數，θ為散射夾角，x為散射點；ω為頻率；r為接收點；l(x)為與散射區域和散射點x有關的標量；為符號函數；為關于x的梯度；ρ⁰(x)為二維背景模型中密度參數；e₃為二維中的單位向量，且e₃＝(0,1)；為二維背景模型中P波速度；e為自然常數；i為虛數單位；λ⁰(x)和μ⁰(x)為二維背景模型中拉梅常數。