[發明專利]快速粘聲波和粘彈性全波場反演

在審

• 申請號: 201580032794.2
• 申請日: 2015-04-22
• 公開（公告）號: CN106662664A
• 公開（公告）日: 2017-05-10
• 發明（設計）人: H·德恩里;A·凱尼伍斯基
• 申請（專利權）人: 埃克森美孚上游研究公司
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 北京紀凱知識產權代理有限公司11245
• 代理人: 徐東升,趙蓉民

1.一種方法，其包括：

利用計算機獲得初始地球物理模型；

利用所述計算機使用粘聲波或粘彈性波方程對前向波場建模；

利用所述計算機使用伴隨粘聲波或伴隨粘彈性波方程對伴隨波場建模；

其中分別地，所述伴隨粘聲波波方程基于是壓力和記憶變量的函數的輔助變量，或所述伴隨粘彈性波方程基于應力和記憶變量的組合；

利用所述計算機基于所述前向波場的模型和所述伴隨波場的模型的組合來獲得成本函數的梯度；以及

利用所述計算機使用所述成本函數的所述梯度來更新所述初始地球物理模型并獲得更新的地球物理模型。

2.根據權利要求1所述的方法，其進一步包括使用所述輔助變量以使記憶變量的空間導數從所述伴隨粘聲波或粘彈性波方程中退出。

3.根據權利要求1或2所述的方法，其中獲得所述地球物理模型包括使用初始流變學模型和聲學或彈性介質參數來獲得針對給定品質因子的松弛參數。

4.根據權利要求1至3中的任一項所述的方法，其中使用所述伴隨粘聲波方程建模所述伴隨波場。

5.根據權利要求1至4中的任一項所述的方法，其中所述輔助變量是其中p是壓力，m_l是機構l的所述記憶變量，其中l是大于或等于1的整數，并且α_l是機構l的反演參數。

6.根據權利要求1至5中的任一項所述的方法，其中所述伴隨粘聲波方程是

$\frac{\∂{\stackrel{\‾}{\mathrm{\ζ}}}_l}{\∂t}+{\mathrm{\φ}}_1{\stackrel{\‾}{\mathrm{\ζ}}}_l-{\mathrm{\φ}}_1{\stackrel{\‾}{r}}_l=0,$

其中

是散度算子，

t是時間，

κ是非松弛體積模量(lim_ω→∞κ(ω)→κ)，

ρ是質量密度，

v是速度(在三維空間中(v＝{v_x v_y v_z}^T)，

是伴隨壓力，

是成本函數，其中θ是模型參數(即，Q或α_l、κ以及ρ的某組合)

是伴隨速度，以及

φ_l是機構l的松弛頻率，

和分別是目標函數關于所述壓力和速度的導數。

7.根據權利要求1至6中的任一項所述的方法，其中使用所述伴隨粘彈性方程建模所述伴隨波場。

8.根據權利要求7所述的方法，其中所述輔助變量是其中σ是應力，m_l是機構l的記憶變量，其中l是大于或等于1的整數，C是各向同性的非松弛系統的彈性本構關系并且D_l是記憶系統的本構關系。