地震波的传播速度是地震勘探中重要参数之一。研究表明, 地下介质具有广泛的各向异性特征^[1], 而各向异性介质中地震波的传播速度受很多因素的影响。随着地震勘探的不断发展以及勘探目标越来越复杂, 各向同性的假设条件已经不能满足勘探精度的需求^[2]。由于地震波在各向同性介质和各向异性介质中传播时存在差异, 所以地震数据的处理、解释和反演等方面需要考虑地下介质的各向异性, 否则可能会导致错误的处理和解释结果^[3-4]。因此, 进行各向异性介质中地震波传播速度的研究具有重要的实际意义^[5]。很多国内外学者对各向异性介质中地震波传播特征、速度分析、数值模拟、参数提取和属性分析等方面进行了研究^[6-18]。

目前油气勘探的主要勘探目标仍然以沉积岩为主。沉积岩石中各向异性是指在横向各向同性而在纵向上为非均质性, 一般称之为横向各向同性(TI)介质, 其对称轴垂直于水平面时称为VTI介质, 对称轴平行于水平面时称为HTI介质。而在页岩油气勘探中, 页岩的各向异性特征也可以用TI介质来刻画^[19-24], 但是页岩一般具有中强各向异性。在褶皱和逆冲作用下, 部分地层发生倾斜时, TI介质的对称轴不再水平或者垂直。当VTI/HTI介质的本构坐标系和观测坐标系不重合时, VTI/HTI介质在观测系统下转变为TTI介质^[25]。关于TTI介质的波场传播、偏移成像、速度分析、AVO等方面的研究成果较多^[26-35]。但利用Christoffel方程得到的TTI介质的相速度理论精确公式比较复杂, 在实际应用中推广比较困难。根据Thomsen理论, 地球介质中绝大部分的地层可以近似为弱各向异性介质^[36]。梁锴等^[25]、吴国忱等^[37-39]基于弱各向异性假设提出了精确相速度的近似公式。基于弱各向异性假设的近似公式在各向异性参数值比较小的情况下, 可以得到较好的结果。然而, 随着勘探目标的不断复杂化, 加上页岩油气勘探中页岩的各向异性程度较大, 因而弱各向异性近似公式的精度不能满足勘探需要, 甚至不再适用。

本文以TTI介质弹性波相速度的精确表达式为出发点, 经过Thomsen参数表征后, 利用近似的配方法将精确公式内的D项近似表示为两项的平方和, 从而实现多项式的开方运算, 进一步推导了TTI介质弹性波相速度的一般性近似表达式, 并与精确值、弱各向异性近似值进行了对比。一般性近似表达式比精确公式简洁, 在一定范围内与精确解吻合较好。

1 TTI介质弹性波相速度 1.1 精确解

TTI介质的刚度矩阵(C_TTI)可以由VTI介质刚度矩阵(C_VTI)经Bond变换^[40]得到, 即:

$ {\mathit{\boldsymbol{C}}_{{\rm{TTI}}}} = {\mathit{\boldsymbol{M}}_{{\theta _R}}}{\mathit{\boldsymbol{C}}_{{\rm{VTI}}}}\mathit{\boldsymbol{M}}_{{\theta _R}}^{\rm{T}} $

(1)

式中: M_θR表示Bond变换矩阵^[37]; θ_R表示TTI介质对称轴的极化角。

根据弹性力学的本构方程、几何方程、运动微分方程和TTI介质的刚度矩阵(C_TTI), 并且忽略体力项, 建立TTI介质弹性波波动方程:

$ \rho \frac{{{\partial ^2}}}{{\partial {t^2}}}\mathit{\boldsymbol{U}} = \mathit{\boldsymbol{L}}\left( {{\mathit{\boldsymbol{C}}_{{\rm{TTI}}}}{\mathit{\boldsymbol{L}}^{\rm{T}}}\mathit{\boldsymbol{U}}} \right) $

(2)

式中: U=(u_x, u_y, u_z)^T为位移矢量; ρ为介质密度; t为时间变量; L为偏导数矩阵^[37]。

假设平面波解为:

$ \mathit{\boldsymbol{U}} = \mathit{\boldsymbol{P}}\exp \left[ {{\rm{i}}k\left( {\mathit{\boldsymbol{n}} \cdot \mathit{\boldsymbol{x}} - vt} \right)} \right] $

(3)

式中: x=(x, y, z)^T为位移矢量; n=(n_x, n_y, n_z)^T= (sinθcosφ, sinθsinφ, cosθ)为平面波传播方向, 其中θ为传播方向的极化角, φ为传播方向的方位角; v为平面波传播速度(相速度); k=ω/v为波数, 其中ω为角频率; P=(p_x, p_y, p_z)^T为平面波偏振方向(在波前满足n·x－vt为常数)。

将平面波解((3)式)代入TTI介质弹性波波动方程((2)式), 得到TTI介质Kelvin-Christoffel方程^[37-39]。Kelvin-Christoffel方程存在非零解的条件是其系数矩阵行列式为零, 求解得到行列式方程的3个根, 即为TTI介质的P波, SV波, SH波相速度的精确表达式:

$ \left\{ \begin{array}{l} {v_{\rm{P}}} = \sqrt {\frac{1}{{2\rho }}\left[ {\left( {{C_{11}} + {C_{44}}} \right)A + \left( {{C_{33}} + {C_{44}}} \right)B + \sqrt D } \right]} \\ {v_{{\rm{SV}}}} = \sqrt {\frac{1}{{2\rho }}\left[ {\left( {{C_{11}} + {C_{44}}} \right)A + \left( {{C_{33}} + {C_{44}}} \right)B - \sqrt D } \right]} \\ {v_{{\rm{SH}}}} = \sqrt {\frac{1}{\rho }\left[ {{C_{66}}A + {C_{44}}B} \right]} \end{array} \right. $

(4)

式中:A=sin²θsin²φ+(sinθcosφcosθ_R+cosθsinθ_R)²; B=(cosθcosθ_R－sinθcosφsinθ_R)²; D=[(C₁₁－C₄₄)A－(C₃₃－C₄₄)B]²+4(C₁₃+C₄₄)²AB; C_ij表示VTI介质的刚度矩阵元素。

1.2 弱各向异性近似

由刚度矩阵元素表示的相速度公式较为复杂, 各向异性的度量不直观。THOMSEN提出了一套表征TI介质弹性性质的参数, 称为Thomsen参数^[36]。因此, VTI介质刚度矩阵的Thomsen参数表征为:

$ {C_{11}} = {C_{22}} = \left( {1 + 2\varepsilon } \right)\rho v_{{{\rm{P0}}}}^2 $

(5a)

$ {C_{12}} = \rho v_{{{\rm{P0}}}}^2\left[ {1 + 2\varepsilon - 2\left( {1 - f} \right)\left( {1 + 2\gamma } \right)} \right] $

(5b)

$ {C_{13}} = {C_{23}} = \rho v_{{{\rm{P0}}}}^2\sqrt {f\left( {f + 2\delta } \right)} - \rho v_{{{\rm{S0}}}}^2 $

(5c)

$ {C_{33}} = \rho v_{{{\rm{P0}}}}^2 $

(5d)

$ {C_{44}} = {C_{55}} = \rho v_{{{\rm{S0}}}}^2 $

(5e)

$ {C_{66}} = \left( {1 + 2\gamma } \right)\rho v_{{{\rm{S0}}}}^2 $

(5f)

$ f = 1 - \frac{{v_{{{\rm{S0}}}}^2}}{{v_{{{\rm{P0}}}}^2}} = 1 - \frac{{{C_{55}}}}{{{C_{33}}}} $

(5g)

式中:v_P0为沿TI介质对称轴方向的纵波速度; v_S0为沿TI介质对称轴方向的横波速度; ε为度量准纵波各向异性强度的参数; γ为度量准横波各向异性强度的参数; δ为连接v_P和v_S的一种过渡性参数; f为用于简化表示的过渡性参数。

将公式(5)代入TTI介质弹性波相速度表达式((4)式)中, 并令φ=0, 则可以得到XOZ平面内TTI介质弹性波相速度的Thomsen表征, 即:

$ \left\{ \begin{array}{l} {v_{\rm{P}}} = {v_{{{\rm{P0}}}}}\sqrt {1 - \frac{f}{2} + \varepsilon {{\sin }^2}\left( {\theta + {\theta _R}} \right) + \frac{1}{{2v_{{{\rm{P0}}}}^2}}\sqrt D } \\ {v_{{\rm{SV}}}} = {v_{{{\rm{P0}}}}}\sqrt {1 - \frac{f}{2} + \varepsilon {{\sin }^2}\left( {\theta + {\theta _R}} \right) - \frac{1}{{2v_{{{\rm{P0}}}}^2}}\sqrt D } \\ {v_{{\rm{SH}}}} = {v_{{{\rm{S0}}}}}\sqrt {1 + 2\gamma {{\sin }^2}\left( {\theta + {\theta _R}} \right)} \end{array} \right. $

(6)

式中:D=v_P0⁴f²{1+4(ε²/f²)sin⁴(θ+θ_R)+(4/f)· sin²(θ+θ_R){2δcos²(θ+θ_R)+ε[sin²(θ+θ_R)－cos²(θ+θ_R)]}}。

在(6)式中, 根据弱各向异性近似方法(即假设ε和δ充分小)可以忽略其高阶项(2次及以上), 且因ε和δ充分小, $ \sqrt {1 + x} \approx 1 + 1/2x$成立, 由此得到;

$ \begin{array}{*{20}{c}} {\sqrt D \approx v_{{{\rm{P0}}}}^2f\left\{ {1 + \frac{2}{f}{{\sin }^2}\left( {\theta + {\theta _R}} \right)\left\{ {2\delta {{\cos }^2}\left( {\theta + {\theta _R}} \right) + } \right.} \right.}\\ {\left. {\left. {\varepsilon \left[ {{{\sin }^2}\left( {\theta + {\theta _R}} \right) - {{\cos }^2}\left( {\theta + {\theta _R}} \right)} \right]} \right\}} \right\}} \end{array} $

(7)

将(7)式代入(6)式, 得到TTI介质Thomsen弱各向异性近似弹性波相速度:

$ \left\{ \begin{array}{l} {v_{\rm{P}}} \approx {v_{{{\rm{P0}}}}} \cdot \\ \sqrt {1 + 2\delta {{\sin }^2}\left( {\theta + {\theta _R}} \right){{\cos }^2}\left( {\theta + {\theta _R}} \right) + 2\varepsilon {{\sin }^4}\left( {\theta + {\theta _R}} \right)} \\ {v_{{\rm{SV}}}} \approx {v_{{{\rm{S0}}}}} \cdot \\ \sqrt {1 + 2\frac{{v_{{{\rm{P0}}}}^2}}{{v_{{{\rm{S0}}}}^2}}\left( {\varepsilon - \delta } \right){{\sin }^2}\left( {\theta + {\theta _R}} \right){{\cos }^2}\left( {\theta + {\theta _R}} \right)} \\ {v_{{\rm{SH}}}} = {v_{{\rm{S0}}}}\sqrt {1 + 2\gamma {{\sin }^2}\left( {\theta + {\theta _R}} \right)} \end{array} \right. $

(8)

对(8)式再次应用$ \sqrt {1 + x} \approx 1 + 1/2x $, 整理后得到:

$ \begin{array}{l} {v_{\rm{P}}} \approx {v_{{{\rm{P0}}}}}\left[ {1 + \delta {{\sin }^2}\left( {\theta + {\theta _R}} \right){{\cos }^2}\left( {\theta + {\theta _R}} \right) + } \right.\\ \;\;\;\;\;\;\;\;\;\left. {\varepsilon {{\sin }^4}\left( {\theta + {\theta _R}} \right)} \right] \end{array} $

(9a)

$ \begin{array}{l} {v_{{\rm{SV}}}} \approx {v_{{{\rm{S0}}}}}\left[ {1 + \frac{{v_{{{\rm{P0}}}}^2}}{{v_{{{\rm{S0}}}}^2}}\left( {\varepsilon - \delta } \right){{\sin }^2}\left( {\theta + {\theta _R}} \right) \cdot } \right.\\ \;\;\;\;\;\;\;\;\left. {{{\cos }^2}\left( {\theta + {\theta _R}} \right)} \right] \end{array} $

(9b)

$ {v_{{\rm{SH}}}} \approx {v_{{{\rm{S0}}}}}\left[ {1 + \gamma {{\sin }^2}\left( {\theta + {\theta _R}} \right)} \right] $

(9c)

(9) 式是TTI介质相速度的弱各向异性近似公式的另一种表征^[37-39], 该形式下弹性波相速度是各向异性参数的线性函数。

1.3 一般性近似

弱各向异性近似的条件是假设各向异性参数数值较小, 各向异性程度较弱, 因此对于页岩等中强各向异性介质, 采用弱各向异性近似计算其相速度, 其精度可能无法满足勘探需要。为了得到更为一般性的近似表达式, 在不限制各向异性参数大小的前提下, 采取近似的配方法对(4)式中的$ \sqrt D $进行近似, 即$ \sqrt D \approx \sqrt {{{\left( {a + b} \right)}^2}} = a + b $, 从而消除相速度精确表达式中的一重根号。具体过程如下:

将公式(5)代入TTI介质弹性波相速度表达式((4)式)的D中, 并令φ=0, 则有:

$ \begin{array}{l} D = {\left[ {\left( {1 + 2\varepsilon } \right)\rho v_{{{\rm{P0}}}}^2 - \rho v_{{{\rm{S0}}}}^2} \right]^2}{\sin ^4}\left( {\theta + {\theta _R}} \right) + \left( {\rho v_{{{\rm{P0}}}}^2 - } \right.\\ \;\;\;\;\;\;{\left. {\rho v_{{{\rm{S0}}}}^2} \right)^2}{\cos ^4}\left( {\theta + {\theta _R}} \right) - 2\left[ {\left( {1 + 2\varepsilon } \right)\rho v_{{{\rm{P0}}}}^2 - \rho v_{{{\rm{S0}}}}^2} \right] \cdot \\ \;\;\;\;\;\;\left( {\rho v_{{{\rm{P0}}}}^2 - \rho v_{{{\rm{S0}}}}^2} \right){\sin ^2}\left( {\theta + {\theta _R}} \right){\cos ^2}\left( {\theta + {\theta _R}} \right) + \\ \;\;\;\;\;\;4\left[ {\left( {1 + 2\delta } \right)\rho v_{{{\rm{P0}}}}^2 - \rho v_{{{\rm{S0}}}}^2} \right]\left( {\rho v_{{{\rm{P0}}}}^2 - \rho v_{{{\rm{S0}}}}^2} \right){\sin ^2}\left( {\theta + } \right.\\ \;\;\;\;\;\;\left. {{\theta _R}} \right){\cos ^2}\left( {\theta + {\theta _R}} \right) = \left\{ {\left[ {\left( {1 + 2\varepsilon } \right)\rho v_{{{\rm{P0}}}}^2 - \rho v_{{{\rm{S0}}}}^2} \right] \cdot } \right.\\ \;\;\;\;\;\;{\left. {{{\sin }^2}\left( {\theta + {\theta _R}} \right) + \left( {\rho v_{{{\rm{P0}}}}^2 - \rho v_{{{\rm{S0}}}}^2} \right){{\cos }^2}\left( {\theta + {\theta _R}} \right)} \right\}^2} + \\ \;\;\;\;\;\;4{\rho ^2}v_{{{\rm{P0}}}}^2\left( {\delta - \varepsilon } \right) \cdot 2\left( {1 - \frac{{v_{{{\rm{P0}}}}^2}}{{v_{{{\rm{S0}}}}^2}}} \right){\sin ^2}\left( {\theta + {\theta _R}} \right) \cdot \\ \;\;\;\;\;\;{\cos ^2}\left( {\theta + {\theta _R}} \right) \end{array} $

(10)

在(10)式中, 对4ρ²v_P0⁴(δ－ε)·2(1－v_S0²/v_P0²)sin²(θ+θ_R)cos²(θ+θ_R)进行以下近似(具体近似过程见附录A), 即:

$ \begin{array}{l} 4{\rho ^2}v_{{{\rm{P0}}}}^2\left( {\delta - \varepsilon } \right) \cdot 2\left( {1 - \frac{{v_{{{\rm{S0}}}}^2}}{{v_{{{\rm{P0}}}}^2}}} \right){\sin ^2}\left( {\theta + {\theta _R}} \right){\cos ^2}\left( {\theta + } \right.\\ \left. {{\theta _R}} \right) \approx 4{\rho ^2}v_{{{\rm{P0}}}}^2\left( {\delta - \varepsilon } \right){\sin ^2}\left( {\theta + {\theta _R}} \right){\cos ^2}\left( {\theta + {\theta _R}} \right) \cdot \\ \left( {1 - \frac{{v_{{{\rm{P0}}}}^2}}{{v_{{{\rm{S0}}}}^2}}} \right) + 4{\rho ^2}v_{{{\rm{P0}}}}^2\left( {\delta - \varepsilon } \right){\sin ^2}\left( {\theta + {\theta _R}} \right){\cos ^2}\left( {\theta + } \right.\\ \left. {{\theta _R}} \right)\left[ {2\varepsilon {{\sin }^2}\left( {\theta + {\theta _R}} \right) + \left( {\delta - \varepsilon } \right){{\sin }^2}\left( {\theta + {\theta _R}} \right) \cdot } \right.\\ \left. {{{\cos }^2}\left( {\theta + {\theta _R}} \right)} \right] \end{array} $

(11)

将(11)式代入(10)式, 得到:

$ \begin{array}{l} D \approx \left\{ {\left[ {\left( {1 + 2\varepsilon } \right)\rho v_{{{\rm{P0}}}}^2 - \rho v_{{{\rm{S0}}}}^2} \right]{{\sin }^2}\left( {\theta + {\theta _R}} \right) + } \right.\\ \;\;\;\;{\left. {\left( {\rho v_{{{\rm{P0}}}}^2 - \rho v_{{{\rm{S0}}}}^2} \right){{\cos }^2}\left( {\theta + {\theta _R}} \right)} \right\}^2} + 4{\rho ^2}v_{{{\rm{P0}}}}^4\left( {\delta - } \right.\\ \;\;\;\;\left. \varepsilon \right){\sin ^2}\left( {\theta + {\theta _R}} \right){\cos ^2}\left( {\theta + {\theta _R}} \right) \cdot \left[ {1 - \frac{{v_{{{\rm{S0}}}}^2}}{{v_{{{\rm{P0}}}}^2}} + 2\varepsilon {{\sin }^2} \cdot } \right.\\ \;\;\;\;\left. {\left( {\theta + {\theta _R}} \right) + \left( {\delta - \varepsilon } \right){{\sin }^2}\left( {\theta + {\theta _R}} \right){{\cos }^2}\left( {\theta + {\theta _R}} \right)} \right]\\ \;\;\;\; = {a^2} + b\left( {2a + b} \right)\\ \;\;\;\; = {\left( {a + b} \right)^2} \end{array} $

(12)

式中:a=(ρv_P0²－ρv_S0²)cos²(θ+θ_R)+[(1+2ε)ρv_P0²－ρv_S0²]sin²(θ+θ_R); b=2ρv_P0²(δ－ε)sin²(θ+θ_R)cos²(θ+θ_R)。

将(12)式代入(6)式, 整理后得:

$ \left\{ \begin{array}{l} {v_{\rm{P}}} = {v_{{{\rm{P0}}}}}\left[ {1 + 2\varepsilon {{\sin }^2}\left( {\theta + {\theta _R}} \right) + \left( {\delta - \varepsilon } \right) \cdot } \right.\\ \;\;\;\;\;\;\;{\left. {{{\sin }^2}\left( {\theta + {\theta _R}} \right){{\cos }^2}\left( {\theta + {\theta _R}} \right)} \right]^{\frac{1}{2}}}\\ {v_{{\rm{SV}}}} \approx {v_{{{\rm{S0}}}}}\left[ {1 + \frac{{v_{{{\rm{P0}}}}^2}}{{v_{{{\rm{S0}}}}^2}}\left( {\varepsilon - \delta } \right){{\sin }^2}\left( {\theta + {\theta _R}} \right) \cdot } \right.\\ \;\;\;\;\;\;\;\;{\left. {{{\cos }^2}\left( {\theta + {\theta _R}} \right)} \right]^{\frac{1}{2}}} \end{array} \right. $

(13)

(13) 式即为基于近似的配方法推导的TTI介质弹性波相速度的一般性近似表达式。

2 数值算例 2.1 近似项数值测试

为了考察一般性近似((11)式)中关键近似项的近似程度, 对精确项:

$ 4{\rho ^2}v_{{{\rm{P0}}}}^2\left( {\delta - \varepsilon } \right){\sin ^2}\left( {\theta + {\theta _R}} \right){\cos ^2}\left( {\theta + {\theta _R}} \right) \cdot \left( {1 - \frac{{v_{{{\rm{S0}}}}^2}}{{v_{{{\rm{P0}}}}^2}}} \right) $

(14)

和其近似项:

$ \begin{array}{*{20}{c}} {4{\rho ^2}v_{{{\rm{P0}}}}^2\left( {\delta - \varepsilon } \right){{\sin }^2}\left( {\theta + {\theta _R}} \right){{\cos }^2}\left( {\theta + {\theta _R}} \right)\left[ {2\varepsilon {{\sin }^2}\left( {\theta + } \right.} \right.}\\ {\left. {\left. {{\theta _R}} \right) + \left( {\delta - \varepsilon } \right){{\sin }^2}\left( {\theta + {\theta _R}} \right){{\cos }^2}\left( {\theta + {\theta _R}} \right)} \right]} \end{array} $

(15)

进行了数值测试。所用模型为一组均匀TTI介质模型, 模型中相同参数为v_P0=3 300 m/s, v_S0=1 905 m/s, ρ=2.440×10³ kg/m³, θ_R=45°, φ=0, 不同的各向异性参数ε和δ如图 1所示。

由图 1可以看到, 近似项在一定精度范围内与精确值吻合较好。当ε=δ, 即椭圆各向异性介质时, 近似值与精确值相等。如果θ_R发生变化, 则精确值与近似值曲线相应地往右(或左)整体移动, 对近似程度不产生影响。

2.2 相速度数值测试

为了比较两种近似的效果, 针对上述均匀TTI介质模型, 利用TTI介质弹性波相速度的精确表达式((4)式)、弱各向异性近似表达式((8)式)和一般性近似表达式((13)式)计算qP波和qSV波的相速度, 结果如图 2和图 3所示。

图 2对比了不同TTI介质模型qP波的精确值、弱各向异性近似值和一般性各向异性近似值。由图 2可以看出, 在大多数角度范围内, 无论是弱各向异性情况还是中强各向异性情况, 一般性各向异性近似值比弱各向异性近似值更接近精确值。弱各向异性近似和一般性各向异性近似均受各向异性参数ε和δ的影响。当ε和δ相等时, 公式(4)、公式(8)和公式(13)相等, 所以3条速度曲线均重叠。近似效果还受|ε-δ|大小的影响, 当ε和δ的差值较小时, 一般性近似的近似程度比ε和δ的差值较大时好。

一般性近似的近似效果还受传播方向(θ)的影响, 并且沿对称轴方向的近似程度等于垂直对称轴方向的近似程度。

图 3对比了qSV波的精确值、弱各向异性近似值和一般性各向异性近似值。由图 3可以看出, 与qP波类似, 对于TTI介质, qSV波相速度的一般性各向异性近似效果较好, 大部分情况下, 一般性各向异性近似值要比弱各向异性近似值更接近精确值, 并且对中强各向异性参数也能有较好的近似效果。对比图 2和图 3可以发现, qP波速度的近似效果更优于qSV波速度。另外, 在图 3中, 当ε=0.1, δ=0.9时, 因弱各向异性近似表达式的根号中出现负值, 不存在弱各向异性近似值, 而一般性各向异性近似值是存在的, 这也从另一个角度说明了一般性各向异性近似表达式的适用范围更广。

3 结论

本文从TTI介质弹性波相速度精确表达式出发, 根据近似的配方法将表达式中的D近似地表示成两项的平方和, 从而消去一重根号, 推导出TTI介质qP波和qSV波相速度的一般性近似。通过相速度精确表达式、弱各向异性近似和一般性近似的对比, 得到如下认识和结论:

1) 一般性近似由于在推导过程中没有直接利用弱各向异性的假设, 因此不仅适用于弱各向异性, 也适用于中强各向异性的情况, 即在较大范围内, 一般性近似的精度高于弱各向异性近似的精度;

2) 对于同一种TTI介质, 由于qSV波相速度中各向异性参数有v_P0²/v_S0²的系数起放大作用, 所以qP波的一般性近似精度高于qSV波的精度。对于弱各向异性近似也有类似结论;

3) 与弱各向异性近似类似, TTI介质弹性波相速度一般性近似表达式仍然是v_P0, v_S, ε, δ, θ, θ_R的函数, 并受这些参数控制。

本文推导的TTI介质弹性波相速度一般性近似表达式, 虽然比精确公式简单, 但是仍然含有开方项, 具有一定的复杂性, 下一步研究可以考虑如何进一步消除开方项从而得到更为简洁的表征形式。

附录A 一般性近似公式的推导过程

Thomsen参数表征的TTI介质弹性波相速度的精确表达式为:

$ {v_{\rm{P}}} = \sqrt {\frac{1}{{2\rho }}\left\{ {\left[ {\left( {1 + 2\varepsilon } \right)\rho v_{{\rm{P0}}}^2 + \rho v_{{\rm{S0}}}^2} \right]{{\sin }^2}\left( {\theta + {\theta _R}} \right) + \left( {\rho v_{{\rm{P0}}}^2 + \rho v_{{\rm{S0}}}^2} \right){{\cos }^2}\left( {\theta + {\theta _R}} \right) + \sqrt D } \right\}} $

(A1)

$ {v_{{\rm{SV}}}} = \sqrt {\frac{1}{{2\rho }}\left\{ {\left[ {\left( {1 + 2\varepsilon } \right)\rho v_{{\rm{P0}}}^2 + \rho v_{{\rm{S0}}}^2} \right]{{\sin }^2}\left( {\theta + {\theta _R}} \right) + \left( {\rho v_{{\rm{P0}}}^2 + \rho v_{{\rm{S0}}}^2} \right){{\cos }^2}\left( {\theta + {\theta _R}} \right) - \sqrt D } \right\}} $

(A2)

其中:

$ \begin{array}{l} D =\\ {\left\{ {\left[ {\left( {1 + 2\varepsilon } \right)\rho v_{{\rm{P0}}}^2 + \rho v_{{\rm{S0}}}^2} \right]{{\sin }^2}\left( {\theta + {\theta _R}} \right) + \left( {\rho v_{{\rm{P0}}}^2 + \rho v_{{\rm{S0}}}^2} \right){{\cos }^2}\left( {\theta + {\theta _R}} \right)} \right\}^2} + 4{\rho ^2}v_{{\rm{P0}}}^2\left( {\delta - \varepsilon } \right) \cdot \\ \;\;\;\;\;\;2\left( {1 - \frac{{v_{{\rm{S0}}}^2}}{{v_{{\rm{P0}}}^2}}} \right){\sin ^2}\left( {\theta + {\theta _R}} \right){\cos ^2}\left( {\theta + {\theta _R}} \right) \end{array} $

(A3)

对比qP波和qSV波相速度精确表达式, 发现二者有一个共同项, 即:

$ \left[ {\left( {1 + 2\varepsilon } \right)\rho v_{{\rm{P0}}}^2 + \rho v_{{\rm{S0}}}^2} \right]{\sin ^2}\left( {\theta + {\theta _R}} \right) + \left( {\rho v_{{\rm{P0}}}^2 + \rho v_{{\rm{S0}}}^2} \right){\cos ^2}\left( {\theta + {\theta _R}} \right) $

(A4)

如果从(A3)式中可以直接提出以下项, 即:

$ \left[ {\left( {1 + 2\varepsilon } \right)\rho v_{{\rm{P0}}}^2 - \rho v_{{\rm{S0}}}^2} \right]{\sin ^2}\left( {\theta + {\theta _R}} \right) + \left( {\rho v_{{\rm{P0}}}^2 - \rho v_{{\rm{S0}}}^2} \right){\cos ^2}\left( {\theta + {\theta _R}} \right) $

(A5)

就能抵消(A4)中的部分项从而简化qP波和qSV波的相速度表达式。

为了使用近似的配方法将(A5)从D中整体提出, 令(A5)式为a, 即:

$ \begin{array}{l} a = \left[ {\left( {1 + 2\varepsilon } \right)\rho v_{{\rm{P0}}}^2 - \rho v_{{\rm{S0}}}^2} \right]{\sin ^2}\left( {\theta + {\theta _R}} \right) + \left( {\rho v_{{\rm{P0}}}^2 - \rho v_{{\rm{S0}}}^2} \right){\cos ^2}\left( {\theta + {\theta _R}} \right)\\ {\rm{ = }}\left( {\rho v_{{\rm{P0}}}^2 - \rho v_{{\rm{S0}}}^2} \right){\sin ^2}\left( {\theta + {\theta _R}} \right) + 2\varepsilon \rho v_{{\rm{P0}}}^2{\sin ^2}\left( {\theta + {\theta _R}} \right) + \left( {\rho v_{{\rm{P0}}}^2 - \rho v_{{\rm{S0}}}^2} \right){\cos ^2}\left( {\theta + {\theta _R}} \right)\\ {\rm{ = }}\rho v_{{\rm{P0}}}^2\left[ {\left( {1 - \frac{{v_{{\rm{P0}}}^2}}{{v_{{\rm{S0}}}^2}}} \right) + 2\varepsilon {{\sin }^2}\left( {\theta + {\theta _R}} \right)} \right] \end{array} $

(A6)

将(A6)式代入(A3)式进行整理可得:

$ \begin{array}{l} D = {a^2} + 4{\rho ^2}v_{{\rm{P0}}}^4\left( {\delta - \varepsilon } \right) \cdot 2\left( {1 - \frac{{v_{{\rm{S0}}}^2}}{{v_{{\rm{P0}}}^2}}} \right) \cdot {\sin ^2}\left( {\theta + {\theta _R}} \right){\cos ^2}\left( {\theta + {\theta _R}} \right)\\ = {a^2} + 2 \cdot \rho v_{{\rm{P0}}}^2\left( {1 - \frac{{v_{{\rm{S0}}}^2}}{{v_{{\rm{P0}}}^2}}} \right) \cdot 2\rho v_{{\rm{P0}}}^2\left( {\delta - \varepsilon } \right){\sin ^2}\left( {\theta + {\theta _R}} \right){\cos ^2}\left( {\theta + {\theta _R}} \right) + \\ 4{\rho ^2}v_{{\rm{P0}}}^4\left( {\delta - \varepsilon } \right){\sin ^2}\left( {\theta + {\theta _R}} \right){\cos ^2}\left( {\theta + {\theta _R}} \right) \cdot \left( {1 - \frac{{v_{{\rm{S0}}}^2}}{{v_{{\rm{P0}}}^2}}} \right) \end{array} $

(A7)

如果令:

$ b = 2\rho v_{{\rm{P0}}}^2\left( {\delta - \varepsilon } \right){\sin ^2}\left( {\theta + {\theta _R}} \right){\cos ^2}\left( {\theta + {\theta _R}} \right) $

(A8)

并试算(a+b)², 有:

$ \begin{array}{l} {\left( {a + b} \right)^2} = {a^2} + 2ab + {b^2}\\ = {a^2} + 2\rho v_{{\rm{P0}}}^2\left[ {\left( {1 - \frac{{v_{{\rm{S0}}}^2}}{{v_{{\rm{P0}}}^2}}} \right) + 2\varepsilon {{\sin }^2}\left( {\theta + {\theta _R}} \right)} \right] \cdot 2\rho v_{{\rm{P0}}}^2\left( {\delta - \varepsilon } \right){\sin ^2}\left( {\theta + {\theta _R}} \right){\cos ^2}\left( {\theta + {\theta _R}} \right) + \\ 4{\rho ^2}v_{{\rm{P0}}}^4{\left( {\delta - \varepsilon } \right)^2}{\sin ^4}\left( {\theta + {\theta _R}} \right){\cos ^4}\left( {\theta + {\theta _R}} \right)\\ = {a^2} + 4{\rho ^2}v_{{\rm{P0}}}^4\left( {\delta - \varepsilon } \right){\sin ^2}\left( {\theta + {\theta _R}} \right){\cos ^2}\left( {\theta + {\theta _R}} \right) \cdot \left( {1 - \frac{{v_{{\rm{S0}}}^2}}{{v_{{\rm{P0}}}^2}}} \right) + \\ 4{\rho ^2}v_{{\rm{P0}}}^4\left( {\delta - \varepsilon } \right){\sin ^2}\left( {\theta + {\theta _R}} \right){\cos ^2}\left( {\theta + {\theta _R}} \right)\\ \left[ {2\varepsilon {{\sin }^2}\left( {\theta + {\theta _R}} \right) + \left( {\delta - \varepsilon } \right){{\sin }^2}\left( {\theta + {\theta _R}} \right){{\cos }^2}\left( {\theta + {\theta _R}} \right)} \right] \end{array} $

(A9)

对比(A7)式与(A9)式, 不难发现, 两个表达式之间存在差异的项为:

$ 4{\rho ^2}v_{{\rm{P0}}}^4\left( {\delta - \varepsilon } \right){\sin ^2}\left( {\theta + {\theta _R}} \right){\cos ^2}\left( {\theta + {\theta _R}} \right) \cdot \left( {1 - \frac{{v_{{\rm{S0}}}^2}}{{v_{{\rm{P0}}}^2}}} \right) $

(A10)

和

$ 4{\rho ^2}v_{{\rm{P0}}}^4\left( {\delta - \varepsilon } \right){\sin ^2}\left( {\theta + {\theta _R}} \right){\cos ^2}\left( {\theta + {\theta _R}} \right)\\ \left[ {2\varepsilon {{\sin }^2}\left( {\theta + {\theta _R}} \right) + \left( {\delta - \varepsilon } \right){{\sin }^2}\left( {\theta + {\theta _R}} \right){{\cos }^2}\left( {\theta + {\theta _R}} \right)} \right] $

(A11)

如果令(A10)式和(A11)近似相等, 即:

$ \begin{array}{l} 4{\rho ^2}v_{{\rm{P0}}}^4\left( {\delta - \varepsilon } \right){\sin ^2}\left( {\theta + {\theta _R}} \right){\cos ^2}\left( {\theta + {\theta _R}} \right)\left( {1 - \frac{{v_{{\rm{S0}}}^2}}{{v_{{\rm{P0}}}^2}}} \right)\\ \approx 4{\rho ^2}v_{p{\rm{0}}}^4\left( {\delta - \varepsilon } \right){\sin ^2}\left( {\theta + {\theta _R}} \right){\cos ^2}\left( {\theta + {\theta _R}} \right)\\ \left[ {2\varepsilon {{\sin }^2}\left( {\theta + {\theta _R}} \right) + \left( {\delta - \varepsilon } \right){{\sin }^2}\left( {\theta + {\theta _R}} \right){{\cos }^2}\left( {\theta + {\theta _R}} \right)} \right] \end{array} $

(A12)

就能使用配方法消除$ \sqrt D $中的根号。此时:

$ \begin{array}{l} \sqrt D \approx \sqrt {{{\left( {a + b} \right)}^2}} = a + b\\ \;\;\;\;\;\;\;\;\; = \left[ {\left( {1 + 2\varepsilon } \right)\rho v_{{\rm{P0}}}^2 - \rho v_{{\rm{S0}}}^2} \right]{\sin ^2}\left( {\theta + {\theta _R}} \right) + \left( {\rho v_{{\rm{P0}}}^2 - \rho v_{{\rm{S0}}}^2} \right){\cos ^2}\left( {\theta + {\theta _R}} \right) + \\ \;\;\;\;\;\;\;\;\;\;\;\;2\rho v_{{\rm{P0}}}^2\left( {\delta - \varepsilon } \right){\sin ^2}\left( {\theta + {\theta _R}} \right){\cos ^2}\left( {\theta + {\theta _R}} \right) \end{array} $

(A13)

由于该方法在进行配方时采用了近似, 所以称该方法为近似的配方法。

将(A13)式代入(A1)式, 可以得到TTI介质qP波一般性近似相速度, 即:

$ {v_{\rm{P}}} \approx \\ {\left\{ {\frac{1}{{2\rho }}\left\{ {\left[ {\left( {1 + 2\varepsilon } \right)\rho v_{{\rm{P0}}}^2 - \rho v_{{\rm{S0}}}^2} \right]{{\sin }^2}\left( {\theta + {\theta _R}} \right) + \left( {\rho v_{{\rm{P0}}}^2 - \rho v_{{\rm{S0}}}^2} \right){{\cos }^2}\\ \left( {\theta + {\theta _R}} \right) + \sqrt {{{\left( {a + b} \right)}^2}} } \right\}} \right\}^{\frac{1}{2}}}\\ \;\;\;\; = \left\{ {\frac{1}{{2\rho }}\left\{ {\left[ {\left( {1 + 2\varepsilon } \right)\rho v_{{\rm{P0}}}^2 - \rho v_{{\rm{S0}}}^2} \right]{{\sin }^2}\left( {\theta + {\theta _R}} \right) + \left( {\rho v_{{\rm{P0}}}^2 - \rho v_{{\rm{S0}}}^2} \right){{\cos }^2}\\ \left( {\theta + {\theta _R}} \right) + \left[ {\left( {1 + 2\varepsilon } \right)\rho v_{{\rm{P0}}}^2 - \rho v_{{\rm{S0}}}^2} \right] \cdot } \right.} \right.\\ \;\;\;\;{\left. {\left. {{{\sin }^2}\left( {\theta + {\theta _R}} \right) + \left( {\rho v_{{\rm{P0}}}^2 - \rho v_{{\rm{S0}}}^2} \right){{\cos }^2}\left( {\theta + {\theta _R}} \right) + 2\rho v_{{\rm{P0}}}^2\left( {\delta - \varepsilon } \right){{\sin }^2}\left( {\theta + {\theta _R}} \right){{\cos }^2}\\ \left( {\theta + {\theta _R}} \right)} \right\}} \right\}^{\frac{1}{2}}}\\ \;\;\;\; = {\left\{ {\frac{1}{{2\rho }}\left[ {2\left( {1 + 2\varepsilon } \right)\rho v_{{\rm{P0}}}^2{{\sin }^2}\left( {\theta + {\theta _R}} \right) + 2\rho v_{{\rm{P0}}}^2{{\cos }^2}\left( {\theta + {\theta _R}} \right) + 2\rho v_{{\rm{P0}}}^2\\ \left( {\delta - \varepsilon } \right){{\sin }^2}\left( {\theta + {\theta _R}} \right){{\cos }^2} \left( {\theta + {\theta _R}} \right)} \right]} \right\}^{\frac{1}{2}}}\\ \;\;\;\; = v_{{\rm{P0}}}{\left[ {\left( {1 + 2\varepsilon } \right){{\sin }^2}\left( {\theta + {\theta _R}} \right) + {{\cos }^2}\left( {\theta + {\theta _R}} \right) + \left( {\delta - \varepsilon } \right){{\sin }^2}\left( {\theta + {\theta _R}} \right){{\cos }^2}\left( {\theta + {\theta _R}} \right)} \right]^{\frac{1}{2}}}\\ \;\;\;\; = {v_{{\rm{P0}}}}{\left[ {1 + 2\varepsilon {{\sin }^2}\left( {\theta + {\theta _R}} \right) + \left( {\delta - \varepsilon } \right){{\sin }^2}\left( {\theta + {\theta _R}} \right){{\cos }^2}\left( {\theta + {\theta _R}} \right)} \right]^{\frac{1}{2}}} $

(A14)

同理, 可以得到TTI介质qSV波一般性近似相速度, 即:

$ {v_{{\rm{SV}}}} \approx {v_{{\rm{S0}}}}\sqrt {1 + \frac{{v_{{\rm{P0}}}^2}}{{v_{{\rm{S0}}}^2}}\left( {\varepsilon - \delta } \right){{\sin }^2}\left( {\theta + {\theta _R}} \right){{\cos }^2}\left( {\theta + {\theta _R}} \right)} $

(A15)