CN104570085B - 一种纵横波射线参数域联合反演方法 - Google Patents
一种纵横波射线参数域联合反演方法 Download PDFInfo
- Publication number
- CN104570085B CN104570085B CN201310520644.9A CN201310520644A CN104570085B CN 104570085 B CN104570085 B CN 104570085B CN 201310520644 A CN201310520644 A CN 201310520644A CN 104570085 B CN104570085 B CN 104570085B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mfrac
- msup
- wave
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 34
- 239000012530 fluid Substances 0.000 claims abstract description 12
- 238000003384 imaging method Methods 0.000 claims description 11
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000004422 calculation algorithm Methods 0.000 claims description 2
- 239000003255 cyclooxygenase 2 inhibitor Substances 0.000 claims description 2
- 230000015572 biosynthetic process Effects 0.000 abstract 2
- 239000003209 petroleum derivative Substances 0.000 abstract 1
- 230000000704 physical effect Effects 0.000 abstract 1
- 238000005516 engineering process Methods 0.000 description 6
- 238000011160 research Methods 0.000 description 5
- 238000011161 development Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000005553 drilling Methods 0.000 description 3
- 238000001514 detection method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 235000019738 Limestone Nutrition 0.000 description 1
- 238000005452 bending Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000006028 limestone Substances 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 239000003607 modifier Substances 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 239000003129 oil well Substances 0.000 description 1
- 239000009671 shengli Substances 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种纵横波射线参数域联合反演方法，通过将纵波、转换波叠前道集资料转换到射线参数域，利用地震振幅随射线参数变化的信息反演得到地层弹性参数，进而对岩性和流体进行预测。本发明是同点物性联合反演，纵波和转换波来自于同一条入射路径，虽然反射路径彼此分离，但是来自于同一深度同一反射点，各参数之间的联系更加紧密，因此能得到更可靠的地层弹性参数，提高利用纵横波地震资料进行油气预测的精度。
Description
技术领域
本发明属于地震资料解释领域，是一种纵波、转换横波地震资料的联合反演方法。
背景技术
在油气检测和油藏描述中，介质的弹性参数十分重要，这些参数与岩性和流体成分有关。描述弹性介质的最主要的3个参数是纵波速度、横波速度和密度。目前地球物理勘探的常规技术是利用纵波AVO(Amplitude Versus Offset，振幅随炮检距变化)叠前道集反演上述3个参数，反演时考虑了纵波反射振幅随入射角变化的因素。但是大量理论研究和实例应用表明，在有限偏移距的纵波资料中，这3个参数的准确求解存在一定的多解性。随着多波多分量地震勘探技术的不断发展，高质量转换波资料的获取，使得纵波、转换波联合反演技术成为随之发展起来的新技术。利用纵波和转换横波地震数据进行联合反演，可以减少反演结果的多解性，提高反演参数的精度，提供更多油藏描述和监测的潜力。随着近年来多波勘探采集和处理技术的发展，纵波、转换波联合反演的研究得到了较大的发展，出现了多种联合反演方法，通常是基于Zeoppritz方程或其线性近似进行反演，这些方法均利用纵波、转换波地震振幅随入射角的变化信息反演得到地层的弹性参数。纵波、转换波联合反演方法的基本原理是基于Smith和Gidlow(1987)提出的加权叠加方法，通过PP数据的加权叠加来计算纵波和横波波阻抗以及其它的相关参数。Fatti等(1994)做了进一步的研究，取消了对于Gardner公式的依赖并且包含了密度项。Stewart等(2001)、Landro等(2006)将加权叠加的思想用于联合反演，给出了实用的纵横波联合反演方法，反演结果得到的是纵波速度、横波速度和密度的变化率ΔVp/Vp、ΔVs/Vs、Δρ/ρ。但这些方法仅使用纵波、转换波地震振幅随入射角的变化信息，没有考虑纵波、转换波地震振幅随射线参数的变化信息，所反演出的结果精度较低，与钻井、测井、试油、生产数据进行对比和统计分析的结果吻合率不高。
纵波、转换波地震振幅随射线参数的变化信息是多波地震勘探的重要参数，用纵波、转换波地震振幅随射线参数的变化信息进行联合反演方面目前还缺乏类似的研究。由于纵波、转换波地震振幅随射线参数的变化信息中包含着更准确的地震能量传递信息，因此纵横波射线参数域联合反演方法的研究和应用对于更充分地利用纵横波地震信息进行油气检测有着重要意义。
发明内容
针对现有技术的不足，本发明的目的在于提供一种利用纵波、转换波地震振幅随射线参数的变化信息进行联合反演的方法。
本发明采用的技术方案是：一种纵横波射线参数域联合反演方法，其特征在于：所述方法通过将纵波、转换波叠前道集资料转换到射线参数域，利用地震振幅随射线参数变化的信息反演得到地层弹性参数，进而对岩性和流体进行预测。
具体包括以下步骤：
步骤一：采用射线追踪法将PP波、PS波共成像点道集数据从偏移距域转换到射线参数域，并进行分射线参数叠加，得到分射线参数叠加剖面；
射线追踪法是指给定发射点和接收点位置及介质的波速，求从发射点到接收点的射线轨迹。射线追踪方法有很多种，本发明采用弯曲射线追踪法，根据最小走时准则迭代求出接收点处的射线路径，然后将偏移距域共成像点道集中某一确定时间和偏移距的反射同相轴，放置到相应射线参数值的相同时间处，即完成共成像点道集的射线参数域转换。利用得到的射线参数域共成像点道集进行分射线参数叠加，即得到用于反演的分射线参数叠加剖面。
步骤二：分别对步骤一得到的纵波、转换波分射线参数叠加剖面进行稀疏约束反射系数序列反演，得到PP波、PS波分射线参数的反射系数剖面。
稀疏约束反射系数序列反演的方法原理是：
地震反射系数序列反演与地震数据处理中提高数据分辨率的反褶积紧密关联，而反褶积的最佳结果正是地层的反射系数序列。从统计学角度认为反射系数序列呈高斯分布，但是传统的反褶积结果通常较光滑，很难作为反射系数序列来进一步使用。实际地震反射系数序列通常是脉冲状的，它们比高斯分布更稀疏，而与柯西分布更接近。为此，本方法中以迭代反演代替反褶积处理，在反演中引入柯西准则，实现柯西稀疏约束的反射系数反演，提高解的稀疏性，得到的反射系数序列具有较好的高频带特性。
柯西约束反射系数序列反演的目标函数：
式中，r：反射系数，d：地震数据，Cw：子波矩阵，CD：地震数据协方差矩阵，λ：柯西分布参数
令得到的解为
这里P是一个斜对角阵，定义为
解(2)式通过迭代方式得到，而柯西分布参数λ则控制了反射系数向量r的稀疏程度。
步骤三：利用步骤二得到的PP波和PS波的分射线参数反射系数剖面进行联合反演，得到纵波速度、横波速度及密度3个参数。
联合反演的方法：利用PP波、PS波射线参数域高阶反射系数的线性近似公式，采用广义线性反演算法同步估算纵波速度、横波速度及密度3个弹性参数。
PP波射线参数域的高阶反射系数：
Rpp(p)≈Rf-A1p2+B1p4 (3)
其中，
PP波射线参数域的线性近似反射系数：
PS波射线参数域的高阶反射系数：
PS波射线参数域的线性近似反射系数：
P是地下第i层介质的射线参数，θi和分别代表了入射角和反射角；αi,βi和ρi分别代表了纵波速度、横波速度和密度，α、β和ρ代表了上述三个参数在相邻介质的平均值。纵波相关的反射系数Rf可表示为其中和是纵波垂直慢度，qα和qβ分别代表了纵波和横波的平均垂直慢度。Δμ＝ρi+1βi+1 2-ρiβi 2是剪切模量差。
广义线性反演方法：通过用泰勒展开把非线性反问题化成线性反问题进行反演求解。将PP、PS波线性近似反射系数公式对弹性参数模型求一阶泰勒展开式，可得
其中RPP/ps是由步骤二求得的实际数据的射线参数域反射系数值，RPP/ps(m)是由PP、PS波线性近似反射系数公式求得的理论反射系数值，m0＝[α0,β0,ρ0]是初始模型向量，m＝[α,β,ρ]是待解的模型向量。
通过迭代计算，不断求解(5)式并修改模型参数，通过一定的准则(本方法采用模型修改量的模达到极小)得到模型的逼近值，即最终的反演结果。
步骤四：利用步骤三得到的纵波速度、横波速度及密度三个参数，计算纵横波速度比、泊松比等对岩性和流体敏感的参数，对储层及流体分布进行预测。
射线追踪的理论基础是,在高频近似条件下,地震波场的主能量沿射线轨迹传播。由于地震波的振幅主要和地震波场的主能量有关，因此本发明利用地震振幅随射线参数变化的信息反演，相比现有技术中利用纵波、转换波地震振幅随入射角的变化信息反演，得到的地层弹性参数更加准确，能有效的提高反演精度。
由于采用的是PP波、PS波共成像点道集数据，因此纵波和转换波来自于同一条入射路径，虽然反射路径彼此分离，但是来自于同一深度同一反射点。因此，射线参数域反演是同点物性联合反演。采用这种同点物性联合反演方法，各参数之间的联系更加紧密，因此能更准确地对纵波、转换波资料进行联合反演，得到更可靠的地层弹性参数，提高纵横波地震资料油气预测的精度。
本发明实现了纵波、转换波射线参数域的联合反演，反演精度较高，为从多波地震信息中提取岩性和流体参数提供了一种新的解决途径。
附图说明
图1罗家地区PP波偏移距域共成像点道集图。
图2罗家地区PP波射线参数域共成像点道集图。
图3罗家地区PS波偏移距域共成像点道集图。
图4罗家地区PS波射线参数域共成像点道集图。
图5罗家地区射线参数值为60的PP波地震剖面图。
图6罗家地区射线参数值为60的PS波地震剖面图。
图7罗家地区PP波射线参数值为60的反射系数剖面图。
图8罗家地区PS波射线参数值为60的反射系数剖面图。
图9罗家地区射线参数域联合反演得到的纵波速度剖面图。
图10罗家地区射线参数域联合反演得到的横波速度剖面图。
图11罗家地区射线参数域联合反演得到的密度剖面图。
图12罗家地区射线参数域联合反演纵横波速度比平面图。
具体实施方式
现在将在下文中详细描述本发明的示例性实施例；然而，本发明的示例性实施例可以以许多不同的形式来实施，且不应该限于这里阐述的示例。相反，提供这些示例使得本公开将是彻底并完整的，并将向本领域技术人员充分地传达本发明的范围。
根据本发明的示例性实施例，可以预先对采集到的地震资料进行处理，然后将处理后的地震数据输入到例如计算机中。例如，可以以本领域公知的方法进行二维多分量或三维多分量野外勘探采集，以得到原始采集的多分量数据。然后，经过二维或三维纵波、转换波地震资料处理，以得到包括能够进行储层纵横波速度比反演的纵波、转换波叠前时间偏移道集数据的地震数据。然后，可以将得到的地震数据输入到计算机中，从而在后面的步骤中对这样的地震数据进行处理。可选择地，也可以直接将采集到的地震资料输入到计算机中并在计算机中进行相应地处理。如此，输入的地震数据可以包括与二维纵波和转换波相关的数据和与三维纵波和转换波相关的数据中的至少一种。
本发明的方法在胜利油田罗家多波工区进行了实验，取得了较好的应用效果。本实验利用研究区三维纵波、转换波地震资料进行射线参数域叠前联合反演，并对目的层进行含油气性预测。
首先是射线追踪及分射线参数剖面的求取。如图1到图4所示，采用射线追踪方法将纵波、转换波共成像点道集从偏移距域转换到射线参数域，并进行分射线参数叠加，其射线参数从0到140ms/km，间隔为20ms/km。图5、图6为得到的射线参数值为60ms/km的PP、PS波地震剖面。
然后分别对纵波、转换波分射线参数叠加剖面进行稀疏约束反射系数序列反演，反射系数反演的主要目的是去除子波的影响，图7、图8为得到的PP波、PS波分射线参数的反射系数剖面，可以看到，稀疏约束反演结果反映了更为精细的地层特征。
利用得到的PP波和PS波的分射线参数反射系数剖面，运用三参数同步反演方法完成联合反演。图9到图11为得到的纵波速度、横波速度及密度剖面图。
利用联合反演得到的纵波速度、横波速度等参数，计算速度比、泊松比等对岩性和流体敏感的参数，对储层及流体分布进行预测。图12为罗家地区目的层的射线参数域联合反演纵横波速度比平面图，图中用不同的井符号代表实际井的含油气情况(其中黑色圆点代表油井，白色圆点代表水井，圆圈代表干井及无生物灰岩井)，对比纵横波速度比参数可以看出，大部分钻遇油层的井位于纵横波速度比低值区(红、黄、绿色区)，大部分钻遇水层和未钻遇该套储层的井位于纵横波速度比高值区(蓝色区)；纵横波速度比的分布规律性较好，低值区主要分布于罗37至罗9-1井一带，与钻井揭示的油气发育区对应。利用区内56口井的钻井、测井、试油、生产数据进行对比和统计分析，吻合率达到85.7％，因此，反演取得了良好的效果。
Claims (1)
1.一种纵横波射线参数域联合反演方法，其特征在于：所述方法通过将纵波、转换波叠前道集资料转换到射线参数域，利用地震振幅随射线参数变化的信息反演得到地层弹性参数，进而对岩性和流体进行预测；
所述方法包括以下步骤：
步骤一：采用射线追踪法将PP波、PS波共成像点道集数据从偏移距域转换到射线参数域，并进行分射线参数叠加，得到分射线参数叠加剖面；
步骤二：分别对步骤一得到的纵波、转换波分射线参数叠加剖面进行稀疏约束反射系数序列反演，得到PP波、PS波分射线参数的反射系数剖面；
步骤三：利用步骤二得到的PP波和PS波的分射线参数反射系数剖面进行联合反演，得到纵波速度、横波速度及密度3个参数；
步骤四：利用步骤三得到的纵波速度、横波速度及密度三个参数，进一步计算得到纵横波速度比、泊松比等对岩性和流体敏感的参数，对储层及流体分布进行预测；
步骤二中所述稀疏约束反射系数序列反演的方法是：以迭代反演代替反褶积处理，在反演中引入柯西准则，实现柯西稀疏约束的反射系数反演，提高解的稀疏性，得到的反射系数序列具有较好的高频带特性；
柯西约束反射系数序列反演的目标函数：
<mrow>
<mi>J</mi>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mi>d</mi>
<mo>-</mo>
<msub>
<mi>C</mi>
<mi>w</mi>
</msub>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<msubsup>
<mi>C</mi>
<mi>D</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>d</mi>
<mo>-</mo>
<msub>
<mi>C</mi>
<mi>w</mi>
</msub>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>N</mi>
<mi> </mi>
<mi>ln</mi>
<mrow>
<mo>(</mo>
<mi>&pi;</mi>
<mi>&lambda;</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<mi>l</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mfrac>
<msubsup>
<mi>r</mi>
<mi>i</mi>
<mn>2</mn>
</msubsup>
<msup>
<mi>&lambda;</mi>
<mn>2</mn>
</msup>
</mfrac>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
式中，r：反射系数，d：地震数据，Cw：子波矩阵，CD：地震数据协方差矩阵，λ：柯西分布参数
令得到的解为
<mrow>
<mi>r</mi>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>C</mi>
<mi>w</mi>
<mi>T</mi>
</msubsup>
<msubsup>
<mi>C</mi>
<mi>D</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<msub>
<mi>C</mi>
<mi>w</mi>
</msub>
<mo>+</mo>
<mfrac>
<mn>2</mn>
<msup>
<mi>&lambda;</mi>
<mn>2</mn>
</msup>
</mfrac>
<mi>P</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msubsup>
<mi>C</mi>
<mi>w</mi>
<mi>T</mi>
</msubsup>
<msubsup>
<mi>C</mi>
<mi>D</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mi>d</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
这里P是一个斜对角阵，定义为
<mrow>
<mi>P</mi>
<mo>=</mo>
<mi>d</mi>
<mi>i</mi>
<mi>a</mi>
<mi>g</mi>
<mo>{</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msubsup>
<mi>r</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
<msup>
<mi>&lambda;</mi>
<mn>2</mn>
</msup>
</mfrac>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>,</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msubsup>
<mi>r</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<msup>
<mi>&lambda;</mi>
<mn>2</mn>
</msup>
</mfrac>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msubsup>
<mi>r</mi>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
<msup>
<mi>&lambda;</mi>
<mn>2</mn>
</msup>
</mfrac>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>}</mo>
</mrow>
解(2)式通过迭代方式得到，而柯西分布参数λ则控制了反射系数向量r的稀疏程度；
步骤三中所述联合反演利用PP波、PS波射线参数域高阶反射系数的线性近似公式，采用广义线性反演算法同步估算纵波速度、横波速度及密度3个弹性参数；
所述PP波射线参数域的高阶反射系数为：
Rpp(p)≈Rf-A1p2+B1p4 (3)
其中，
<mrow>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mn>2</mn>
<mfrac>
<mrow>
<mi>&Delta;</mi>
<mi>&mu;</mi>
</mrow>
<mi>&rho;</mi>
</mfrac>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>R</mi>
<mi>f</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>q</mi>
<mi>&alpha;</mi>
</msub>
<msub>
<mi>q</mi>
<mi>&beta;</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>&Delta;</mi>
<mi>&mu;</mi>
</mrow>
<mi>&rho;</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<msub>
<mi>B</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mn>2</mn>
<msub>
<mi>R</mi>
<mi>f</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>&Delta;</mi>
<mi>&mu;</mi>
</mrow>
<mi>&rho;</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<mn>2</mn>
<msub>
<mi>q</mi>
<mi>&alpha;</mi>
</msub>
<msub>
<mi>q</mi>
<mi>&beta;</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>&Delta;</mi>
<mi>&mu;</mi>
</mrow>
<mi>&rho;</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>R</mi>
<mi>f</mi>
</msub>
<mo>)</mo>
</mrow>
<msubsup>
<mi>q</mi>
<mi>&alpha;</mi>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>q</mi>
<mi>&beta;</mi>
<mn>2</mn>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>&Delta;</mi>
<mi>&mu;</mi>
</mrow>
<mi>&rho;</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mn>4</mn>
</msup>
</mrow>
所述PP波射线参数域的线性近似反射系数为：
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>p</mi>
<mi>p</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>p</mi>
<mo>)</mo>
</mrow>
<mo>&ap;</mo>
<msub>
<mi>R</mi>
<mi>f</mi>
</msub>
<mo>-</mo>
<mn>2</mn>
<mfrac>
<mrow>
<mi>&Delta;</mi>
<mi>&mu;</mi>
</mrow>
<mi>&rho;</mi>
</mfrac>
<msup>
<mi>p</mi>
<mn>2</mn>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
所述PS波射线参数域的高阶反射系数为：
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>p</mi>
<mi>s</mi>
</mrow>
</msub>
<mo>&ap;</mo>
<mo>-</mo>
<mfrac>
<msub>
<mi>&alpha;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>&beta;</mi>
<mn>1</mn>
</msub>
</mfrac>
<mfrac>
<mrow>
<mi>&Delta;</mi>
<mi>&mu;</mi>
</mrow>
<mi>&rho;</mi>
</mfrac>
<msub>
<mi>q</mi>
<mi>&alpha;</mi>
</msub>
<mi>p</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>A</mi>
<mn>2</mn>
</msub>
<msup>
<mi>p</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>B</mi>
<mn>2</mn>
</msub>
<msup>
<mi>p</mi>
<mn>4</mn>
</msup>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，
<mrow>
<msub>
<mi>A</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msub>
<mi>q</mi>
<mi>&alpha;</mi>
</msub>
<msub>
<mi>q</mi>
<mi>&beta;</mi>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mn>2</mn>
<mfrac>
<mrow>
<mi>&Delta;</mi>
<mi>&mu;</mi>
</mrow>
<mi>&rho;</mi>
</mfrac>
<mo>+</mo>
<msub>
<mi>q</mi>
<mi>&alpha;</mi>
</msub>
<msub>
<mi>q</mi>
<mi>&beta;</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>&Delta;</mi>
<mi>&mu;</mi>
</mrow>
<mi>&rho;</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<msub>
<mi>B</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mfrac>
<mn>2</mn>
<mrow>
<msub>
<mi>q</mi>
<mi>&alpha;</mi>
</msub>
<msub>
<mi>q</mi>
<mi>&beta;</mi>
</msub>
</mrow>
</mfrac>
<mfrac>
<mrow>
<mi>&Delta;</mi>
<mi>&mu;</mi>
</mrow>
<mi>&rho;</mi>
</mfrac>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>&Delta;</mi>
<mi>&mu;</mi>
</mrow>
<mi>&rho;</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<mn>2</mn>
<msub>
<mi>q</mi>
<mi>&alpha;</mi>
</msub>
<msub>
<mi>q</mi>
<mi>&beta;</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>&Delta;</mi>
<mi>&mu;</mi>
</mrow>
<mi>&rho;</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
<mo>+</mo>
<msubsup>
<mi>q</mi>
<mi>&alpha;</mi>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>q</mi>
<mi>&beta;</mi>
<mn>2</mn>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>&Delta;</mi>
<mi>&mu;</mi>
</mrow>
<mi>&rho;</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mn>4</mn>
</msup>
</mrow>
所述PS波射线参数域的线性近似反射系数为：
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>p</mi>
<mi>s</mi>
</mrow>
</msub>
<mo>&ap;</mo>
<mo>-</mo>
<mfrac>
<msub>
<mi>&alpha;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>&beta;</mi>
<mn>1</mn>
</msub>
</mfrac>
<mfrac>
<mrow>
<mi>&Delta;</mi>
<mi>&mu;</mi>
</mrow>
<mi>&rho;</mi>
</mfrac>
<msub>
<mi>q</mi>
<mi>&alpha;</mi>
</msub>
<mi>p</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<msup>
<mi>p</mi>
<mn>2</mn>
</msup>
<mrow>
<msub>
<mi>q</mi>
<mi>&alpha;</mi>
</msub>
<msub>
<mi>q</mi>
<mi>&beta;</mi>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
p是地下第i层介质的射线参数，θi和分别代表了入射角和反射角；αi,βi和ρi分别代表了纵波速度、横波速度和密度，α、β和ρ代表了上述三个参数在相邻介质的平均值,纵波相关的反射系数Rf可表示为其中和是纵波垂直慢度，qα和qβ分别代表了纵波和横波的平均垂直慢度；Δμ＝ρi+1βi+1 2-ρiβi 2是剪切模量差；
所述广义线性反演方法为：通过用泰勒展开把非线性反问题化成线性反问题进行反演求解；将PP、PS波线性近似反射系数公式对弹性参数模型求一阶泰勒展开式，可得
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>P</mi>
<mi>P</mi>
<mo>/</mo>
<mi>P</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>P</mi>
<mi>P</mi>
<mo>/</mo>
<mi>P</mi>
<mi>S</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mo>&ap;</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>P</mi>
<mi>P</mi>
<mo>/</mo>
<mi>P</mi>
<mi>S</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>m</mi>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>-</mo>
<msub>
<mi>m</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
其中RPP/ps是由步骤二求得的实际数据的射线参数域反射系数值，RPP/ps(m)是由PP、PS波线性近似反射系数公式求得的理论反射系数值，m0＝[α0,β0,ρ0]是初始模型向量，m＝[α,β,ρ]是待解的模型向量；
通过迭代计算，不断求解(5)式并修改模型参数，得到模型的逼近值，即最终的反演结果。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310520644.9A CN104570085B (zh) | 2013-10-29 | 2013-10-29 | 一种纵横波射线参数域联合反演方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310520644.9A CN104570085B (zh) | 2013-10-29 | 2013-10-29 | 一种纵横波射线参数域联合反演方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104570085A CN104570085A (zh) | 2015-04-29 |
| CN104570085B true CN104570085B (zh) | 2017-10-03 |
Family
ID=53086599
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310520644.9A Active CN104570085B (zh) | 2013-10-29 | 2013-10-29 | 一种纵横波射线参数域联合反演方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104570085B (zh) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104965224B (zh) * | 2015-06-03 | 2017-10-27 | 北京多分量地震技术研究院 | 用平均入射角道集进行pp波与ps波联合avo反演方法 |
| CN106842291B (zh) * | 2015-12-04 | 2020-07-14 | 中国石油化工股份有限公司 | 一种基于叠前地震射线阻抗反演的不整合圈闭储层岩性预测方法 |
| CN107092029B (zh) * | 2017-04-26 | 2019-01-18 | 中国石油大学(北京) | 一种地震反演方法和装置 |
| CN108287365B (zh) * | 2018-01-16 | 2019-10-29 | 中国石油大学(华东) | 一种基于波动方程的vsp三参数同步反演方法及装置 |
| CN110879411B (zh) * | 2018-09-05 | 2021-09-17 | 中国石油化工股份有限公司 | 一种基于两参数的多波联合叠前属性提取方法及系统 |
| CN109270574B (zh) * | 2018-10-18 | 2019-09-03 | 国家海洋局第二海洋研究所 | 一种基于多种震源采集数据的联合反褶积方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102129086A (zh) * | 2010-12-10 | 2011-07-20 | 中国石油天然气集团公司 | 用pp、ps波叠前数据联合反演纵、横波阻抗的方法 |
| CN102692645A (zh) * | 2012-06-01 | 2012-09-26 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 利用纵波、转换波数据联合反演储层纵横波速度比的方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6611764B2 (en) * | 2001-06-08 | 2003-08-26 | Pgs Americas, Inc. | Method and system for determining P-wave and S-wave velocities from multi-component seismic data by joint velocity inversion processing |
-
2013
- 2013-10-29 CN CN201310520644.9A patent/CN104570085B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102129086A (zh) * | 2010-12-10 | 2011-07-20 | 中国石油天然气集团公司 | 用pp、ps波叠前数据联合反演纵、横波阻抗的方法 |
| CN102692645A (zh) * | 2012-06-01 | 2012-09-26 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 利用纵波、转换波数据联合反演储层纵横波速度比的方法 |
Non-Patent Citations (3)
| Title |
|---|
| P-SV波反射系数计算方法综述;马昭军 等;《勘探地球物理进展》;20071231;第30卷(第6期);第733-439页 * |
| Zeoppritz-based joint AVO inversion of PP and PS waves in Ray parameters domain;Tiansheng Chen et al.;《SEG Las Vegas 2012 Annual Meeting》;20121231;第1-5页 * |
| 修正柯西约束地震盲反褶积方法;张繁昌 等;《石油地球物理勘探》;20080831;第43卷(第4期);第391-396页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104570085A (zh) | 2015-04-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104570085B (zh) | 一种纵横波射线参数域联合反演方法 | |
| CN113759425B (zh) | 井震联合评价深层古岩溶储层充填特征的方法与系统 | |
| CN104597490B (zh) | 基于精确Zoeppritz方程的多波AVO储层弹性参数反演方法 | |
| CN102937721B (zh) | 利用初至波走时的有限频层析成像方法 | |
| CN102033242B (zh) | 一种深层倾斜裂缝储层地震振幅预测方法 | |
| CN102692645B (zh) | 利用纵波、转换波数据联合反演储层纵横波速度比的方法 | |
| CN101872024B (zh) | 一种利用时移地震进行井位部署的方法 | |
| US6721661B2 (en) | Method of distinguishing types of geologic sedimentation | |
| CN109669212B (zh) | 地震数据处理方法、地层品质因子估算方法与装置 | |
| CN108845351A (zh) | 一种vsp地震资料转换波全波形反演方法 | |
| CN103675911B (zh) | 一种基于纵波和转换横波联合反演截距和梯度的方法 | |
| CN109188511A (zh) | 一种砂泥岩薄互层介质多波avo联合反演方法 | |
| CN103777245A (zh) | 基于地震资料的油气成藏条件定量评价方法 | |
| CN108508489B (zh) | 一种基于波形微变化匹配的地震反演方法 | |
| CN104991272A (zh) | 一种针对无井地震反演的地震速度扰动建模方法 | |
| CN103792578B (zh) | 一种频散avo相关属性反演的流体识别方法 | |
| ITMI951653A1 (it) | Metodo per l'aggiornamento continuo dell'immagine sismica bidimensionale e tridimensionale in profondita' mediante | |
| CN104216013B (zh) | 基于宽方位角资料的c3相干体的方法 | |
| CN103472482B (zh) | 基于基因排序体系的多波地震资料时间域匹配方法及系统 | |
| CN106842291B (zh) | 一种基于叠前地震射线阻抗反演的不整合圈闭储层岩性预测方法 | |
| CN106226814B (zh) | 利用转换横波地震数据反演储层横波阻抗和密度的方法 | |
| Wills et al. | Time-lapse measurements of shallow horizontal wave velocity over a compacting field | |
| CN110703326B (zh) | 一种基于中小偏移距道集的fvo反演方法 | |
| CN104111475A (zh) | 自适应高精度协方差规则化地震数据叠加速度分析方法 | |
| CN104297781A (zh) | 一种射线参数域无约束弹性参数反演方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |