CN108398720A - 一种基于杨氏模量、泊松比的两项式地震叠前反演方法 - Google Patents
一种基于杨氏模量、泊松比的两项式地震叠前反演方法 Download PDFInfo
- Publication number
- CN108398720A CN108398720A CN201810186328.5A CN201810186328A CN108398720A CN 108398720 A CN108398720 A CN 108398720A CN 201810186328 A CN201810186328 A CN 201810186328A CN 108398720 A CN108398720 A CN 108398720A
- Authority
- CN
- China
- Prior art keywords
- modulus
- poisson
- young
- ratio
- formula
- 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 41
- 238000002922 simulated annealing Methods 0.000 claims abstract description 8
- 238000010276 construction Methods 0.000 claims abstract description 5
- 238000000605 extraction Methods 0.000 claims abstract description 5
- 102000008297 Nuclear Matrix-Associated Proteins Human genes 0.000 claims description 7
- 108010035916 Nuclear Matrix-Associated Proteins Proteins 0.000 claims description 7
- 239000011159 matrix material Substances 0.000 claims description 7
- 210000000299 nuclear matrix Anatomy 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000002310 reflectometry Methods 0.000 claims description 5
- 238000012937 correction Methods 0.000 claims description 3
- 230000004927 fusion Effects 0.000 claims description 3
- 238000013213 extrapolation Methods 0.000 claims description 2
- 238000007796 conventional method Methods 0.000 abstract description 2
- 238000005516 engineering process Methods 0.000 description 4
- 239000010410 layer Substances 0.000 description 4
- 239000011435 rock Substances 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 239000003079 shale oil Substances 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 208000035126 Facies Diseases 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- RAFZYSUICBQABU-HMMYKYKNSA-N Phytal Chemical compound CC(C)CCCC(C)CCCC(C)CCC\C(C)=C\C=O RAFZYSUICBQABU-HMMYKYKNSA-N 0.000 description 1
- RAFZYSUICBQABU-QYLFUYDXSA-N Phytal Natural products CC(C)CCC[C@@H](C)CCC[C@@H](C)CCC\C(C)=C/C=O RAFZYSUICBQABU-QYLFUYDXSA-N 0.000 description 1
- 241000219000 Populus Species 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- UZVHFVZFNXBMQJ-UHFFFAOYSA-N butalbital Chemical compound CC(C)CC1(CC=C)C(=O)NC(=O)NC1=O UZVHFVZFNXBMQJ-UHFFFAOYSA-N 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 235000021185 dessert Nutrition 0.000 description 1
- 239000000686 essence Substances 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000002365 multiple layer Substances 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- RAFZYSUICBQABU-UHFFFAOYSA-N phytenal Natural products CC(C)CCCC(C)CCCC(C)CCCC(C)=CC=O RAFZYSUICBQABU-UHFFFAOYSA-N 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 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
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)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明属于地震勘探技术领域，公开了一种基于杨氏模量、泊松比的两项式地震叠前反演方法，所述基于杨氏模量、泊松比的两项式地震叠前反演方法包括以下步骤：叠前地震资料预处理；构造层位解释；测井资料外推与模型构建；多角度子波提取；最终利用模拟退火反演方法获得杨氏模量和泊松比数据体。本发明通过发展两项AVO反射系数方程，将密度项以合理的方式与纵波速度建立数学关系，消除密度项，可以在缺少大角度叠前地震资料的情况下进行叠前直接反演，得到目标弹性参数，可以减小传统方法带来的累积误差。
Description
技术领域
本发明属于地震勘探技术领域，尤其涉及一种基于杨氏模量、泊松比的两项式地震叠前反演方法。
背景技术
当前全球地震勘探技术发展的几个主要趋势包括：从常规油气向常规与非常规油气并重局面发展、从简单地表和浅水区向复杂地表和深水区发展、从叠后地震反演向叠前弹性反演发展。页岩油气的开采需压裂页岩生储层，因此岩区的易压裂程度成为评价优质页岩油气储层非常重要的指标。岩石物理研究表明，杨氏模量和泊松比能够较好的表征岩体的造缝难易程度，标识页岩气的“甜点”区；由于地震测线长度与储层埋深的比值通常很小，故通常很难获得超深储层的大反射倾角叠前地震数据。叠前AVO反演可以充分利用叠前地震资料包含的地质、岩性和流体信息，获取丰富的地下介质弹性参数。AVO：(Amplitudevariation with offset)利用地震反射波振幅随偏移距的变化特征来探讨反射系数响应随入射角的变化,进而确定反射界面上覆、下伏介质的岩性特征及物性参数的技术。主要用于石油天然气勘探。叠前AVO反演的核心是反射系数方程。通常反射系数方程由三项构成，其中之一为密度项。密度项通常不作为直接用来预测油气储层的指标，原因是具体工区的密度变化范围不大，而各种岩性的理论密度范围交集较广，且其对油气储层特征的直接相关性差，单由密度项无法进行岩性、流体的判断。在缺乏大角度地震资料或者地震资料数据信噪比较低的情况下，大角度反演的结果品质较差，提取参数的过程中，参数矩阵的稳定性较差，造成最终提取的弹性参数误差较大。另外，由于密度项的参数权值在角度有限情况下较小，反演的难度更加大。以YPD三项反射系数方程为基础，构建的YPD叠前地震反演方法，同时反演了杨氏模量、泊松比与密度。该方法为目前较新的杨氏模量、泊松比反演方法(2014年建立)。具体做法如下：假设待反演杨氏模量、泊松比及密度反射系数服从柯西分布，假设似然函数服从高斯分布，同时在反演目标函数中加入初始模型约束，并通过初始模型建立各道去相关矩阵，消除待反演参数间的互相关性，在贝叶斯反演框架下进行反演。该方法在叠前地震数据远角道集缺失的情况下，大角度入射角的引入会使参数矩阵的稳定性降低，造成最终提取的弹性参数存在较大误差。
综上所述，现有技术存在的问题是：针对深层、超深层页岩气储层的地震勘探，由于远角道集数据的缺失，使得利用AVO叠前反演较难获取地下介质准确的杨氏模量、泊松比等岩石物理弹性参数。因此，提高参数矩阵稳定性，降低计算带来的累积误差成为亟需解决的问题。
发明内容
针对现有技术存在的问题，本发明提供了一种基于杨氏模量、泊松比的两项式地震叠前反演方法。
本发明是这样实现的，一种基于杨氏模量、泊松比的两项式地震叠前反演方法，所述基于杨氏模量、泊松比的两项式地震叠前反演方法包括以下步骤：
步骤一，叠前地震资料预处理；
步骤二，构造层位解释；
步骤三，测井资料外推与模型构建；
步骤四，多角度子波提取；
步骤五，最终利用模拟退火反演方法获得杨氏模量和泊松比数据体。
进一步，所述基于杨氏模量、泊松比的两项式地震叠前反演方法具体包括以下步骤：
(1)基于Aki-Richards近似方程推导得到YP两项式反射系数近似方程，作为深层杨氏模量、泊松比叠前反演方法的理论基础；
Aki-Richards近似方程表述为：
其中，Vp、Vs分别为纵、横波速度，ρ为密度，θ为入射角。分别是纵、横波速度反射系数与密度反射系数。
(2)采用幂指数拟合的方式，引入参数可变的关系式消除密度项，由密度与纵波速度的微商关系：
dρ＝FLVp L-1dVp (2)
可得
纵横波模量反射系数与纵横波速度、密度反射系数的关系如下所示：
将(3)带入(4)(5)得
令
(3)整理杨氏模量、泊松比关于纵横波速度的关系，并融合以上各式，得
将(9)、(10)、(11)带入Aki-Richard方程(1)，得到新的两项式近似方程
其中，k为横、纵波速度比的平方，L为幂指数。
进一步，所述基于杨氏模量、泊松比的两项式地震叠前反演方法在进行杨氏模量与泊松比反演计算以前，需从先验信息中求得参数L、k；在实施操作过程中，其中k为测井资料中横、纵波速度比的平方，L则通过密度与纵波速度的幂指数拟合计算得到。
待反演参数的计算可简化为求解：
其中，CE和Cσ与角度有关；RE和Rσ分别是杨氏模量、泊松比的反射系数时间序列
反射系数与某一角度子波卷积形成一个固定角度的地震合成记录。将式(13)等号两边同时与子波卷积可推导为式
SPP(t，θi)是入射角为θi的角道集，W(t，θi)是一个子波核矩阵。
卷积运算之后，可得
若设D＝Gm，SPP表示包含向量D的叠前角度道集，WE和Wσ构成小波核矩阵G，RE和Rσ是包含m的反射率项。
反演的目标函数为
f(V)＝||S-D||→min (21)
式中，V为待反演弹性参数构成的参数向量；S(t)＝W*R为初始低频模型在各角度下的合成地震记录，W为一组不同角度地震子波，R为计算得到的低频地震模型的反射系数；D为实际部分叠加道集地震数据；
构造f(V)的Jacobian矩阵J，通过模拟退火算法对J进行求解，得到对V的修正量ΔV0，将V0＝(V+ΔV0)作为新的反演参数向量带入式(21)，反复迭代求解f(V)的极小值，直至求得使f(V)达到允许误差的ΔVn，此时认为Vn＝(Vn-1+ΔVn)为最优解；对m的多次迭代最终值即为反演得到的最终杨氏模量、泊松比数据体。
本发明的优点及积极效果为：通过发展两项AVO反射系数方程，将密度项以合理的方式与纵波速度建立数学关系，消除密度项，可以在缺少大角度叠前地震资料的情况下进行叠前直接反演，得到目标弹性参数，这样可以减小传统方法带来的累积误差；实际数据来自CGG Hampson-Russell公司开发软件的AVO模块的测试数据。对地震资料进行预处理后抽取角道集数据，利用HRS软件依据测井资料对角道集数据提取五组不同角度的子波，并依据所建立的不同子波核矩阵与YP近似方程计算的反射系数合成出不同角度域地震剖面，带入反演流程进行迭代计算；得到最终反演结果的杨氏模量、泊松比反演剖面，反演结果与实际情况较为吻合，显示了本发明较好的应用潜力。
附图说明
图1是本发明实施例提供的基于杨氏模量、泊松比的两项式地震叠前反演方法流程图。
图2是本发明实施例提供的YP近似方程地震叠前反演流程图。
图3是本发明实施例提供的正波阻抗界面反射系数对比示意图；
图中：(a)不同方程反射系数对比；(b)不同方程反射系数差值对比。
图4是本发明实施例提供的负波阻抗界面反射系数对比示意图；
图中：(a)不同方程反射系数对比；(b)不同方程反射系数差值对比。
图5是本发明实施例提供的层状介质模型反演结果示意图；
图中：(a)角度域合成地震记录；(b)杨氏模量、泊松比反演曲线对比。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
下面结合附图对本发明的应用原理作详细的描述。
如图1所示，本发明实施例提供的基于杨氏模量、泊松比的两项式地震叠前反演方法包括以下步骤：
S101：叠前地震资料预处理；
S102：地质层位解释；
S103：测井资料外推与模型构建；
S104：多角度子波提取；
S105：最终利用模拟退火反演方法获得杨氏模量和泊松比数据体。
本发明实施例提供的基于杨氏模量、泊松比的两项式地震叠前反演方法具体包括以下步骤：
首先，基于Aki-Richards近似方程推导得到YP两项式反射系数近似方程，作为深层杨氏模量、泊松比叠前反演方法的理论基础。
Aki-Richards近似方程表述为：
其中，Vp、Vs分别为纵、横波速度，ρ为密度，θ为入射角。 分别是纵、横波速度反射系数与密度反射系数。
由于可认为密度项是一个与泊松比等信息相关的非独立参数，且相对于其他参数，其变化范围不大，因此可以与纵波速度建立关系。Gardner针对砂泥岩建立了密度与纵波速度的数学联系，而后Gardner公式成为了将岩体密度拟合纵波速度的经典公式。采用幂指数拟合的方式，引入参数可变的关系式消除密度项，以此保证对实际工区拟合的准确性。
由密度与纵波速度的微商关系
dρ＝FLVp L-1dVp (2)
可得
纵横波模量反射系数与纵横波速度、密度反射系数的关系如下所示：
将(3)带入(4)(5)得
令
整理杨氏模量、泊松比关于纵横波速度的关系，并融合以上各式，得
将(9)、(10)、(11)带入Aki-Richard方程(1)，得到新的两项式近似方程
其中，k为横、纵波速度比的平方，L为幂指数。可将上式称为YP近似方程，该式建立了杨氏模量反射率、泊松比反射率与反射系数的直接线性关系。
本发明开展基于YP近似方程的两项AVO反演，其主要流程包含叠前地震资料预处理、地质层位解释、测井资料外推与模型构建、多角度子波提取等环节，最终利用模拟退火反演方法获得杨氏模量和泊松比数据体。主要反演流程如图所示：在进行杨氏模量与泊松比反演计算以前，需从先验信息中求得参数L、k。在实施操作过程中，为确保参数拟合结果准确可靠，针对目的层段可采用测井资料的纵、横波速度与密度确定L、k，其中k为测井资料中横、纵波速度比的平方，L则通过密度与纵波速度的幂指数拟合计算得到。
待反演参数的计算可简化为求解
其中，CE和Cσ与角度有关；RE和Rσ分别是杨氏模量、泊松比的反射系数时间序列
反射系数与某一角度子波卷积形成一个固定角度的地震合成记录。将式(13)等号两边同时与子波卷积可推导为式
SPP(t，θi)是入射角为θi的角道集，W(t，θi)是一个子波核矩阵。
卷积运算之后，可得
若设D＝Gm，SPP表示包含向量D的叠前角度道集，WE和Wσ构成小波核矩阵G，RE和Rσ是包含m的反射率项。
反演的目标函数为
f(V)＝||S-D||→min (21)
式中，V为待反演弹性参数构成的参数向量；S(t)＝W*R为初始低频模型在各角度下的合成地震记录，W为一组不同角度地震子波，R为计算得到的低频地震模型的反射系数；D为实际部分叠加道集地震数据；
构造f(V)的Jacobian矩阵J，通过模拟退火算法对J进行求解，得到对V的修正量ΔV0，将V0＝(V+ΔV0)作为新的反演参数向量带入式(21)，反复迭代求解f(V)的极小值，直至求得使f(V)达到允许误差的ΔVn，此时认为Vn＝(Vn-1+ΔVn)为最优解；对m的多次迭代最终值即为反演得到的最终杨氏模量、泊松比数据体。下面结合具体应用实施例对本发明的应用效果作详细的描述。
位于四川的元坝气田是迄今为止中国发现的埋藏最深的海相气田，选择元坝气田实测井资料构造地质模型对精确Zoeppritz方程、Aki-Richards近似方程与YP近似方程的计算精度进行实验对比分析，模型参数如表1所示。
表1含水饱和度(2014)地震模型
波速与密度的乘积被称为波阻抗，正波阻抗界面的入射介质波阻抗小于投射介质波阻抗，负波阻抗界面则相反。在本模型中，上覆地层与气层的分界面为负波阻抗界面，气层与下伏地层的分界面为正波阻抗界面。分别用精确的Zoeppritz方程、Aki-Richard近似方程、YP近似方程计算上述模型不同界面处的反射系数及近似方程与精确方程的残差。图3和图4分别为正、负波阻抗界面对比分析，图3a、图4a为分别用精确的Zoeppritz方程、Aki-Richard近似方程、YP近似方程计算得到的反射系数随入射角的变化，图3b、图4b为Aki-Richard近似方程、YP近似方程计算得到的反射系数与精确Zoeppritz方程计算得到的反射系数的差值随入射角的变化。
由图3a、图4a可见，基于YP近似方程计算的反射系数与精确Zoeppritz方程有一定误差，但精度值稳定优于Aki-Richard近似方程的计算结果。由图3b、图4b可见，Aki-Richard近似方程对Zoeppritz方程的残差在入射角为50°时开始急剧增大。相对而言，YP近似方程的误差在入射角为60°以内时相当稳定而微小，证明了YP近似方程具有的实际可行性。
为验证YP两项AVO反演的可行性与稳定性，本发明采用多种模型进行实验。首先设置一组具有一定地层厚度、共16层的多层水平层状介质模型，赋予其纵、横波速度与密度。子波采用主频为40Hz的雷克子波，采用精确Zoeppritz方程进行正演得到角度域叠前角道集，采用本发明所述YP近似方程进行杨氏模量、泊松比的直接反演，结果如图5所示。
图5(a)为正演所得角度域叠前角道集，图5(b)蓝色曲线为地层参数计算所得的杨氏模量、泊松比曲线，红色曲线为反演所得杨氏模量、泊松比曲线。由反演结果可见，利用YP近似反演能够获取与真实值基本吻合的杨氏模量和泊松比。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
Claims (3)
1.一种基于杨氏模量、泊松比的两项式地震叠前反演方法，其特征在于，所述基于杨氏模量、泊松比的两项式地震叠前反演方法包括以下步骤：
步骤一，叠前地震资料预处理；
步骤二，构造层位解释；
步骤三，测井资料外推与模型构建；
步骤四，多角度子波提取；
步骤五，最终利用模拟退火反演方法获得杨氏模量和泊松比数据体。
2.如权利要求1所述的基于杨氏模量、泊松比的两项式地震叠前反演方法，其特征在于，所述基于杨氏模量、泊松比的两项式地震叠前反演方法具体包括以下步骤：
(1)基于Aki-Richards近似方程推导得到YP两项式反射系数近似方程，作为深层杨氏模量、泊松比叠前反演方法的理论基础；
Aki-Richards近似方程表述为：
其中，Vp、Vs分别为纵、横波速度，ρ为密度，θ为入射角， 分别是纵、横波速度反射系数与密度反射系数；
(2)采用幂指数拟合的方式，引入参数可变的关系式消除密度项，由密度与纵波速度的微商关系：
dρ＝FLVp L-1dVp (2)
可得
纵横波模量反射系数与纵横波速度、密度反射系数的关系如下所示：
将(3)带入(4)(5)得
令
(3)整理杨氏模量、泊松比关于纵横波速度的关系，并融合以上各式，得
将(9)、(10)、(11)带入Aki-Richard方程(1)，得到新的两项式近似方程
其中，k为横、纵波速度比的平方，L为幂指数。
3.如权利要求1所述的基于杨氏模量、泊松比的两项式地震叠前反演方法，其特征在于，所述基于杨氏模量、泊松比的两项式地震叠前反演方法在进行杨氏模量与泊松比反演计算以前，需从先验信息中求得参数L、k；在实施操作过程中，其中k为测井资料中横、纵波速度比的平方，L则通过密度与纵波速度的幂指数拟合计算得到；
待反演参数的计算可简化为求解：
其中，CE和Cσ与角度有关；RE和Rσ分别是杨氏模量、泊松比的反射系数时间序列
反射系数与某一角度子波卷积形成一个固定角度的地震合成记录，将式(13)等号两边同时与子波卷积可推导为式
SPP(t，θi)＝∫0 tRPP(τ，θi)*wavelet(t-τ)dτ (18)
W(t，θi)＝∫0 tC(τ，θi)*wavelet(t-τ)dτ (19)
SPP(t，θi)是入射角为θi的角道集，W(t，θi)是一个子波核矩阵；
卷积运算之后，可得
若设D＝Gm，SPP表示包含向量D的叠前角度道集，WE和Wσ构成小波核矩阵G，RE和Rσ是包含m的反射率项；
反演的目标函数为
f(V)＝||S-D||→min (21)
式中，V为待反演弹性参数构成的参数向量；S(t)＝W*R为初始低频模型在各角度下的合成地震记录，W为一组不同角度地震子波，R为计算得到的低频地震模型的反射系数；D为实际部分叠加道集地震数据；
构造f(V)的Jacobian矩阵J，通过模拟退火算法对J进行求解，得到对V的修正量ΔV0，将V0＝(V+ΔV0)作为新的反演参数向量带入式(21)，反复迭代求解f(V)的极小值，直至求得使f(V）达到允许误差的ΔVn，此时认为Vn＝(Vn-1+ΔVn)为最优解；对m的多次迭代最终值即为反演得到的最终杨氏模量、泊松比数据体。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810186328.5A CN108398720A (zh) | 2018-03-07 | 2018-03-07 | 一种基于杨氏模量、泊松比的两项式地震叠前反演方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810186328.5A CN108398720A (zh) | 2018-03-07 | 2018-03-07 | 一种基于杨氏模量、泊松比的两项式地震叠前反演方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108398720A true CN108398720A (zh) | 2018-08-14 |
Family
ID=63092086
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810186328.5A Pending CN108398720A (zh) | 2018-03-07 | 2018-03-07 | 一种基于杨氏模量、泊松比的两项式地震叠前反演方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108398720A (zh) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109471165A (zh) * | 2018-12-03 | 2019-03-15 | 中国石油化工股份有限公司 | 基于包含敏感岩性识别因子为变量的avo近似式叠前反演方法 |
| CN109490964A (zh) * | 2018-11-12 | 2019-03-19 | 同济大学 | 一种改进的高精度avo弹性参数快速反演方法 |
| CN110941009A (zh) * | 2018-09-21 | 2020-03-31 | 中国石油化工股份有限公司 | 基于时频联合域反演的杨氏模量和泊松比预测方法及装置 |
| CN111239833A (zh) * | 2020-03-06 | 2020-06-05 | 中海石油(中国)有限公司 | 一种基于泊松比分解的k值鲁棒ypd叠前同时反演方法 |
| CN112904430A (zh) * | 2019-12-03 | 2021-06-04 | 中国石油化工股份有限公司 | 非线性直接叠前地震泊松阻抗反演的计算机实现的方法 |
| CN113376692A (zh) * | 2020-03-10 | 2021-09-10 | 中国石油天然气集团有限公司 | 致密砂岩气水平井压裂改造方案优化方法及装置 |
-
2018
- 2018-03-07 CN CN201810186328.5A patent/CN108398720A/zh active Pending
Non-Patent Citations (2)
| Title |
|---|
| 张瑞: "深层叠前反演与储层评价", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
| 张进: "地震叠前数据的弹性阻抗非线性反演方法研究", 《中国优秀博士学位论文全文数据库 基础科学辑》 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110941009A (zh) * | 2018-09-21 | 2020-03-31 | 中国石油化工股份有限公司 | 基于时频联合域反演的杨氏模量和泊松比预测方法及装置 |
| CN109490964A (zh) * | 2018-11-12 | 2019-03-19 | 同济大学 | 一种改进的高精度avo弹性参数快速反演方法 |
| CN109471165A (zh) * | 2018-12-03 | 2019-03-15 | 中国石油化工股份有限公司 | 基于包含敏感岩性识别因子为变量的avo近似式叠前反演方法 |
| CN112904430A (zh) * | 2019-12-03 | 2021-06-04 | 中国石油化工股份有限公司 | 非线性直接叠前地震泊松阻抗反演的计算机实现的方法 |
| US11493658B2 (en) | 2019-12-03 | 2022-11-08 | China Petroleum & Chemical Corporation | Computer-implemented method and system employing nonlinear direct prestack seismic inversion for poisson impedance |
| CN111239833A (zh) * | 2020-03-06 | 2020-06-05 | 中海石油(中国)有限公司 | 一种基于泊松比分解的k值鲁棒ypd叠前同时反演方法 |
| CN113376692A (zh) * | 2020-03-10 | 2021-09-10 | 中国石油天然气集团有限公司 | 致密砂岩气水平井压裂改造方案优化方法及装置 |
| CN113376692B (zh) * | 2020-03-10 | 2023-05-26 | 中国石油天然气集团有限公司 | 致密砂岩气水平井压裂改造方案优化方法及装置 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108398720A (zh) | 一种基于杨氏模量、泊松比的两项式地震叠前反演方法 | |
| Regone et al. | Geologic model building in SEAM Phase II—Land seismic challenges | |
| CN101329405B (zh) | 一种简单的多参数地震反演方法 | |
| US7493241B2 (en) | 3D velocity modeling, with calibration and trend fitting using geostatistical techniques, particularly advantageous for curved for curved-ray prestack time migration and for such migration followed by prestack depth migration | |
| US6269310B1 (en) | System for eliminating headwaves in a tomographic process | |
| US10234582B2 (en) | Joint inversion of seismic data | |
| CN108363100A (zh) | 基于封存单元与粗糙集的煤层气甜点区地震地质识别方法 | |
| CN104516018A (zh) | 一种地球物理勘探中岩性约束下的孔隙度反演方法 | |
| CN106405651A (zh) | 一种基于测井匹配的全波形反演初始模型构建方法 | |
| Huang et al. | Geological structure-guided initial model building for prestack AVO/AVA inversion | |
| US20160320504A1 (en) | System and method for processing microseismic data | |
| CN113740901B (zh) | 基于复杂起伏地表的陆上地震数据全波形反演方法及装置 | |
| CN111722284B (zh) | 一种基于道集数据建立速度深度模型的方法 | |
| CN105089652A (zh) | 一种拟声波曲线重构与稀疏脉冲联合反演方法 | |
| CN107817526A (zh) | 叠前地震道集分段式振幅能量补偿方法及系统 | |
| CN113031068B (zh) | 一种基于反射系数精确式的基追踪叠前地震反演方法 | |
| CN109696704B (zh) | 一种基于纵波阻抗约束的地震各向异性δ建模方法 | |
| CN109143352B (zh) | 一种各向异性介质地震反射特征方程建立方法 | |
| CN107238859A (zh) | 计算虚拟井合成地震记录的方法及系统 | |
| Velásquez et al. | Depth-conversion techniques and challenges in complex sub-Andean provinces | |
| CN104316959B (zh) | 一种基于等效流体声波阻抗的流体识别方法 | |
| Huang et al. | Reservoir characterization by integration of seismic and dynamic data | |
| CN112925022B (zh) | 页岩vti介质各向异性参数的预测方法 | |
| US20230266491A1 (en) | Method and system for predicting hydrocarbon reservoir information from raw seismic data | |
| CN116338785A (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: 20180814 |
|
| RJ01 | Rejection of invention patent application after publication |