CN107861153A - 薄层pp波反射系数的快速求解方法 - Google Patents
薄层pp波反射系数的快速求解方法 Download PDFInfo
- Publication number
- CN107861153A CN107861153A CN201711009284.0A CN201711009284A CN107861153A CN 107861153 A CN107861153 A CN 107861153A CN 201711009284 A CN201711009284 A CN 201711009284A CN 107861153 A CN107861153 A CN 107861153A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mfrac
- msubsup
- msup
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 18
- 239000011159 matrix material Substances 0.000 claims abstract description 53
- 230000005540 biological transmission Effects 0.000 claims abstract description 40
- 230000014509 gene expression Effects 0.000 claims abstract description 30
- 239000010409 thin film Substances 0.000 claims abstract description 16
- AYEKOFBPNLCAJY-UHFFFAOYSA-O thiamine pyrophosphate Chemical compound CC1=C(CCOP(O)(=O)OP(O)(O)=O)SC=[N+]1CC1=CN=C(C)N=C1N AYEKOFBPNLCAJY-UHFFFAOYSA-O 0.000 claims description 10
- 230000008859 change Effects 0.000 claims description 3
- 238000001615 p wave Methods 0.000 claims description 2
- 230000010363 phase shift Effects 0.000 claims description 2
- 239000010410 layer Substances 0.000 description 115
- 239000003245 coal Substances 0.000 description 5
- 239000011435 rock Substances 0.000 description 4
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 238000010304 firing Methods 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 239000011229 interlayer Substances 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 238000009738 saturating Methods 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 235000009470 Theobroma cacao Nutrition 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 244000240602 cacao Species 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 230000011514 reflex Effects 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/63—Seismic attributes, e.g. amplitude, polarity, instant phase
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Pharmaceuticals Containing Other Organic And Inorganic Compounds (AREA)
Abstract
本发明公开了一种薄层PP波反射系数的快速求解方法,包括:根据薄层模型,并假设薄层模型的薄层厚度远小于波数的倒数,以取得极薄层反射透射矩阵方程;对入射角取负值,并对极薄层反射透射矩阵方程进行三角函数奇偶性变换得到系数矩阵方程,且将系数矩阵方程与极薄层反射透射矩阵方程对比,以建立反透射系数与入射角正弦函数的奇偶关系式;根据奇偶关系式、极薄层反射透射矩阵方程的常数项、入射角正弦函数一次项及入射角正弦函数二次项,以获得入射角正弦函数幂级数的系数;根据入射角正弦函数幂级数的系数,以获得薄层PP波反射系数的解析解近似式。通过本发明,以解决现有技术存在的薄层反射系数计算复杂、难以用于AVA反演的问题。
Description
技术领域
本发明涉及地震勘探的技术领域,尤其涉及一种薄层PP波反射系数的快速求解方法。
背景技术
随着我国油气田、煤田等资源开发的不断深入,薄(互)储层、岩性、碳酸盐岩、非常规油气等复杂油气藏及煤储层成为勘探的重点。在我国东部地区,绝大多数的中、新生代陆相含油盆地大都以薄层砂、泥岩沉积为主,并夹有少量的碳酸盐岩、页岩及膏盐层,这些地层岩性和厚度横向变化均较大,而且厚度远远低于常规地震勘探的垂向分辨率,以薄层的形式存在。我国煤系地层绝大部分煤层及围岩物性差异较大,煤层顶、底岩层为良好的波阻抗分界面,可以形成较强的反射波检测。由于煤层厚度多在2~10m之间,且常以薄互层的形式存在,反射响应非单一界面产生,而是顶底反射、层间多次波叠加而成的复合波。因此,基于Zoeppritz方程的常规AVO分析和反演方法,对于薄层问题将不再适用。薄层反射地震理论的研究成为推动薄层AVO技术发展的基础。
针对薄层AVO理论的研究,Thomson(1950)、Haskell(1953)、Brekhovskikh(1960)给出了层状弹性介质的位移位反透射矩阵方程,Pan等人(2013)将Brekhovskikh的公式应用于单薄层情况,但其复杂的传播矩阵形式使得以上公式难以应用于实际。Meissner等人(1969)、Juhlin等人(1993)给出了基于时延的薄层反透射系数公式,将层间多次波依次叠加到薄层顶界面的反射波上,且表达式以界面处的反透射系数表示。Yang等人(2017)在Juhlin等人研究的基础上给出三个忽略多次波影响的PP波反射系数解析近似式。Kennett(1983)建立了针对于多层介质中球面波入射情况下反射系数的层间递归关系。Liu等人(2003)在忽略层间多次波及转换波的情况下给出了单薄层的反射系数,Rubino等人(2009)将其推广到弹性介质,并利用振幅谱实现了薄层谱反演。
然而,现有薄层反射系数及近似不是过于复杂难以实际应用,就是忽略了层间多次波及转换波影响了近似精度。
发明内容
本发明的主要目的在于提供一种薄层PP波反射系数的快速求解方法,以解决现有技术存在的薄层反射系数计算复杂、难以用于AVA反演的问题。
为解决上述问题,本发明实施例提供一种薄层PP波反射系数的快速求解方法,包括:根据薄层模型,并假设所述薄层模型的薄层厚度远小于波数的倒数,以取得极薄层反射透射矩阵方程;对入射角取负值,并对所述极薄层反射透射矩阵方程进行三角函数奇偶性变换得到系数矩阵方程,且将所述系数矩阵方程与所述极薄层反射透射矩阵方程对比,以建立反透射系数与入射角正弦函数的奇偶关系式;根据所述奇偶关系式、所述极薄层反射透射矩阵方程的常数项、入射角正弦函数一次项及入射角正弦函数二次项,以获得入射角正弦函数幂级数的系数;根据所述入射角正弦函数幂级数的系数,以获得薄层PP波反射系数的解析解近似式。
根据本发明的技术方案,通过假设所述薄层模型的薄层厚度远小于波数的倒数,以取得极薄层反射透射矩阵方程;对入射角取负值,并对所述极薄层反射透射矩阵方程进行三角函数奇偶性变换得到系数矩阵方程,且将所述系数矩阵方程与所述极薄层反射透射矩阵方程对比,以建立反透射系数与入射角正弦函数的奇偶关系式;根据所述奇偶关系式、所述极薄层反射透射矩阵方程的常数项、入射角正弦函数一次项及入射角正弦函数二次项,以获得入射角正弦函数幂级数的系数;根据所述入射角正弦函数幂级数的系数,以获得薄层PP波反射系数的解析解近似式。从而,薄层PP波反射系数计算更为简单且形式较为简洁,可以更有效地实现薄层正演分析,也利于薄层AVO反演的实现。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1是根据本发明实施例的薄层PP波反射系数的快速求解方法的流程图;
图2是根据本发明实施例的薄层模型的示意图;
图3是根据本发明实施例的极薄层近似、解析解近似对精准值逼近程度对比图;
图4是根据本发明实施例的模型1的反射系数振幅和相位近似误差图;
图5是根据本发明实施例的模型2的反射系数振幅和相位近似误差图;
图6是根据本发明实施例的模型3的反射系数振幅和相位近似误差图;
图7是根据本发明实施例的模型4的反射系数振幅和相位近似误差图。
具体实施方式
本发明的主要思想在于,基于假设所述薄层模型的薄层厚度远小于波数的倒数,以取得极薄层反射透射矩阵方程;对入射角取负值,并对所述极薄层反射透射矩阵方程进行三角函数奇偶性变换得到系数矩阵方程,且将所述系数矩阵方程与所述极薄层反射透射矩阵方程对比,以建立反透射系数与入射角正弦函数的奇偶关系式;根据所述奇偶关系式、所述极薄层反射透射矩阵方程的常数项、入射角正弦函数一次项及入射角正弦函数二次项,以获得入射角正弦函数幂级数的系数;根据所述入射角正弦函数幂级数的系数,以获得薄层PP波反射系数的解析解近似式。从而,薄层PP波反射系数计算更为简单且形式较为简洁,可以更有效地实现薄层正演分析,也利于薄层AVO反演的实现。
为使本发明的目的、技术方案和优点更加清楚,以下结合附图及具体实施例,对本发明做进一步地详细说明。
图1是根据本发明实施例的薄层PP波反射系数的快速求解方法的流程图。
在步骤S102中,根据薄层模型,并假设所述薄层模型的薄层厚度远小于波数的倒数,以取得极薄层反射透射矩阵方程。假设本实施例的薄层模型如图2所示,其由三层各向同性弹性介质构成了单薄层模型,其中,P为入射P波,PP为反射PP波、PS为反射PS波,vPi,vSi,ρi分别为各层的纵波、横波速度及密度,i=1,2,3,h为薄层厚度,θ为入射角。由图2的薄层模型,并假设中间层厚度远小于波数的倒数时,可获得极薄层反射透射矩阵方程,如式(1)所示:
M[RPP RPS TPP TPS]T=N, (1)
其中:RPP、RPS、TPP、TPS分别为薄层反射透射系数,M为4*4矩阵,N为4*1矩阵,并且矩阵M、N中的参数分别如(1a)~(1t)所示:
vPi、vSi、ρi分别为各层的纵波速度、横波速度及密度,且它们的下角标(即i)取值为1、2、3(即i=1,2,3),分别对应单薄层的三层介质、θ为入射角、h为薄层厚度,ω为圆频率,ω=2πf,其中f为入射波的频率,
由上述内容可知,本实施例所提供的极薄层反射透射矩阵方程具有较为简洁的、类似于Zoeppritz方程的形式,该近似方程避免了传播矩阵的运算,有效地简化了薄层反透射系数的求解。同时,式(1)中的角度项都以入射角θ表征,更有效地反映了反透射系数随入射角的变化规律。
在步骤S104中,对入射角取负值,并对所述极薄层反射透射矩阵方程进行三角函数奇偶性变换得到系数矩阵方程,且将所述系数矩阵方程与所述极薄层反射透射矩阵方程对比,以建立反透射系数与入射角正弦函数的奇偶关系式。
当式(1)中的sinθ取值为-sinθ时,sin2θ相应取值为-sin2θ,而cosθ取值不变,并设定此时极薄层的反透射系数分别为R'PP、R'PS、T'PP、T'PS,将上述参数代入式(1)并整理可得系数矩阵方程,如式(2)所示:
M[R'PP -R'PS T'PP -T'PS]T=N, (2)
比较式(1)和式(2),可以获得式(3),如下所示:
R'PP=RPP、R'PS=-RPS、T'PP=TPP、T'PS=-TPS, (3)
由此可知,极薄层非转换波的反射和透射系数是sinθ(或射线参数)的偶函数,而转换波的反射和透射系数是sinθ(或射线参数)的奇函数。
并且,利用反射和透射系数的奇偶性可以将它们表示成sinθ的幂级数形式,即反透射系数与入射角正弦函数的奇偶关系式,如式(4)所示:
其中,A2n、A2n+1、B2n、B2n+1分别为RPP、RPS、TPP、TPS反透射系数幂级数表达式的系数。当入射角较小时,sinnθ随着n的增加迅速衰减,则近似求解出系数A2n、A2n+1、B2n、B2n+1,根据精度需求对(4)式进行高阶舍弃便可获得相应的薄层近似公式。
在步骤S106中,根据所述奇偶关系式、所述极薄层反射透射矩阵方程的常数项、入射角正弦函数一次项及入射角正弦函数二次项,以获得入射角正弦函数幂级数的系数。
式(1)中考虑垂直入射情况,即θ=0°,对应极薄层反射透射矩阵方程的常数项。此时假设RPP=A0、TPP=B0、RPS=0、TPS=0,式(1)可化简为式(5),如下所示:
其中,τ=ωh/vP2为薄层内纵波的单程相移,zP=ρvP为纵波阻抗,它们的下角标取值为1、2、3,分别对应单薄层的三层介质。并且,对式(5)计算可可得式(6a)和式(6b),如下所示:
m0B0=2(1+τ2), (6b)
其中,为系数A0、B0的分母。
另外,式(1)中只考虑到sinθ一次项,即sinθ以上高次项舍弃,对应极薄层反射透射矩阵方程的入射角正弦函数一次项。此时假设RPP=A0、TPP=B0、RPS=A1sinθ、TPS=B1sinθ,式(1)可化简为式(7),如下所示:
其中,zS=ρvS为横波阻抗,r=vP/vS为纵横波速度比,它们的下角标取值为1、2、3,分别对应薄层的三层介质;lk=vPk+1/vPk为相邻层纵波速度的比值,其中k取值为1或2。并且,对(7)计算可得式(8a)和式(8b),如下所示:
其中,为系数A1、B1的分母,和为简化替换参量。
此外,式(1)中只考虑到sinθ二次项,对应极薄层反射透射矩阵方程的入射角正弦函数二次项。此时假设RPP=A0+A2sin2θ、TPP=B0+B2sin2θ、RPS=A1sinθ、TPS=B1sinθ,式(1)可化简为式(9),如下所示:
并且,对(9)计算可得式(10),如下所示:
由上述的计算方式,可以获得入射角正弦函数幂级数的系数A0、B0、A1、B1、A2。
在步骤S108中,根据所述入射角正弦函数幂级数的系数,以获得薄层PP波反射系数的解析解近似式。其中,将式(6a)、式(10)整合,以获得薄层PP波反射系数的解析解近似式,如式(10)所示:
RPP=A0+A2sin2θ。 (10)
将RPP的AVA曲线转换为RPP随sin2θ的变化曲线时,参数A0为该曲线的截距,而参数A2为该曲线的斜率,由此从曲线上提取截距、斜率等参数,使得薄层弹性参数的反演成为可能。
上述已说明了薄层PP波反射系数的快速求解方法,以下将提供一些实例来验证上述方法的适用性及有效性。
为了验证RPP解析解近似式的适用性及有效性,设定四个不同类型的薄层模型如表1所示,分别代表着高阻抗薄层、阻抗递增薄层、低阻抗薄层、阻抗递减薄层;同时模型1和模型3为反极性薄层,模型2和模型4为同极性的薄层。
表1薄层模型参数
在表1中,速度的单位为:m/s,密度的单位为:g/cm3。
遵循Widess(1973)关于薄层的定义:薄层厚度小于λ/8的单层为薄层,其中λ为地震波在薄层中主频对应的波长。由于单频讨论,因此取薄层厚度由λ/8变化到λ/60,其中λ为单频下PP波波长。以模型1的λ/10、λ/20为例,对比分析极薄层近似、解析解近似对精准PP波反射系数的近似逼近程度。
图3是根据本发明实施例的极薄层近似、解析解近似对精准值逼近程度对比图。在图3中,标号310表示精准值,标号320表示极薄层近似,标号330表示解析解近似,实线表示λ/10,虚线表示λ/20。如图3所示,在小角度范围内,解析解近似完全可以达到极薄层近似对精准反射值的近似逼近程度,且薄层厚度越小,解析解近似对精准PP波反射系数的近似程度越高。类似地,模型2~4也反映了类似的规律及效果,故未附图。
由解析解近似产生的近似误差如图4~7所示,分别对应四种薄层模型的情况。
模型1的顶、底界面处的阻抗差分别为0.3333、–0.4186,且模型1为极性相反(+,–)的高阻抗薄层,并存在30°临界角。如图4所示,解析解近似误差随着薄层厚度的减少(即波长/厚度的增加)而减小。以10%为最大阀值,薄层厚度为λ/8情况下在角度小于23°时满足条件,且随着厚度的减少,满足条件的角度范围有扩大的趋势。对于相位而言,除了λ/8情况以外,近似误差基本都在10%以内。
模型2的顶、底界面处的阻抗差分别为0.1586、0.1845,且模型2为极性相同(+,+)的阻抗递增薄层,并存在30°临界角。如图5所示,对于不同薄层厚度的情况,入射角小于20°时,振幅近似误差都小于5%;入射角大于20°时,近似误差随着入射角的增加而迅速增大。
模型3的顶、底界面处的阻抗差分别为–0.2105、0.3695,且模型3为极性相反(–,+)的低阻抗薄层,并存在46.16°临界角。如图6所示,随着入射角的增加,解析解近似误差增大,以10%为最大阀值,则在入射角小于25°的情况下满足精度要求。对于不同薄层厚度的情况(除λ/8),在临界角小于25°的情况下,振幅近似误差都小于10%;对于薄层厚度为λ/8的情况,在入射角小于20°情况下,近似误差小于10%。振幅近似误差随着入射角的增加而增大,而相位近似误差在不同厚度及入射角情况下都小于10%。
模型4的顶、底界面处的阻抗差分别为–0.1845、–0.1836,且模型4为极性相同(–,–)的阻抗递减薄层。如图7所示,解析解近似对模型4的PP波反射系数近似精度很高。在讨论的薄层厚度及角度范围内,振幅近似误差都小于10%,且在薄层厚度小于λ/8情况下,近似误差小于5%;在讨论的薄层厚度及角度范围内,相位近似误差都小于5%。
综上所述,根据本发明的技术方案,通过假设所述薄层模型的薄层厚度远小于波数的倒数,以取得极薄层反射透射矩阵方程;对入射角取负值,并对所述极薄层反射透射矩阵方程进行三角函数奇偶性变换得到系数矩阵方程,且将所述系数矩阵方程与所述极薄层反射透射矩阵方程对比,以建立反透射系数与入射角正弦函数的奇偶关系式;根据所述奇偶关系式、所述极薄层反射透射矩阵方程的常数项、入射角正弦函数一次项及入射角正弦函数二次项,以获得入射角正弦函数幂级数的系数;根据所述入射角正弦函数幂级数的系数,以获得薄层PP波反射系数的解析解近似式。从而,薄层PP波反射系数计算更为简单且形式较为简洁,可以更有效地实现薄层正演分析,也利于薄层AVO反演的实现。
以上所述仅为本发明的实施例,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的权利要求范围之内。
Claims (5)
1.一种薄层PP波反射系数的快速求解方法,其特征在于,包括以下步骤:
根据薄层模型,并假设所述薄层模型的薄层厚度远小于波数的倒数,以取得极薄层反射透射矩阵方程;
对入射角取负值,并对所述极薄层反射透射矩阵方程进行三角函数奇偶性变换得到系数矩阵方程,且将所述系数矩阵方程与所述极薄层反射透射矩阵方程对比,以建立反透射系数与入射角正弦函数的奇偶关系式;
根据所述奇偶关系式、所述极薄层反射透射矩阵方程的常数项、入射角正弦函数一次项及入射角正弦函数二次项,以获得入射角正弦函数幂级数的系数;
根据所述入射角正弦函数幂级数的系数,以获得薄层PP波反射系数的解析解近似式。
2.根据权利要求1所述的薄层PP波反射系数的快速求解方法,其特征在于,所述极薄层反射透射矩阵方程满足如下公式:
M[RPP RPS TPP TPS]T=N,
其中:RPP、RPS、TPP、TPS分别为薄层反射透射系数,M为4*4矩阵,N为4*1矩阵,且有:
<mrow>
<msub>
<mi>m</mi>
<mn>11</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&theta;</mi>
<mo>+</mo>
<mi>j</mi>
<mfrac>
<mrow>
<mi>&omega;</mi>
<mi>h</mi>
</mrow>
<msub>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mo>-</mo>
<mfrac>
<msub>
<mi>&rho;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>&rho;</mi>
<mn>2</mn>
</msub>
</mfrac>
<mfrac>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<mo>)</mo>
</mrow>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mn>2</mn>
<mi>&theta;</mi>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>m</mi>
<mn>21</mn>
</msub>
<mo>=</mo>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&theta;</mi>
<mo>+</mo>
<mi>j</mi>
<mfrac>
<mrow>
<mi>&omega;</mi>
<mi>h</mi>
</mrow>
<msub>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mo>&lsqb;</mo>
<mfrac>
<msub>
<mi>&rho;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>&rho;</mi>
<mn>2</mn>
</msub>
</mfrac>
<mfrac>
<msubsup>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mn>2</mn>
<mfrac>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<mo>+</mo>
<mn>2</mn>
<mfrac>
<msub>
<mi>&rho;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>&rho;</mi>
<mn>2</mn>
</msub>
</mfrac>
<mfrac>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<mo>)</mo>
</mrow>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
<mo>&rsqb;</mo>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>m</mi>
<mn>31</mn>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mn>2</mn>
<mfrac>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>j</mi>
<mfrac>
<mrow>
<mi>&omega;</mi>
<mi>h</mi>
</mrow>
<msub>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>&rho;</mi>
<mn>2</mn>
</msub>
<msub>
<mi>&rho;</mi>
<mn>1</mn>
</msub>
</mfrac>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&theta;</mi>
<mo>-</mo>
<mfrac>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mn>2</mn>
<mi>&theta;</mi>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>m</mi>
<mn>41</mn>
</msub>
<mo>=</mo>
<mi>sin</mi>
<mn>2</mn>
<mi>&theta;</mi>
<mo>+</mo>
<mi>j</mi>
<mfrac>
<mrow>
<mi>&omega;</mi>
<mi>h</mi>
</mrow>
<msub>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mfrac>
<msub>
<mi>&rho;</mi>
<mn>2</mn>
</msub>
<msub>
<mi>&rho;</mi>
<mn>1</mn>
</msub>
</mfrac>
<mfrac>
<msubsup>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<mo>+</mo>
<mn>4</mn>
<mfrac>
<msub>
<mi>&rho;</mi>
<mn>2</mn>
</msub>
<msub>
<mi>&rho;</mi>
<mn>1</mn>
</msub>
</mfrac>
<mfrac>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<mrow>
<mo>(</mo>
<mrow>
<mfrac>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
</mrow>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mrow>
<mfrac>
<mrow>
<mn>2</mn>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
<msubsup>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mrow>
<mfrac>
<msubsup>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<mo>-</mo>
<mn>2</mn>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mi>sin</mi>
<mi>&theta;</mi>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>m</mi>
<mn>12</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mrow>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
<mo>-</mo>
<mi>j</mi>
<mfrac>
<mrow>
<mi>&omega;</mi>
<mi>h</mi>
</mrow>
<msub>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mfrac>
<msub>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mo>&lsqb;</mo>
<mfrac>
<msub>
<mi>&rho;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>&rho;</mi>
<mn>2</mn>
</msub>
</mfrac>
<mfrac>
<msubsup>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mn>2</mn>
<mfrac>
<msub>
<mi>&rho;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>&rho;</mi>
<mn>2</mn>
</msub>
</mfrac>
<mfrac>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
</mrow>
<mo>)</mo>
</mrow>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
<mo>&rsqb;</mo>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>m</mi>
<mn>22</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<msub>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&theta;</mi>
<mo>-</mo>
<mi>j</mi>
<mfrac>
<mrow>
<mi>&omega;</mi>
<mi>h</mi>
</mrow>
<msub>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mn>2</mn>
<mfrac>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<mo>+</mo>
<mn>2</mn>
<mfrac>
<msub>
<mi>&rho;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>&rho;</mi>
<mn>2</mn>
</msub>
</mfrac>
<mfrac>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mrow>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&theta;</mi>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>m</mi>
<mn>32</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<mn>2</mn>
<mfrac>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mrow>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&theta;</mi>
<mo>-</mo>
<mi>j</mi>
<mfrac>
<mrow>
<mi>&omega;</mi>
<mi>h</mi>
</mrow>
<msub>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mfrac>
<msub>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mfrac>
<msub>
<mi>&rho;</mi>
<mn>2</mn>
</msub>
<msub>
<mi>&rho;</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>-</mo>
<mn>2</mn>
<mfrac>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&theta;</mi>
<mo>,</mo>
</mrow>
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>m</mi>
<mn>42</mn>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mn>2</mn>
<mfrac>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
</mrow>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>j</mi>
<mfrac>
<mrow>
<mi>&omega;</mi>
<mi>h</mi>
</mrow>
<msub>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mo>&lsqb;</mo>
<mfrac>
<msub>
<mi>&rho;</mi>
<mn>2</mn>
</msub>
<msub>
<mi>&rho;</mi>
<mn>1</mn>
</msub>
</mfrac>
<mfrac>
<msubsup>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<mo>+</mo>
<mn>4</mn>
<mfrac>
<msub>
<mi>&rho;</mi>
<mn>2</mn>
</msub>
<msub>
<mi>&rho;</mi>
<mn>1</mn>
</msub>
</mfrac>
<mfrac>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<mrow>
<mo>(</mo>
<mrow>
<mfrac>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
</mrow>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mrow>
<mfrac>
<mrow>
<mn>2</mn>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
<msubsup>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
</mrow>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
<mo>&rsqb;</mo>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>m</mi>
<mn>13</mn>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>3</mn>
</mrow>
</msub>
<msub>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&theta;</mi>
<mo>,</mo>
<msub>
<mi>m</mi>
<mn>23</mn>
</msub>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<msubsup>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>3</mn>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mrow>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
<mo>,</mo>
<msub>
<mi>m</mi>
<mn>33</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<msub>
<mi>&rho;</mi>
<mn>3</mn>
</msub>
<msub>
<mi>&rho;</mi>
<mn>1</mn>
</msub>
</mfrac>
<mfrac>
<msub>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>3</mn>
</mrow>
</msub>
<msub>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mn>2</mn>
<mfrac>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>3</mn>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>m</mi>
<mn>43</mn>
</msub>
<mo>=</mo>
<mn>2</mn>
<mfrac>
<msub>
<mi>&rho;</mi>
<mn>3</mn>
</msub>
<msub>
<mi>&rho;</mi>
<mn>1</mn>
</msub>
</mfrac>
<mfrac>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>3</mn>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<msubsup>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>3</mn>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mrow>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&theta;</mi>
<mo>,</mo>
<msub>
<mi>m</mi>
<mn>14</mn>
</msub>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>3</mn>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mrow>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
<mo>,</mo>
<msub>
<mi>m</mi>
<mn>24</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<msub>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>3</mn>
</mrow>
</msub>
<msub>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&theta;</mi>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>m</mi>
<mn>34</mn>
</msub>
<mo>=</mo>
<mn>2</mn>
<mfrac>
<msub>
<mi>&rho;</mi>
<mn>3</mn>
</msub>
<msub>
<mi>&rho;</mi>
<mn>1</mn>
</msub>
</mfrac>
<mfrac>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>3</mn>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>3</mn>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
<mi>sin</mi>
<mi>&theta;</mi>
<mo>,</mo>
<msub>
<mi>m</mi>
<mn>44</mn>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>&rho;</mi>
<mn>3</mn>
</msub>
<msub>
<mi>&rho;</mi>
<mn>1</mn>
</msub>
</mfrac>
<mfrac>
<mrow>
<msub>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>3</mn>
</mrow>
</msub>
<msub>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
</msub>
</mrow>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mn>2</mn>
<mfrac>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>3</mn>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
</mrow>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>n</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&theta;</mi>
<mo>+</mo>
<mi>j</mi>
<mfrac>
<mrow>
<mi>&omega;</mi>
<mi>h</mi>
</mrow>
<msub>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mo>-</mo>
<mfrac>
<msub>
<mi>&rho;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>&rho;</mi>
<mn>2</mn>
</msub>
</mfrac>
<mfrac>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<mo>)</mo>
</mrow>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mn>2</mn>
<mi>&theta;</mi>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>n</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&theta;</mi>
<mo>+</mo>
<mi>j</mi>
<mfrac>
<mrow>
<mi>&omega;</mi>
<mi>h</mi>
</mrow>
<msub>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mn>2</mn>
<mfrac>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<mo>+</mo>
<mn>2</mn>
<mfrac>
<msub>
<mi>&rho;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>&rho;</mi>
<mn>2</mn>
</msub>
</mfrac>
<mfrac>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<mo>)</mo>
</mrow>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
<mo>-</mo>
<mfrac>
<msub>
<mi>&rho;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>&rho;</mi>
<mn>2</mn>
</msub>
</mfrac>
<mfrac>
<msubsup>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<mo>&rsqb;</mo>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>n</mi>
<mn>3</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mn>2</mn>
<mfrac>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>j</mi>
<mfrac>
<mrow>
<mi>&omega;</mi>
<mi>h</mi>
</mrow>
<msub>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>&rho;</mi>
<mn>2</mn>
</msub>
<msub>
<mi>&rho;</mi>
<mn>1</mn>
</msub>
</mfrac>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&theta;</mi>
<mo>+</mo>
<mfrac>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mn>2</mn>
<mi>&theta;</mi>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>n</mi>
<mn>4</mn>
</msub>
<mo>=</mo>
<mi>sin</mi>
<mn>2</mn>
<mi>&theta;</mi>
<mo>-</mo>
<mi>j</mi>
<mfrac>
<mrow>
<mi>&omega;</mi>
<mi>h</mi>
</mrow>
<msub>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mn>4</mn>
<mfrac>
<msub>
<mi>&rho;</mi>
<mn>2</mn>
</msub>
<msub>
<mi>&rho;</mi>
<mn>1</mn>
</msub>
</mfrac>
<mfrac>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<mrow>
<mo>(</mo>
<mrow>
<mfrac>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
</mrow>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
<mo>+</mo>
<mfrac>
<msub>
<mi>&rho;</mi>
<mn>2</mn>
</msub>
<msub>
<mi>&rho;</mi>
<mn>1</mn>
</msub>
</mfrac>
<mfrac>
<msubsup>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mrow>
<mn>2</mn>
<mfrac>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mrow>
<mfrac>
<msubsup>
<mi>v</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>v</mi>
<mrow>
<mi>S</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<mo>-</mo>
<mn>2</mn>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mi>sin</mi>
<mi>&theta;</mi>
<mo>,</mo>
</mrow>
vPi、vSi、ρi分别为各层的纵波速度、横波速度及密度,i=1,2,3,分别对应单薄层的三层介质、θ为入射角、h为薄层厚度,ω为圆频率,ω=2πf,其中f为入射波的频率,
3.根据权利要求2所述的薄层PP波反射系数的快速求解方法,其特征在于,所述系数矩阵方程满足如下公式:
M[R'PP -R'PS T′PP -T′PS]T=N,
其中,R'PP、R'PS、T′PP、T′PS分别为入射角取负值情况下极薄层的反透射系数,
所述系数矩阵方程与所述极薄层反射透射矩阵方程对比满足如下公式:
R'PP=RPP、R'PS=-RPS、T′PP=TPP、T′PS=-TPS,
所述反透射系数与入射角正弦函数的奇偶关系式满足如下公式:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>P</mi>
<mi>P</mi>
</mrow>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>&infin;</mi>
</munderover>
<msub>
<mi>A</mi>
<mrow>
<mn>2</mn>
<mi>n</mi>
</mrow>
</msub>
<msup>
<mi>sin</mi>
<mrow>
<mn>2</mn>
<mi>n</mi>
</mrow>
</msup>
<mi>&theta;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>P</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>&infin;</mi>
</munderover>
<msub>
<mi>A</mi>
<mrow>
<mn>2</mn>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<msup>
<mi>sin</mi>
<mrow>
<mn>2</mn>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msup>
<mi>&theta;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>T</mi>
<mrow>
<mi>P</mi>
<mi>P</mi>
</mrow>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>&infin;</mi>
</munderover>
<msub>
<mi>B</mi>
<mrow>
<mn>2</mn>
<mi>n</mi>
</mrow>
</msub>
<msup>
<mi>sin</mi>
<mrow>
<mn>2</mn>
<mi>n</mi>
</mrow>
</msup>
<mi>&theta;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>T</mi>
<mrow>
<mi>P</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>&infin;</mi>
</munderover>
<msub>
<mi>B</mi>
<mrow>
<mn>2</mn>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<msup>
<mi>sin</mi>
<mrow>
<mn>2</mn>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msup>
<mi>&theta;</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
</mrow>
其中A2n、A2n+1、B2n、B2n+1分别为反透射系数幂级数表达式的系数。
4.根据权利要求3所述的薄层PP波反射系数的快速求解方法,其特征在于,所述入射角正弦函数幂级数的系数满足如下公式:
<mrow>
<msub>
<mi>m</mi>
<mn>0</mn>
</msub>
<msub>
<mi>A</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>z</mi>
<mrow>
<mi>P</mi>
<mn>3</mn>
</mrow>
</msub>
<msub>
<mi>z</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>j</mi>
<mi>&tau;</mi>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>z</mi>
<mrow>
<mi>P</mi>
<mn>2</mn>
</mrow>
</msub>
<msub>
<mi>z</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mo>-</mo>
<mfrac>
<msub>
<mi>z</mi>
<mrow>
<mi>P</mi>
<mn>3</mn>
</mrow>
</msub>
<msub>
<mi>z</mi>
<mrow>
<mi>P</mi>
<mn>2</mn>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
m0B0=2(1+τ2),
其中,为系数A0、B0的分母,τ=ωh/vP2为薄层内纵波的单程相移,zP=ρvP为纵波阻抗,它们的下角标取值为1、2、3,分别对应单薄层的三层介质;
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>m</mi>
<mn>1</mn>
</msub>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>b</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<mfrac>
<msub>
<mi>z</mi>
<mrow>
<mi>S</mi>
<mn>3</mn>
</mrow>
</msub>
<msub>
<mi>z</mi>
<mrow>
<mi>S</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mo>-</mo>
<msub>
<mi>j&tau;r</mi>
<mn>2</mn>
</msub>
<mfrac>
<msub>
<mi>z</mi>
<mrow>
<mi>S</mi>
<mn>2</mn>
</mrow>
</msub>
<msub>
<mi>z</mi>
<mrow>
<mi>S</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<msub>
<mi>A</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>l</mi>
<mn>1</mn>
</msub>
<msub>
<mi>l</mi>
<mn>2</mn>
</msub>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mn>2</mn>
<msub>
<mi>r</mi>
<mn>3</mn>
</msub>
</mfrac>
</mrow>
<mo>)</mo>
<mfrac>
<msub>
<mi>z</mi>
<mrow>
<mi>s</mi>
<mn>3</mn>
</mrow>
</msub>
<msub>
<mi>z</mi>
<mrow>
<mi>S</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<msub>
<mi>B</mi>
<mn>0</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mo>&lsqb;</mo>
<mfrac>
<mn>2</mn>
<msub>
<mi>r</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>+</mo>
<mi>j</mi>
<mi>&tau;</mi>
<mo>(</mo>
<mrow>
<mo>-</mo>
<msub>
<mi>l</mi>
<mn>1</mn>
</msub>
<mfrac>
<msub>
<mi>z</mi>
<mrow>
<mi>S</mi>
<mn>3</mn>
</mrow>
</msub>
<msub>
<mi>z</mi>
<mrow>
<mi>S</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mo>+</mo>
<mn>2</mn>
<mfrac>
<msub>
<mi>r</mi>
<mn>2</mn>
</msub>
<msub>
<mi>r</mi>
<mn>1</mn>
</msub>
</mfrac>
<mfrac>
<msub>
<mi>z</mi>
<mrow>
<mi>S</mi>
<mn>3</mn>
</mrow>
</msub>
<msub>
<mi>z</mi>
<mrow>
<mi>S</mi>
<mn>2</mn>
</mrow>
</msub>
</mfrac>
</mrow>
<mo>)</mo>
<mo>&rsqb;</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>A</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>,</mo>
</mrow>
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>m</mi>
<mn>1</mn>
</msub>
<msub>
<mi>B</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mo>&lsqb;</mo>
<msub>
<mi>b</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<msub>
<mi>j&tau;r</mi>
<mn>2</mn>
</msub>
<mfrac>
<msub>
<mi>z</mi>
<mrow>
<mi>S</mi>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>z</mi>
<mrow>
<mi>S</mi>
<mn>2</mn>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<msub>
<mi>A</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mo>&lsqb;</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msub>
<mi>b</mi>
<mn>2</mn>
</msub>
</mrow>
<msub>
<mi>r</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>+</mo>
<msub>
<mi>j&tau;l</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>+</mo>
<msub>
<mi>j&tau;r</mi>
<mn>2</mn>
</msub>
<mfrac>
<msub>
<mi>z</mi>
<mrow>
<mi>S</mi>
<mn>2</mn>
</mrow>
</msub>
<msub>
<mi>z</mi>
<mrow>
<mi>S</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
</mrow>
<mo>)</mo>
<mo>&rsqb;</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>A</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>l</mi>
<mn>1</mn>
</msub>
<msub>
<mi>l</mi>
<mn>2</mn>
</msub>
<mo>&lsqb;</mo>
<mn>1</mn>
<mo>+</mo>
<mfrac>
<mn>2</mn>
<msub>
<mi>r</mi>
<mn>3</mn>
</msub>
</mfrac>
<mfrac>
<msub>
<mi>z</mi>
<mrow>
<mi>S</mi>
<mn>3</mn>
</mrow>
</msub>
<msub>
<mi>z</mi>
<mrow>
<mi>S</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mo>+</mo>
<msub>
<mi>j&tau;r</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mfrac>
<msub>
<mi>z</mi>
<mrow>
<mi>S</mi>
<mn>2</mn>
</mrow>
</msub>
<msub>
<mi>z</mi>
<mrow>
<mi>S</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mo>+</mo>
<mfrac>
<mn>2</mn>
<msub>
<mi>r</mi>
<mn>3</mn>
</msub>
</mfrac>
<mfrac>
<msub>
<mi>z</mi>
<mrow>
<mi>S</mi>
<mn>3</mn>
</mrow>
</msub>
<msub>
<mi>z</mi>
<mrow>
<mi>S</mi>
<mn>2</mn>
</mrow>
</msub>
</mfrac>
</mrow>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<msub>
<mi>B</mi>
<mn>0</mn>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>,</mo>
</mrow>
其中,zS=ρvS为横波阻抗,r=vP/vS为纵横波速度比,它们的下角标取值为1、2、3,分别对应薄层的三层介质,lk=vPk+1/vPk为相邻层纵波速度的比值,其中k取值为1或2,为系数A1、B1的分母,和为简化替换参量;
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>m</mi>
<mn>0</mn>
</msub>
<msub>
<mi>A</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mrow>
<mo>(</mo>
<mrow>
<mfrac>
<msub>
<mi>z</mi>
<mrow>
<mi>P</mi>
<mn>3</mn>
</mrow>
</msub>
<msub>
<mi>z</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mo>+</mo>
<mi>j</mi>
<mi>&tau;</mi>
<mfrac>
<msub>
<mi>z</mi>
<mrow>
<mi>P</mi>
<mn>2</mn>
</mrow>
</msub>
<msub>
<mi>z</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mn>4</mn>
<msub>
<mi>j&tau;l</mi>
<mn>1</mn>
</msub>
</mrow>
<msubsup>
<mi>r</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
</mfrac>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>A</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<mfrac>
<mn>2</mn>
<msub>
<mi>r</mi>
<mn>1</mn>
</msub>
</mfrac>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mo>&lsqb;</mo>
<mfrac>
<mn>2</mn>
<msubsup>
<mi>r</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
</mfrac>
<mo>+</mo>
<mi>j</mi>
<mi>&tau;</mi>
<mo>(</mo>
<mrow>
<msub>
<mi>l</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msub>
<mi>l</mi>
<mn>1</mn>
</msub>
</mrow>
<msubsup>
<mi>r</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
</mfrac>
<mo>+</mo>
<mfrac>
<mn>2</mn>
<msubsup>
<mi>r</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
</mfrac>
<mfrac>
<msub>
<mi>z</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>z</mi>
<mrow>
<mi>P</mi>
<mn>2</mn>
</mrow>
</msub>
</mfrac>
</mrow>
<mo>)</mo>
<mfrac>
<msub>
<mi>z</mi>
<mrow>
<mi>P</mi>
<mn>3</mn>
</mrow>
</msub>
<msub>
<mi>z</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mo>&rsqb;</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<msub>
<mi>A</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>j&tau;l</mi>
<mn>1</mn>
</msub>
</mrow>
<msub>
<mi>r</mi>
<mn>1</mn>
</msub>
</mfrac>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mn>4</mn>
<msubsup>
<mi>r</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
</mfrac>
</mrow>
<mo>)</mo>
</mrow>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msubsup>
<mi>l</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>l</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mrow>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mo>-</mo>
<mfrac>
<mn>2</mn>
<msubsup>
<mi>r</mi>
<mn>3</mn>
<mn>2</mn>
</msubsup>
</mfrac>
</mrow>
<mo>)</mo>
</mrow>
<mfrac>
<msub>
<mi>z</mi>
<mrow>
<mi>P</mi>
<mn>3</mn>
</mrow>
</msub>
<msub>
<mi>z</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<msub>
<mi>B</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>l</mi>
<mn>1</mn>
</msub>
<msub>
<mi>l</mi>
<mn>2</mn>
</msub>
</mrow>
<msub>
<mi>r</mi>
<mn>3</mn>
</msub>
</mfrac>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mn>2</mn>
<msub>
<mi>r</mi>
<mn>3</mn>
</msub>
</mfrac>
</mrow>
<mo>)</mo>
</mrow>
<mfrac>
<msub>
<mi>z</mi>
<mrow>
<mi>P</mi>
<mn>3</mn>
</mrow>
</msub>
<msub>
<mi>z</mi>
<mrow>
<mi>P</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<msub>
<mi>B</mi>
<mn>1</mn>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>,</mo>
</mrow>
其中,A0、B0、A1、B1、A2为所述入射角正弦函数幂级数的系数。
5.根据权利要求4所述的薄层PP波反射系数的快速求解方法,其特征在于,所述薄层PP波反射系数的解析解近似式满足如下公式:
RPP=A0+A2sin2θ,
其中,RPP为所述薄层PP波反射系数的解析解近似式。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711009284.0A CN107861153A (zh) | 2017-10-25 | 2017-10-25 | 薄层pp波反射系数的快速求解方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711009284.0A CN107861153A (zh) | 2017-10-25 | 2017-10-25 | 薄层pp波反射系数的快速求解方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107861153A true CN107861153A (zh) | 2018-03-30 |
Family
ID=61697678
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711009284.0A Pending CN107861153A (zh) | 2017-10-25 | 2017-10-25 | 薄层pp波反射系数的快速求解方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107861153A (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109188511A (zh) * | 2018-08-27 | 2019-01-11 | 中国地质大学(北京) | 一种砂泥岩薄互层介质多波avo联合反演方法 |
CN109613608A (zh) * | 2018-12-10 | 2019-04-12 | 中国石油大学(华东) | 一种任意复杂介质地震波传播矩阵模拟方法 |
CN110542924A (zh) * | 2019-09-02 | 2019-12-06 | 成都理工大学 | 一种高精度的纵横波阻抗反演方法 |
CN110673212A (zh) * | 2019-10-25 | 2020-01-10 | 北京多分量地震技术研究院 | 一种模型约束的薄层多波ava联合反演方法 |
CN110764145A (zh) * | 2019-10-10 | 2020-02-07 | 淮南矿业(集团)有限责任公司 | 薄层顶底界面反射系数的反演方法及装置 |
CN112505751A (zh) * | 2019-09-16 | 2021-03-16 | 中国石油化工股份有限公司 | 球面波ps反射系数计算方法及系统 |
CN112649845A (zh) * | 2019-10-10 | 2021-04-13 | 中国石油化工股份有限公司 | 球面波反射系数计算方法及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140297188A1 (en) * | 2013-03-29 | 2014-10-02 | Cgg Services Sa | Time-frequency representations of seismic traces using wigner-ville distributions |
CN104965224A (zh) * | 2015-06-03 | 2015-10-07 | 北京多分量地震技术研究院 | 用平均入射角道集进行pp波与ps波联合avo反演方法 |
CN105629301A (zh) * | 2016-03-29 | 2016-06-01 | 中国地质大学(北京) | 薄层弹性波反射系数快速求解方法 |
CN106556862A (zh) * | 2015-09-29 | 2017-04-05 | 中国石油化工股份有限公司 | 用于页岩气勘探avo技术的pp波反射系数计算方法 |
-
2017
- 2017-10-25 CN CN201711009284.0A patent/CN107861153A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140297188A1 (en) * | 2013-03-29 | 2014-10-02 | Cgg Services Sa | Time-frequency representations of seismic traces using wigner-ville distributions |
CN104965224A (zh) * | 2015-06-03 | 2015-10-07 | 北京多分量地震技术研究院 | 用平均入射角道集进行pp波与ps波联合avo反演方法 |
CN106556862A (zh) * | 2015-09-29 | 2017-04-05 | 中国石油化工股份有限公司 | 用于页岩气勘探avo技术的pp波反射系数计算方法 |
CN105629301A (zh) * | 2016-03-29 | 2016-06-01 | 中国地质大学(北京) | 薄层弹性波反射系数快速求解方法 |
Non-Patent Citations (1)
Title |
---|
石瑛,等: "平均入射角道集PP波与PS波联合反演", 《地球物理学报》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109188511A (zh) * | 2018-08-27 | 2019-01-11 | 中国地质大学(北京) | 一种砂泥岩薄互层介质多波avo联合反演方法 |
CN109613608A (zh) * | 2018-12-10 | 2019-04-12 | 中国石油大学(华东) | 一种任意复杂介质地震波传播矩阵模拟方法 |
CN110542924A (zh) * | 2019-09-02 | 2019-12-06 | 成都理工大学 | 一种高精度的纵横波阻抗反演方法 |
CN112505751A (zh) * | 2019-09-16 | 2021-03-16 | 中国石油化工股份有限公司 | 球面波ps反射系数计算方法及系统 |
CN110764145A (zh) * | 2019-10-10 | 2020-02-07 | 淮南矿业(集团)有限责任公司 | 薄层顶底界面反射系数的反演方法及装置 |
CN112649845A (zh) * | 2019-10-10 | 2021-04-13 | 中国石油化工股份有限公司 | 球面波反射系数计算方法及系统 |
CN110764145B (zh) * | 2019-10-10 | 2021-07-23 | 淮南矿业(集团)有限责任公司 | 薄层顶底界面反射系数的反演方法及装置 |
CN110673212A (zh) * | 2019-10-25 | 2020-01-10 | 北京多分量地震技术研究院 | 一种模型约束的薄层多波ava联合反演方法 |
CN110673212B (zh) * | 2019-10-25 | 2021-06-18 | 北京多分量地震技术研究院 | 一种模型约束的薄层多波ava联合反演方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107861153A (zh) | 薄层pp波反射系数的快速求解方法 | |
Chaljub et al. | Quantitative comparison of four numerical predictions of 3D ground motion in the Grenoble Valley, France | |
CN103389513B (zh) | 应用声波测井资料约束反演提高地震资料分辨率的方法 | |
Ladopoulos | Elastodynamics for Non-linear Seismic Wave Motion in Real-Time Expert Seismology | |
Zhang et al. | Robust streamline tracing using inter‐cell fluxes in locally refined and unstructured grids | |
CN103412333B (zh) | 一种静校正基准面确定方法 | |
CN103439739B (zh) | 地球物理勘探用纵横波匹配方法及匹配装置 | |
CN104950332B (zh) | 一种弹性多层介质中平面波反射系数的计算方法 | |
CN109444959B (zh) | 全频高精度层速度场建立方法 | |
CN105629301A (zh) | 薄层弹性波反射系数快速求解方法 | |
CN109470187A (zh) | 基于地震三属性的储层厚度预测方法 | |
CN107193046A (zh) | 一种基于谱反演的砂体厚度预测方法及系统 | |
CN104820242B (zh) | 一种面向叠前反演的道集振幅分频补偿方法 | |
CN106199704B (zh) | 一种三维三分量海底电缆地震资料速度建模方法 | |
CN104237937A (zh) | 叠前地震反演方法及其系统 | |
CN111722283B (zh) | 一种地层速度模型建立方法 | |
CN106569260B (zh) | 针对高陡构造地区二维地震变速空校成图方法 | |
CN105093279A (zh) | 针对山前带的三维地震初至波菲涅尔体层析反演方法 | |
CN106547020A (zh) | 一种地震数据的保幅处理方法 | |
CN102768367B (zh) | 基于三重条件约束的双相介质avo正演方法 | |
CN103984018B (zh) | 多波联合振幅随入射角变化反演方法 | |
CN108279437A (zh) | 变密度声波方程时间高阶精度交错网格有限差分方法 | |
CN113156505B (zh) | 一种“基底构造-古地貌-地震相”三要素递进约束的断陷湖盆礁滩储层识别方法 | |
CN104007465B (zh) | 基于弹性反演的纯纵波数据刻画储层发育细节方法 | |
CN104678439A (zh) | 一种多波连续介质速度分析方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180330 |
|
RJ01 | Rejection of invention patent application after publication |