CN106291677A - 一种基于匹配追踪方法的叠后声波阻抗反演方法 - Google Patents
一种基于匹配追踪方法的叠后声波阻抗反演方法 Download PDFInfo
- Publication number
- CN106291677A CN106291677A CN201510268461.1A CN201510268461A CN106291677A CN 106291677 A CN106291677 A CN 106291677A CN 201510268461 A CN201510268461 A CN 201510268461A CN 106291677 A CN106291677 A CN 106291677A
- Authority
- CN
- China
- Prior art keywords
- impedance
- reflection coefficient
- wavelet
- match tracing
- data
- 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.)
- Granted
Links
Abstract
本发明提供了一种基于匹配追踪方法的叠后声波阻抗反演方法,属于油气及煤层气地震勘探与开发领域。本方法包括:S1,输入叠后地震数据,由叠后地震数据进行构造解释;S2,输入测井数据,提取或给定地震子波,并对叠后地震数据进行标定;同时计算出过井波阻抗:S3,以构造解释结果作为约束,对S2得到的过井波阻抗进行内插外推,获得初始波阻抗体,进而获得初始反射系数序列;S4,对于每一道地震数据,由初始反射系数及子波构建目标函数;S5,用匹配追踪方法求解目标函数,获得反演后的反射系数;S6,由反演后的反射系数以及标定结果计算最终的波阻抗。
Description
技术领域
本发明属于油气及煤层气地震勘探与开发领域,具体涉及一种基于匹配追踪方法的叠后声波阻抗反演方法。
背景技术
随着油气勘探开发的需要,储层预测和精细描述越来越引起重视。围绕这个目的而开展的探索和研究也越来越多,地震反演就是最重要的一点。当反演的目标是波阻抗时,我们称之为波阻抗反演。目前的地震波阻抗反演按照所使用的资料分为叠后反演和叠前反演两大类,按照反演方法又可以分为直接反演和间接反演。直接反演就是直接从地震资料出发,进行运算得到波阻抗。间接反演是从一个初始模型出发,合成地震记录,然后构造目标函数,最后求得目标函数在一定的范数意义下极值。无论哪种反演,其实质都是想去除子波的影响,从而将地震剖面转化为能与钻井、地质等资料直接对比的形式,因此反演在许多情况下提高了常规地震的分辨率并提高了油藏参数研究的水平。
叠后声波阻抗反演概括起来不外乎有两大类:基于反射系数逆公式的直接反演和基于正演模型的迭代反演。
基于反射系数逆公式的直接反演有道积分包括递推反演。地震道积分近似等于对数波阻抗,该方法无法求得地层绝对波阻抗,而且使用时无法用地质或测井资料进行约束。以稀疏脉冲反演为代表的递推反演虽然更加精确和稳定,但是在实现和应用上都更复杂。
基于模型的反演需要由测井、地质以及地震资料给出一个初始模型,然后迭代反演,得到与地震资料最佳匹配的地震波阻抗模型。在现实的实际应用中,这种基于模型的反演其结果往往对初始模型具有一定的依赖性,并且反演结果 具有非唯一性。为了减少反演结果的非唯一性,产生了井约束的波阻抗反演。井约束的反演可以在一定程度上减少反演结果的非唯一性,但是毕竟井给出的是一个点上的信息,这种约束作用在空间上具有一定的局限性。
这两类反演,最终都表示为一个目标函数的最优化问题。通常都采用最小二乘或共轭梯度法来解此最优化问题。
发明内容
本发明的目的在于解决上述现有技术中存在的难题,提供一种基于匹配追踪方法的叠后声波阻抗反演方法,用匹配追踪法来解反射系数估计中的最优化问题,进而获得声波阻抗。
本发明是通过以下技术方案实现的:
一种基于匹配追踪方法的叠后声波阻抗反演方法,包括:
S1,输入叠后地震数据,由叠后地震数据进行构造解释;
S2,输入测井数据,提取或给定地震子波,并对叠后地震数据进行标定;同时计算出过井波阻抗:
S3,以构造解释结果作为约束,对S2得到的过井波阻抗进行内插外推,获得初始波阻抗体,进而获得初始反射系数序列;
S4,对于每一道地震数据,由初始反射系数序列及子波构建目标函数;
S5,用匹配追踪方法求解目标函数,获得反演后的反射系数;
S6,由反演后的反射系数以及标定结果计算最终的波阻抗。
所述S2是这样实现的:
输入测井数据中的速度和密度,利用合成记录对叠后地震数据进行层位标定和子波提取;
合成记录F(t)是地震子波S(t)与反射系数R(t)褶积的结果,即F(t)=S(t)*R(t),初始合成地震记录是由反射系数R(t)与标准地震子波S(t)进行褶积得到;
在初始标定的基础上通过井旁地震道和速度密度测井曲线联合提取子波,子波提取和合成记录制作是一个迭代的过程,经过多次迭代,即得到合适的子波和高精度的合成记录;然后对F(t)=S(t)*R(t)进行反褶积,得到反射系数R(t),再利用反射系数与波阻抗P的关系,即P(n)=[(1+R(n))/(1-R(n))]xP(n-1),计算得到过井波阻抗。
所述S3是这样实现的:
利用构造解释得到的某一层位的过井波阻抗值(通过S2得到的),用插值方法获得该层位上其他点的波阻抗值,然后用获得初始反射系数序列。
所述S4中的目标函数如下:
在约束||dobs-W·R||p≤γ下求解min||R||q其中dobs是观测数据,W为子波,R为反射系数模型,γ为任意小的一个数。这不同于现有技术的通常思路,即在约束||R||q下求min||dobs-W·R||p。目标函数可以在解此问题的过程中根据需要形成,也可以不用形成目标函数。
所述S5是这样实现的:
在约束||dobs-W·R||p≤γ下求解min||R||q,取p=2,q=1,变成二次规划问题,然后采用拉格朗日乘子等方法来解此问题。
所述S6是这样实现的:
利用S5得到的反射系数,然后用获得声波阻抗。
与现有技术相比,本发明的有益效果是:
目前的声波阻抗反演方法都最终表示为一个目标函数最优化的问题,其出发点都是基于数据拟合的思想,即观测数据与合成数据的误差最小,数学上表示为min||dobs-W·R||p,为了得到合适的解再加入先验信息的约束项:λ||R||q,形成 目标函数min||dobs-W·R||p+λ||R||q,然后诸如Newton法或模拟退火算法解此最优化问题,获得反射系数序列,进而获得声波阻抗。本发明借鉴匹配追踪方法的思想而不是数据拟合的思想,即在约束||dobs-W·R||p≤γ下求解min||R||q,这样就将问题转化为二次规划问题,然后采用拉格朗日乘子等方法来解此问题,从而拓宽了问题的解决方法。
附图说明
图1本发明方法的步骤框图。
具体实施方式
下面结合附图对本发明作进一步详细描述:
目前叠后波阻抗反演都基于如下褶积模型:即地震信号可以表示为d=W*R+n。其中d表示地震道数据,W表示地震子波,R表示反射系数序列,*表示褶积,n表示噪声。反演的目的就是由观测到的地震数据d获得反射系数序列R,进而获得波阻抗Z。由反射系数获得波阻抗的过程比较简单,目前的反演方法都聚焦在如何获得反射系数。而且目前通常都是将反演问题表示为如下形式的最优化问题:min||dobs-W·R||p+λ||R||q。对于这个最优化问题,通常的做法是令p=2,q=2,用最小二乘法或者共轭梯度法解此问题。商业软件中稀疏脉冲反演是令p=2,q=1,然后用稀疏脉冲反褶积来解此问题。本发明是令p=2,q=0,然后采用匹配追踪方法来解此最优化问题。
具体步骤如图1所示,包括:
输入叠后地震数据,由叠后地震数据进行构造解释;
输入测井数据,提取或给定地震子波,并对叠后地震数据进行标定;同时计算出过井波阻抗:
输入测井数据中的速度和密度,利用合成记录对叠后地震数据进行层位标 定和子波提取。标定过程目前都是通过人机交互手段,人工完成。合成记录的制作是一个简化的一维正演的过程,合成记录F(t)是地震子波S(t)与反射系数R(t)褶积的结果即:F(t)=S(t)*R(t),初始合成地震记录是由反射系数R(t)与标准地震子波S(t)如雷克子波进行褶积得到,为了使得合成记录与地震更匹配,在初始标定的基础上可通过井旁地震道和速度密度测井曲线联合提取子波,子波提取和合成记录制作是一个迭代的过程,经过多次迭代,即能够得到合适的子波和高精度的合成记录;然后对F(t)=S(t)*R(t)进行反褶积,可得到反射系数R(t),再利用反射系数与波阻抗P的关系,即P(n)=[(1+R(n))/(1-R(n))]x P(n-1),可计算出过井波阻抗;;
以构造解释结果作为约束,对过井波阻抗进行内插外推,获得初始波阻抗体,进而获得初始反射系数序列:
这实际上是一个空间插值过程,即利用构造解释得到的某一层位的过井波阻抗值,用插值程序获得该层位上其他点的波阻抗值,然后用其中Ri为各层反射系数序,Zi,Zi-1分别为层位上下的波阻抗值,获得初始反射系数序列;
对于每一道地震数据,由初始反射系数及子波构建目标函数:
借鉴匹配追踪方法的思想而不是数据拟合的思想,即在约束||dobs-W·R||p≤γ下求解min||R||q,其中dobs是观测数据,W为子波,R为反射系数模型,γ为任意小的一个数;
用匹配追踪方法求解目标函数,获得反演后的反射系数:
对于上述最优化问题:在约束||dobs-W·R||p≤γ下求解min||R||q,取p=2,q=1,变成二次规划问题,然后采用拉格朗日乘子等方法来解此问题;
由反演后的反射系数以及标定结果计算最终的波阻抗:
解上述二次规划问题以后,得到反射系数以后,然后用其中,Ri为反演得到各层反射系数序,Z0为第一层的绝对波阻抗,Zn为第n层的绝对波阻抗,获得声波阻抗。
上述技术方案只是本发明的一种实施方式,对于本领域内的技术人员而言,在本发明公开了应用方法和原理的基础上,很容易做出各种类型的改进或变形,而不仅限于本发明上述具体实施方式所描述的方法,因此前面描述的方式只是优选的,而并不具有限制性的意义。
Claims (6)
1.一种基于匹配追踪方法的叠后声波阻抗反演方法,其特征在于:所述方法包括:
S1,输入叠后地震数据,由叠后地震数据进行构造解释;
S2,输入测井数据,提取或给定地震子波,并对叠后地震数据进行标定;同时计算出过井波阻抗:
S3,以构造解释结果作为约束,对S2得到的过井波阻抗进行内插外推,获得初始波阻抗体,进而获得初始反射系数序列;
S4,对于每一道地震数据,由初始反射系数序列及子波构建目标函数;
S5,用匹配追踪方法求解目标函数,获得反演后的反射系数;
S6,由反演后的反射系数以及标定结果计算最终的波阻抗。
2.根据权利要求1所述的基于匹配追踪方法的叠后声波阻抗反演方法,其特征在于:所述S2是这样实现的:
输入测井数据中的速度和密度,利用合成记录对叠后地震数据进行层位标定和子波提取;
合成记录F(t)是地震子波S(t)与反射系数R(t)褶积的结果,即F(t)=S(t)*R(t),初始合成地震记录是由反射系数R(t)与标准地震子波S(t)进行褶积得到;
在初始标定的基础上通过井旁地震道和速度密度测井曲线联合提取子波,子波提取和合成记录制作是一个迭代的过程,经过多次迭代,即得到合适的子波和高精度的合成记录;然后对F(t)=S(t)*R(t)进行反褶积,得到反射系数R(t),再利用反射系数与波阻抗P的关系,即P(n)=[(1+R(n))/(1-R(n))]xP(n-1),计算得到过井波阻抗。
3.根据权利要求2所述的基于匹配追踪方法的叠后声波阻抗反演方法,其特征在于:所述S3是这样实现的:
利用构造解释得到的某一层位的过井波阻抗值,用插值方法获得该层位上其他点的波阻抗值,然后用获得初始反射系数序列。
4.根据权利要求3所述的基于匹配追踪方法的叠后声波阻抗反演方法,其特征在于:所述S4中的目标函数如下:
在约束||dobs-W·R||p≤γ下求解min||R||q,其中dobs是观测数据,W为子波,R为反射系数模型,γ为任意小的一个数。
5.根据权利要求4所述的基于匹配追踪方法的叠后声波阻抗反演方法,其特征在于:所述S5是这样实现的:
在约束||dobs-W·R||p≤γ下求解min||R||q,取p=2,q=1,变成二次规划问题,然后采用拉格朗日乘子等方法来解此问题。
6.根据权利要求5所述的基于匹配追踪方法的叠后声波阻抗反演方法,其特征在于:所述S6是这样实现的:
利用S5得到的反射系数,然后用获得声波阻抗。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510268461.1A CN106291677B (zh) | 2015-05-22 | 2015-05-22 | 一种基于匹配追踪方法的叠后声波阻抗反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510268461.1A CN106291677B (zh) | 2015-05-22 | 2015-05-22 | 一种基于匹配追踪方法的叠后声波阻抗反演方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106291677A true CN106291677A (zh) | 2017-01-04 |
CN106291677B CN106291677B (zh) | 2018-10-02 |
Family
ID=57633256
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510268461.1A Active CN106291677B (zh) | 2015-05-22 | 2015-05-22 | 一种基于匹配追踪方法的叠后声波阻抗反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106291677B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107390269A (zh) * | 2017-07-05 | 2017-11-24 | 西安交通大学 | 一种测井资料统计特性约束的波阻抗反演粒子群优化算法 |
CN107621654A (zh) * | 2017-08-29 | 2018-01-23 | 电子科技大学 | 一种基于最大相关熵的地震叠后波阻抗反演方法 |
CN107894612A (zh) * | 2017-10-23 | 2018-04-10 | 中国地质大学(武汉) | 一种q吸收衰减补偿的声波阻抗反演方法及系统 |
CN108594304A (zh) * | 2018-07-25 | 2018-09-28 | 中国石油化工股份有限公司胜利油田分公司勘探开发研究院 | 基于多极子库利用线性规划求解l1范数的阻抗反演方法 |
CN110749923A (zh) * | 2018-07-24 | 2020-02-04 | 中国石油化工股份有限公司 | 一种基于范数方程提高分辨率的反褶积方法 |
CN112363222A (zh) * | 2020-10-28 | 2021-02-12 | 中国石油天然气集团有限公司 | 叠后自适应宽带约束波阻抗反演方法及装置 |
CN112485826A (zh) * | 2020-11-12 | 2021-03-12 | 中国地质大学(武汉) | 绝对波阻抗反演成像方法、装置、设备及存储介质 |
CN112630835A (zh) * | 2020-12-03 | 2021-04-09 | 重庆三峡学院 | 一种高分辨率的叠后地震波阻抗反演的方法 |
CN113534250A (zh) * | 2020-04-18 | 2021-10-22 | 中国石油化工股份有限公司 | 一种基于快速匹配追踪的多尺度地震反演方法 |
CN113589386A (zh) * | 2021-09-15 | 2021-11-02 | 中国石油大学(北京) | 一种基于对比函数的块状声波阻抗反演方法、装置及设备 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080162049A1 (en) * | 2006-12-28 | 2008-07-03 | Rdsp I, L.P. | Method for interpreting seismic data and controlled source electromagnetic data to estimate subsurface reservoir properties |
CN102393532A (zh) * | 2011-09-06 | 2012-03-28 | 电子科技大学 | 地震信号反演方法 |
CN103293551A (zh) * | 2013-05-24 | 2013-09-11 | 中国石油天然气集团公司 | 一种基于模型约束的阻抗反演方法及系统 |
CN103487835A (zh) * | 2012-06-12 | 2014-01-01 | 中国石油化工股份有限公司 | 一种基于模型约束的多分辨率波阻抗反演方法 |
-
2015
- 2015-05-22 CN CN201510268461.1A patent/CN106291677B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080162049A1 (en) * | 2006-12-28 | 2008-07-03 | Rdsp I, L.P. | Method for interpreting seismic data and controlled source electromagnetic data to estimate subsurface reservoir properties |
CN102393532A (zh) * | 2011-09-06 | 2012-03-28 | 电子科技大学 | 地震信号反演方法 |
CN103487835A (zh) * | 2012-06-12 | 2014-01-01 | 中国石油化工股份有限公司 | 一种基于模型约束的多分辨率波阻抗反演方法 |
CN103293551A (zh) * | 2013-05-24 | 2013-09-11 | 中国石油天然气集团公司 | 一种基于模型约束的阻抗反演方法及系统 |
Non-Patent Citations (3)
Title |
---|
RUI ZHANG ET AL.: "Seismic sparse-layer reflectivity inversion using basis pursuit decomposition", 《GEOPHYSICS》 * |
XIAOTAO WEN ET AL.: "Impedance inversion based on basis pursuit decomposition and contourlet transform", 《SEG DENVER 2014 ANNUAL MEETING》 * |
周东勇等: "双极子分解匹配追踪算法在薄层反演中的应用", 《石油地球物理勘探》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107390269A (zh) * | 2017-07-05 | 2017-11-24 | 西安交通大学 | 一种测井资料统计特性约束的波阻抗反演粒子群优化算法 |
CN107621654A (zh) * | 2017-08-29 | 2018-01-23 | 电子科技大学 | 一种基于最大相关熵的地震叠后波阻抗反演方法 |
CN107894612A (zh) * | 2017-10-23 | 2018-04-10 | 中国地质大学(武汉) | 一种q吸收衰减补偿的声波阻抗反演方法及系统 |
CN110749923A (zh) * | 2018-07-24 | 2020-02-04 | 中国石油化工股份有限公司 | 一种基于范数方程提高分辨率的反褶积方法 |
CN108594304A (zh) * | 2018-07-25 | 2018-09-28 | 中国石油化工股份有限公司胜利油田分公司勘探开发研究院 | 基于多极子库利用线性规划求解l1范数的阻抗反演方法 |
CN113534250A (zh) * | 2020-04-18 | 2021-10-22 | 中国石油化工股份有限公司 | 一种基于快速匹配追踪的多尺度地震反演方法 |
CN112363222A (zh) * | 2020-10-28 | 2021-02-12 | 中国石油天然气集团有限公司 | 叠后自适应宽带约束波阻抗反演方法及装置 |
CN112485826A (zh) * | 2020-11-12 | 2021-03-12 | 中国地质大学(武汉) | 绝对波阻抗反演成像方法、装置、设备及存储介质 |
CN112630835A (zh) * | 2020-12-03 | 2021-04-09 | 重庆三峡学院 | 一种高分辨率的叠后地震波阻抗反演的方法 |
CN112630835B (zh) * | 2020-12-03 | 2023-10-17 | 重庆三峡学院 | 一种高分辨率的叠后地震波阻抗反演的方法 |
CN113589386A (zh) * | 2021-09-15 | 2021-11-02 | 中国石油大学(北京) | 一种基于对比函数的块状声波阻抗反演方法、装置及设备 |
Also Published As
Publication number | Publication date |
---|---|
CN106291677B (zh) | 2018-10-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106291677A (zh) | 一种基于匹配追踪方法的叠后声波阻抗反演方法 | |
CN104516018B (zh) | 一种地球物理勘探中岩性约束下的孔隙度反演方法 | |
CN103713315B (zh) | 一种地震各向异性参数全波形反演方法及装置 | |
US10422922B2 (en) | Method for predicting rock strength by inverting petrophysical properties | |
CN106291682A (zh) | 一种基于基追踪方法的叠后声波阻抗反演方法 | |
CN104570067B (zh) | 一种地球物理勘探中相控地震反演方法 | |
Wu et al. | Deep learning for characterizing paleokarst collapse features in 3‐D seismic images | |
CN104614763A (zh) | 基于反射率法的多波avo储层弹性参数反演方法及系统 | |
CN104122581B (zh) | 一种叠后声波阻抗反演方法 | |
CN107329171A (zh) | 深度域储层地震反演方法及装置 | |
CN104297785A (zh) | 岩相约束储层物性参数反演方法及装置 | |
CN105089652A (zh) | 一种拟声波曲线重构与稀疏脉冲联合反演方法 | |
CN104155687A (zh) | 一种相控叠后声波阻抗反演方法 | |
CN103439740B (zh) | 基于偶极地震子波多重积分的相对阻抗预测的方法及装置 | |
CN105005080A (zh) | 一种利用振幅比率属性识别地层圈闭尖灭线的方法 | |
CN105116449B (zh) | 一种弱反射储层的识别方法 | |
CN105527653A (zh) | 一种基于地质信息的虚拟井构建方法 | |
CN110927791B (zh) | 基于深度学习利用地震数据进行流体预测的方法及装置 | |
CN104991270A (zh) | 一种利用地震正演模拟的叠置砂体识别方法 | |
CN103969685A (zh) | 一种薄互层地震信号的处理方法 | |
CN106338764A (zh) | 生物启发计算地层圈闭油气藏超剥线识别方法 | |
CN109358365A (zh) | 一种基于井震信息匹配的薄层识别方法 | |
US20120099396A1 (en) | System and method for characterization with non-unique solutions of anisotropic velocities | |
CN111273346B (zh) | 去除沉积背景的方法、装置、计算机设备及可读存储介质 | |
Romero et al. | Seeing through the CO2 plume: Joint inversion-segmentation of the Sleipner 4D seismic data set |
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 |