(19)中华人民共和国国家知识产权局
(12)发明专利说明书 | ||
(10)申请公布号 CN 103257361 A (43)申请公布日 2013.08.21 | ||
(21)申请号 CN201310198497.8
(22)申请日 2013.05.24
(71)申请人 中国石油天然气集团公司;中国石油大学(北京)
地址 100007 北京市东城区东直门北大街9号
(72)发明人 袁三一 柴新涛 施佩东 王尚旭
正则化反演(74)专利代理机构 北京三友知识产权代理有限公司
代理人 任默闻
(51)Int.CI
G01V1/28
权利要求说明书 说明书 幅图 |
(54)发明名称
基于Zoeppritz方程近似式的油气预测方法及系统 | |
(57)摘要
本发明提供一种基于Zoeppritz方程近似式的油气预测方法及系统，方法包括：采集当前目标区的测井数据以及叠前角道集数据；根据测井数据以及叠前角道集数据提取地震子波；根据叠前角道集数据、地震子波以及Zoeppritz方程近似式得到反演目标函数；根据第一关系、第二关系以及Zoeppritz方程近似式确定目标Zoeppritz方程近似式；根据反演目标函数以及目标Zoeppritz方程近似式确定纵波模量以及横波模量；根据当前目标区地层的纵波模量以及横波模量确定流体因子；对流体因子进行综合解释，得到当前目标区地层的油气藏预测结果。可稳定求取纵波、横波模量信息，提升了流体因子的精度，进而可降低油气藏预测的风险。 | |
法律状态
法律状态公告日 | 法律状态信息 | 法律状态 |
权 利 要 求 说 明 书
1.一种基于佐普利兹Zoeppritz方程近似式的油气预测方法，其特征是， 所述的方法包括：
采集当前目标区的测井数据以及叠前角道集数据；
根据所述的测井数据以及叠前角道集数据提取地震子波；
根据所述的叠前角道集数据、所述的地震子波以及Zoeppritz方程近似式 得到反演目标函数；
根据纵波速度、纵波模量以及密度的弹性参数关系推导纵波模量反射系 数、纵波速度
反射系数与密度反射系数的关系，称为第一关系；
根据横波速度、横波模量以及密度的弹性参数关系推导横波模量反射系 数、横波速度反射系数与密度反射系数的关系，称为第二关系；
根据所述的第一关系、第二关系以及Zoeppritz方程近似式确定目标 Zoeppritz方程近似式；
根据所述的反演目标函数以及目标Zoeppritz方程近似式确定当前目标区 地层的纵波模量以及横波模量；
根据当前目标区地层的纵波模量以及横波模量确定流体因子；
对所述的流体因子进行综合解释，得到当前目标区地层的油气藏预测结 果。
2.根据权利要求1所述的方法，其特征是，所述的Zoeppritz方程近似式 为
<maths><math><mrow><mfrac><mn>1</mn><mn>2</mn></mfrac><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msup><mi>tan</mi><mn>2</mn></msup><mover><mi>θ</mi
><mo>‾</mo></mover><mo>)</mo></mrow><mfrac><mrow><mi>Δ</mi><msub><mi>V</mi><mi>p</mi></msub></mrow><mover><msub><mi>V</mi><mi>p</mi></msub><mo>‾</mo></mover></mfrac><mo>-</mo><mn>4</mn><msup><mi>sin</mi><mn>2</mn></msup><mover><mi>θ</mi><mo>‾</mo></mover><mfrac><msubsup><mover><mi>V</mi><mo>‾</mo></mover><mi>s</mi><mn>2</mn></msubsup><msubsup><mover><mi>V</mi><mo>‾</mo></mover><mi>p</mi><mn>2</mn></msubsup></mfrac><mfrac><mrow><mi>Δ</mi><msub><mi>V</mi><mi>s</mi></msub></mrow><mover><msub><mi>V</mi><mi>s</mi></msub><mo>‾</mo></mover></mfrac><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mn>4</mn><msup><mi>sin</mi><mn>2</mn></msup><mover><mi>θ</mi><mo>‾</mo></mover><mfrac><msubsup><mover><mi>V</mi><mo>‾</mo></mover><mi>s</mi><mn>2</mn></msubsup><msubsup><mover><mi>V</mi><mo>‾</mo></mover><mi>p</mi><mn>2</mn></msubsup></mfrac><mo>)</mo></mrow><mfrac><mi>Δρ</mi><mover><mi>ρ</mi><mo>‾</mo></mover></mfrac><mo>=</mo><msub><mi>r</mi><mi>pp</mi></msub><mrow><mo>(</mo><mover><mi>θ</mi><mo>‾</mo></mover><mo>)</mo></mrow></mrow></math></maths>
其中，为地层分界面处上下两弹性介质的平均纵波速度； 为地层分界面处上下两弹性介质的平均横波速度；为地层分界面处上下两弹性介质的平均密度；ΔV<sub>p</sub>=V<sub>p2</sub>-V<sub>p1</sub>为上下两弹性介质 纵波速度的差值；ΔV<sub>s</sub>=V<sub>s2</sub>-V<sub>s1</sub>为上下两弹性介质横波速度的差值；Δρ=ρ<sub>2</sub>-ρ<sub>1</sub>为上下两弹性介质密度的差值；为入射角θ<sub>i</sub>和透射角θ<sub>t</sub>的平均值； 为随角度变化的反射纵波反射系数。
3.根据权利要求1所述的方法，其特征是，根据所述的叠前角道集数据、 所述的地震子波以及Zoeppritz方程近似式得到反演目标函数具体包括：
对所述的叠前角道集数据以及所述的地震子波进行傅里叶变换；
对傅里叶变换后的叠前角道集数据以及地震子波进行频谱分析，并选定 优势频带范围；
根据所述的优势频带范围及Zoeppritz方程近似式确定频率域反演模型；
对所述的频率域反演模型加入1范数正则化，得到反演目标函数。
4.根据权利要求3所述的方法，其特征是，根据所述的优势频带范围以 及Zoeppritz方程近似式确定的频率域反演模型为：
WCm=d
其中，W为所述的地震子波对应的频谱矩阵，C为Zoeppritz方程近似式 中所有入射角θ<sub>i</sub>对应的三参数系数矩阵，m为模型参数，d为所有入射角θ<sub>i</sub>角 道集地震记录的复数谱信息。
5.根据权利要求4所述的方法，其特征是，对所述的频率域反演模型加 入1范数正则化，得到的反演目标函数为：
min‖d-WCm‖<sub>2</sub>+λ‖m‖<sub>1</sub>
其中，λ为一调节因子。
6.根据权利要求1所述的方法，其特征是，所述的第一关系为：
<maths><math><mrow><mfrac><mrow><mi>Δ</mi><msub><mi>V</mi><mi>p</mi></msub></mrow><mover><msub><mi>V</mi><mi>p</mi></msub><mo>‾</mo></mover></mfrac><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mrow><mo>(</mo><mfrac><mi>ΔM</mi><mover><mi>M</mi><mo>‾</mo></mover></mfrac><mo>-</mo><mfrac><mi>Δρ</mi><mover><mi>ρ</mi><mo>‾</mo></mover></mfrac><mo>)</mo></mrow></mrow></math></maths>
版权声明：本站内容均来自互联网，仅供演示用，请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198，我们将在24小时内删除。
发表评论