CN104769458A - 一种基于柯西分布的叠后波阻抗反演方法 - Google Patents
一种基于柯西分布的叠后波阻抗反演方法 Download PDFInfo
- Publication number
- CN104769458A CN104769458A CN201480002799.6A CN201480002799A CN104769458A CN 104769458 A CN104769458 A CN 104769458A CN 201480002799 A CN201480002799 A CN 201480002799A CN 104769458 A CN104769458 A CN 104769458A
- Authority
- CN
- China
- Prior art keywords
- data
- seismic
- wave impedance
- impedance
- sampled point
- 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
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. analysis, for interpretation, for correction
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
Abstract
一种基于柯西分布的叠后波阻抗反演方法,采用常规的地震勘探方法采集地震数据对其常规处理得到叠后地震数据;对叠后地震数据进行层位拾取得到层位数据,对确定的目的层层位进行检验和校正以及内插和平滑;采用常规的测井方法得到测井声波时差曲线和密度曲线;根据叠后地震数据、层位数据和已知的钻井分层数据,把深度域的声波时差曲线和密度曲线标定为时间域的曲线,同时生成井中的时间域波阻抗曲线数据,并提取地震子波;生成初始波阻抗模型;对子波数据分别进行归一化;通过反演得到某道的反射系数序列;通过第n个采样点处的相对波阻抗的定义,通过推导得到波阻抗与反射系数的关系对所有的地震道重复上述过程,得到所有道的波阻抗反演结果。
Description
技术领域
本发明涉及地球物理勘探技术,属于储层预测反演技术类,是一种基于柯西分布的叠后波阻抗反演方法。
背景技术
地震勘探是通过人工激发地震波,在地表或地下用单分量或多分量传感器记录地震波的地层响应,研究它们在地层中的传播规律,通过地震数据处理及反演等方法以查明地下的地质构造岩性特征,进而寻找油气等矿产资源的地球物理勘探方法。地震勘探从认识地下的构造形态开始,发展到直接应用地震信息判断岩性、分析岩相、定量计算岩层的物性参数及寻找油气显示等。地震反演技术正是这一发展过程的产物。
地震反演的基本目的是利用地震波在地下介质中的传播规律,通过数据采集、处理与解释等流程,来推测地下岩层结构和物性参数的空间分布,为勘探开发提供重要依据。在地震反演研究中,有多种参数的反演方法,如波阻抗、速度、密度、孔隙度、渗透率、泊松比等。由于波阻抗信息是联系地质和地球物理的一座桥梁,在叠后计算数据量相对要小,在实际生产中应用方便而且效果明显,因此波阻抗反演在地震反演中具有特殊的地位,地震反演通常是指波阻抗反演。
常规的地震波阻抗反演就是指利用地震贤料反演地层/岩层波阻抗的地震特殊处理技术。与地震多参数模式识别预测储层油气、神经网络预测地层参数、振幅拟合预测储层厚度等统计性方法相比,波阻抗反演具有明确的物理意义,它是储层岩性预测、油藏特征描述的确定性方法,在实际应用中取得了显著的地质效果。
目前的反演方法多数是以模型为基础的方法,这些方法一般都依据测井及地质资料建立初始模型,通过广义线性反演方法进行迭代求取岩性参数。由于该问题的非线性,所以除了要求精细的子波外,还要求初始模型接近真实模型,才能达到可靠的结果,即反演结果强烈依赖于初始模型的选择。除此类方法外,全局优化的反演方法(如遗传算法和模拟退火算法等)虽然克服了基于模型方法对初始模型依赖性强的缺陷,但是由于其得到全局最优的反演结果,所以反演速度很慢。
发明内容
本发明目的是提供一种基于柯西分布的叠后波阻抗反演方法。
本发明通过如下技术方案实现:
1)采用常规的地震勘探方法采集地震数据,对地震数据进行常规处理得到叠后地震数据;
2)对叠后地震数据进行层位拾取得到层位数据,对确定的目的层层位进行检验和校正以及内插和平滑;
3)采用常规的测井方法得到测井数据,得到测井声波时差曲线和密度曲线;
4)根据叠后地震数据、层位数据和已知的钻井分层数据,把深度域的声波时差曲线和密度曲线标定为时间域的曲线,同时生成井中的时间域波阻抗曲线数据,并在井旁地震道上提取地震子波;
所述的标定为利用测井曲线和地震子波模拟井旁地震记录,实现测井分层到地震层位的标定和映射,由此得到时深关系曲线,由此时深关系可以将深度域的测井曲线转换为时间域曲线。
5)利用步骤2)的层位数据和步骤4)得到的时间域波阻抗曲线,生成初始波阻抗模型;
6)对地震数据和步骤4)得到的子波数据分别进行归一化,归一到范围[-1,1]之间;
7)把步骤1)读入的一道地震数据、步骤5)生成的初始波阻抗模型的一道数据和步骤4)提取的子波数据,代入以下公式中,通过反演得到该道的反射系数序列:
r=(GTG+μQ+ρCTC)-1(GTd+ρCTξ) (1)
其中d=[d1,d2,...,dN]T是叠后地震数据,N为地震数据的总采样点数;r=[r1,r2,...,rN]T是反射系数序列;G是N×N维子波褶积矩阵,上角标T代表矩阵的转置;μ是稀疏约束因子,控制反射系数的稀疏程度;矩阵Q的对角元素为:Qnn表示矩阵Q的第n行第n列元素的值,n为Q矩阵的行列号,矩阵Q除对角元素外其余元素均为零,σr代表噪声的标准方差,Rn是由初始波阻抗模型计算得到的第n个采样点位置处的初始反射系数;ρ为模型约束因子,控制反演结果对初始模型的依赖程度;C为积分算子矩阵,其离散形式表示为:
公式(1)中上角标-1为对矩阵求逆;为第n个采样点处的相对波阻抗值,I0为反演时窗内第一个采样点对应的初始波阻抗值,In为反演时窗内第n个采样点处的初始波阻抗值,ln为自然对数符号,ri为反演时窗内第i个采样点的反射系数值,表示对从第1个采样点到第n个采样点的ri进行求和计算。
8)通过步骤7)中第n个采样点处的相对波阻抗的定义,通过推导可以得到波阻抗与反射系数的关系:
其中:In为反演时窗内第n个采样点的波阻抗值,I0为反演时窗内第一个采样点对应的初始波阻抗值,ri为步骤7)中反演得到的第i个采样点反射系数值,e表示自然对数的底,表示对从第1个采样点到第n个采样点的ri进行求和计算。
通过公式(2)将步骤7)反演得到的一道的反射系数系列转换为波阻抗序列,便可得到该道的波阻抗反演结果;
9)对所有的地震道重复步骤7)至8)过程,得到所有道的波阻抗反演结果。
本发明具有如下特点,主要表现为:
(1)通过假设反射系数的先验概率服从柯西分布来实现,柯西分布属于长尾巴分布,相对于高斯分布来说其峰值处分布更窄,且逼近零的速度也较缓,从而得到少量的非零值和大量的零值来实现稀疏脉冲反演;
(2)通过添加波阻抗模型约束可以控制反演结果的准确性和稳定性;
(3)稀疏脉冲反演是以地震道为主的反演方法,反演结果的分辨率、信噪比以及可靠程度主要依赖于地震资料本身的品质,地震噪音对反演结果敏感,且影响大,因此用于稀疏脉冲反演的地震资料应具有较宽的频带、较低的噪声、相对振幅保持和成像准确等特征。测井资料,尤其是声波测井和密度测井资料,是地震横向预测的对比标准和解释依据,在反演处理之前应进行仔细的编辑和校正,使其能够正确反映岩层的物理特征;
(4)稀疏脉冲反演得到的波阻抗结果在井位处与井曲线吻合度高,且运算效率高。
本发明的反演方法不仅具有一般递推反演方法的特点,即反演结果忠实于地震资料,能反映储层的横向变化。而且,在迭代过程中引入地质和测井资料参与反演约束,增加了部分低频和高频成分,一定程度拓宽了反演频带。该方法对初始模型依赖较小,反演结果的唯一性较好,不易出现假象。
附图说明
图1本发明的应用实例1模型反演结果的叠后地震剖面;
图2是本发明的应用实例1模型反演结果的波阻抗剖面;
图3是本大港岐口地区本发明反演波阻抗结果与Jason软件反演结果对比。
具体实施方式
本发明所提供的方法,是一种基于柯西分布的叠后波阻抗反演方法,而且运算效率非常高。
本发明通过如下技术方案实现:
1)采用常规的地震勘探方法采集地震数据,对地震数据进行常规处理得到叠后地震数据;
2)对叠后地震数据进行层位拾取得到层位数据,对确定的目的层层位进行检验和校正以及内插和平滑;
3)采用常规的测井方法得到测井数据,得到测井声波时差曲线和密度曲线;
4)根据叠后地震数据、层位数据和已知的钻井分层数据,把深度域的声波时差曲线和密度曲线标定为时间域的曲线,同时生成井中的时间域波阻抗曲线数据,并在井旁地震道上提取地震子波;
所述的标定为利用测井曲线和地震子波模拟井旁地震记录,实现测井分层到地震层位的标定和映射,由此得到时深关系曲线,由此时深关系可以将深度域的测井曲线转换为时间域曲线。
5)利用步骤2)的层位数据和步骤4)得到的时间域波阻抗曲线,生成初始波阻抗模型;
6)对地震数据和步骤4)得到的子波数据分别进行归一化,归一到范围[-1,1]之间;
7)把步骤1)读入的一道地震数据、步骤5)生成的初始波阻抗模型的一道数据和步骤4)提取的子波数据,代入以下公式中,通过反演得到该道的反射系数序列:
r=(GTG+μQ+ρCTC)-1(GTd+ρCTξ) (1)
其中d=[d1,d2,...K,dN]T是叠后地震数据,N为地震数据的总采样点数;r=[r1,r2,...,rN]T是反射系数序列;G是N×N维子波褶积矩阵,上角标T代表矩阵的转置;μ是稀疏约束因子,控制反射系数的稀疏程度;矩阵Q的对角元素为:Qnn表示矩阵Q的第n行第n列元素的值,n为Q矩阵的行列号,矩阵Q除对角元素外其余元素均为零,σr代表噪声的标准方差,Rn是由初始波阻抗模型计算得到的第n个采样点位置处的初始反射系数;ρ为模型约束因子,控制反演结果对初始模型的依赖程度;C为积分算子矩阵,其离散形式表示为:
公式(1)中上角标-1为对矩阵求逆;为第n个采样点处的相对波阻抗值,I0为反演时窗内第一个采样点对应的初始波阻抗值,In为反演时窗内第n个采样点处的初始波阻抗值,ln为自然对数符号,ri为反演时窗内第i个采样点的反射系数值,表示对从第1个采样点到第n个采样点的ri进行求和计算。
8)通过步骤7)中第n个采样点处的相对波阻抗的定义,通过推导可以得到波阻抗与反射系数的关系:
其中:In为反演时窗内第n个采样点的波阻抗值,I0为反演时窗内第一个采样点对应的初始波阻抗值,ri为步骤7)中反演得到的第i个采样点反射系数值,e表示自然对数的底,表示对从第1个采样点到第n个采样点的ri进行求和计算。
通过公式(2)将步骤7)反演得到的一道的反射系数系列转换为波阻抗序列,便可得到该道的波阻抗反演结果;
9)对所有的地震道重复步骤7)至8)过程,得到所有道的波阻抗反演结果。
以下通过具体实例说明本发明的效果。
发明应用实施例1:
使用本专利方法对一个二维模型数据(来自Jason公司)进行试验,该模型垂直方向是三个菱形砂体,每个砂体内部从上到下所含流体分别为气、油和水。图1为该模型的无噪音叠后地震剖面,依据前面的步骤1)-9)进行基于柯西分布的稀疏脉冲反演,得到图2的波阻抗剖面。由图2可以看出反演结果在井位处与井曲线吻合的很好,并且从图上可以清晰的反映出砂体内部的阻抗差异,说明反演结果很准确。
反演时窗为1250毫秒,采样间隔为4毫秒,共200道,用本发明在HP2单机上进行反演耗时约110秒。
发明应用实施例2:
使用本专利方法对大港岐口地区进行反演,依据前面的步骤1)-9)进行基于柯西分布的稀疏脉冲反演,图3(a)为本发明得到的波阻抗剖面,图3(b)为使用国外软件Jason系统进行反演得到的波阻抗剖面,从图中可看出总体上本发明反演结果与Jason系统反演结果相当。
反演目的层时窗约900毫秒,采样间隔为4毫秒,共300道,用本发明在HP2单机上进行反演耗时约340秒。
理论模型和实际资料的结果表明,柯西分布不仅能较好地恢复反射系数序列的稀疏性,同时能在提高地震资料的分辨率和减小对弱反射信息的压制之间达到一个平衡。由该方法得到的反演结果不仅准确,精度高,而且运算效率很高,利用此结果可以进行准确的储层预测。
Claims (1)
1.一种基于柯西分布的叠后波阻抗反演方法,步骤如下:
1)采用常规的地震勘探方法采集地震数据,对地震数据进行常规处理得到叠后地震数据;
2)对叠后地震数据进行层位拾取得到层位数据,对确定的目的层层位进行检验和校正以及内插和平滑;
3)采用常规的测井方法得到测井数据,得到测井声波时差曲线和密度曲线;
4)根据叠后地震数据、层位数据和已知的钻井分层数据,把深度域的声波时差曲线和密度曲线标定为时间域的曲线,同时生成井中的时间域波阻抗曲线数据,并在井旁地震道上提取地震子波;
步骤4)所述的标定为利用测井曲线和地震子波模拟井旁地震记录,实现测井分层到地震层位的标定和映射,由此得到时深关系曲线,由此时深关系可以将深度域的测井曲线转换为时间域曲线。
5)利用步骤2)的层位数据和步骤4)得到的时间域波阻抗曲线,生成初始波阻抗模型;
6)对地震数据和步骤4)得到的子波数据分别进行归一化,归一到范围[-1,1]之间;
7)把步骤1)读入的一道地震数据、步骤5)生成的初始波阻抗模型的一道数据和步骤4)提取的子波数据,代入以下公式中,通过反演得到该道的反射系数序列:
r=(GTG+μQ+ρCTC)-1(GTd+ρCTξ) (1)
其中d=[d1,d2,...,dN]T是叠后地震数据,N为地震数据的总采样点数;r=[r1,r2,...,rN]T是反射系数序列;G是N×N维子波褶积矩阵,上角标T代表矩阵的转置;μ是稀疏约束因子,控制反射系数的稀疏程度;矩阵Q的对角元素为:Qnn表示矩阵Q的第n行第n列元素的值,n为Q矩阵的行列号,矩阵Q除对角元素外其余元素均为零,σr代表噪声的标准方差,Rn是由初始波阻抗模型计算得到的第n个采样点位置处的初始反射系数;ρ为模型约束因子,控制反演结果对初始模型的依赖程度;C为积分算子矩阵,其离散形式表示为:
公式(1)中上角标-1为对矩阵求逆;为第n个采样点处的相对波阻抗值,I0为反演时窗内第一个采样点对应的初始波阻抗值,In为反演时窗内第n个采样点处的初始波阻抗值,ln为自然对数符号,ri为反演时窗内第i个采样点的反射系数值,表示对从第1个采样点到第n个采样点的ri进行求和计算;
8)通过步骤7)中第n个采样点处的相对波阻抗的定义,通过推导得到波阻抗与反射系数的关系:
其中:In为反演时窗内第n个采样点的波阻抗值,I0为反演时窗内第一个采样点对应的初始波阻抗值,ri为步骤7)中反演得到的第i个采样点反射系数值,e表示自然对数的底,表示对从第1个采样点到第n个采样点的ri进行求和计算;
通过公式(2)将步骤7)反演得到的一道的反射系数系列转换为波阻抗序列,便可得到该道的波阻抗反演结果;
9)对所有的地震道重复步骤7)至8)过程,得到所有道的波阻抗反演结果。
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/CN2014/082265 WO2016008105A1 (zh) | 2014-07-15 | 2014-07-15 | 一种基于柯西分布的叠后波阻抗反演方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104769458A true CN104769458A (zh) | 2015-07-08 |
Family
ID=53649880
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201480002799.6A Pending CN104769458A (zh) | 2014-07-15 | 2014-07-15 | 一种基于柯西分布的叠后波阻抗反演方法 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN104769458A (zh) |
WO (1) | WO2016008105A1 (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105223617A (zh) * | 2015-09-21 | 2016-01-06 | 中国石油化工股份有限公司江汉油田分公司物探研究院 | 利用地震属性控制初始波阻抗的建模方法及系统 |
CN105954803A (zh) * | 2016-07-08 | 2016-09-21 | 中国石油天然气股份有限公司 | 叠后地震反演方法及装置 |
CN106443770A (zh) * | 2015-08-07 | 2017-02-22 | 中国石油化工股份有限公司 | 一种页岩气地质甜点的预测方法 |
CN108535775A (zh) * | 2018-03-30 | 2018-09-14 | 中国石油大学(北京) | 非平稳地震资料声波阻抗反演方法及装置 |
CN110954953A (zh) * | 2019-12-16 | 2020-04-03 | 中国地质大学(武汉) | 一种基于柯西分布的叠后波阻抗反演方法 |
CN111427085A (zh) * | 2020-04-01 | 2020-07-17 | 中国石油天然气股份有限公司 | 一种碳酸盐岩层间岩溶储层预测方法及装置 |
CN112558158A (zh) * | 2019-09-25 | 2021-03-26 | 中国石油化工股份有限公司 | 一种基于测井曲线的地震数据双向拓频方法及系统 |
CN113589386A (zh) * | 2021-09-15 | 2021-11-02 | 中国石油大学(北京) | 一种基于对比函数的块状声波阻抗反演方法、装置及设备 |
CN117310802A (zh) * | 2023-09-13 | 2023-12-29 | 成都捷科思石油天然气技术发展有限公司 | 一种深度域储层地震反演方法 |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108072901B (zh) * | 2016-11-10 | 2020-12-01 | 中国石油化工股份有限公司 | 一种获得准确的静水压力和上覆地层压力的方法及系统 |
CN112363222A (zh) * | 2020-10-28 | 2021-02-12 | 中国石油天然气集团有限公司 | 叠后自适应宽带约束波阻抗反演方法及装置 |
CN113296150B (zh) * | 2021-07-02 | 2022-06-10 | 清华大学 | 一种测井约束下的高维闭环网络地震反演方法 |
CN115494547B (zh) * | 2022-10-21 | 2023-04-28 | 成都理工大学 | 基于对数全变分稀疏约束的地震波阻抗反演方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6324477B1 (en) * | 2000-03-01 | 2001-11-27 | Apache Corporation | System for processing well log data |
CN1797037A (zh) * | 2004-12-29 | 2006-07-05 | 中国石油天然气集团公司 | 一种地震波波阻抗反演的方法 |
US20100296367A1 (en) * | 2009-05-20 | 2010-11-25 | Florence Delprat-Jannaud | Method of imaging a target area of the subsoil from walkaway type data |
CN103792571A (zh) * | 2012-10-26 | 2014-05-14 | 中国石油化工股份有限公司 | 点约束贝叶斯稀疏脉冲反演方法 |
CN104280767A (zh) * | 2013-07-12 | 2015-01-14 | 中国石油天然气集团公司 | 一种基于柯西分布的稀疏脉冲反演方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005004335A2 (en) * | 2003-06-25 | 2005-01-13 | Georgia Tech Research Corporation | Cauchy-distribution based coding system and method |
CN101520516B (zh) * | 2008-02-25 | 2012-02-29 | 中国石油集团东方地球物理勘探有限责任公司 | 一种三维地震记录的叠前似二维变换的方法 |
CN103376464B (zh) * | 2012-04-13 | 2016-04-06 | 中国石油天然气集团公司 | 一种地层品质因子反演方法 |
US9581707B2 (en) * | 2012-12-04 | 2017-02-28 | OOO Rock Flow Dynamics | System for automated identification of surfaces for building of geologic hydrodynamic model of oil and gas deposit by seismic data |
-
2014
- 2014-07-15 WO PCT/CN2014/082265 patent/WO2016008105A1/zh active Application Filing
- 2014-07-15 CN CN201480002799.6A patent/CN104769458A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6324477B1 (en) * | 2000-03-01 | 2001-11-27 | Apache Corporation | System for processing well log data |
CN1797037A (zh) * | 2004-12-29 | 2006-07-05 | 中国石油天然气集团公司 | 一种地震波波阻抗反演的方法 |
US20100296367A1 (en) * | 2009-05-20 | 2010-11-25 | Florence Delprat-Jannaud | Method of imaging a target area of the subsoil from walkaway type data |
CN103792571A (zh) * | 2012-10-26 | 2014-05-14 | 中国石油化工股份有限公司 | 点约束贝叶斯稀疏脉冲反演方法 |
CN104280767A (zh) * | 2013-07-12 | 2015-01-14 | 中国石油天然气集团公司 | 一种基于柯西分布的稀疏脉冲反演方法 |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106443770A (zh) * | 2015-08-07 | 2017-02-22 | 中国石油化工股份有限公司 | 一种页岩气地质甜点的预测方法 |
CN105223617A (zh) * | 2015-09-21 | 2016-01-06 | 中国石油化工股份有限公司江汉油田分公司物探研究院 | 利用地震属性控制初始波阻抗的建模方法及系统 |
CN105954803A (zh) * | 2016-07-08 | 2016-09-21 | 中国石油天然气股份有限公司 | 叠后地震反演方法及装置 |
CN105954803B (zh) * | 2016-07-08 | 2018-02-02 | 中国石油天然气股份有限公司 | 叠后地震反演方法及装置 |
CN108535775A (zh) * | 2018-03-30 | 2018-09-14 | 中国石油大学(北京) | 非平稳地震资料声波阻抗反演方法及装置 |
CN108535775B (zh) * | 2018-03-30 | 2019-08-23 | 中国石油大学(北京) | 非平稳地震资料声波阻抗反演方法及装置 |
CN112558158A (zh) * | 2019-09-25 | 2021-03-26 | 中国石油化工股份有限公司 | 一种基于测井曲线的地震数据双向拓频方法及系统 |
CN112558158B (zh) * | 2019-09-25 | 2024-04-09 | 中国石油化工股份有限公司 | 一种基于测井曲线的地震数据双向拓频方法及系统 |
CN110954953A (zh) * | 2019-12-16 | 2020-04-03 | 中国地质大学(武汉) | 一种基于柯西分布的叠后波阻抗反演方法 |
CN111427085B (zh) * | 2020-04-01 | 2023-02-10 | 中国石油天然气股份有限公司 | 一种碳酸盐岩层间岩溶储层预测方法及装置 |
CN111427085A (zh) * | 2020-04-01 | 2020-07-17 | 中国石油天然气股份有限公司 | 一种碳酸盐岩层间岩溶储层预测方法及装置 |
CN113589386A (zh) * | 2021-09-15 | 2021-11-02 | 中国石油大学(北京) | 一种基于对比函数的块状声波阻抗反演方法、装置及设备 |
CN117310802A (zh) * | 2023-09-13 | 2023-12-29 | 成都捷科思石油天然气技术发展有限公司 | 一种深度域储层地震反演方法 |
Also Published As
Publication number | Publication date |
---|---|
WO2016008105A1 (zh) | 2016-01-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104769458A (zh) | 一种基于柯西分布的叠后波阻抗反演方法 | |
Chadwick et al. | Quantitative analysis of time-lapse seismic monitoring data at the Sleipner CO 2 storage operation | |
Bosch et al. | Seismic inversion for reservoir properties combining statistical rock physics and geostatistics: A review | |
Vardy et al. | State‐of‐the‐art remote characterization of shallow marine sediments: the road to a fully integrated solution | |
CN106368691A (zh) | 基于岩石物理地震信息三维异常孔隙压力预测方法 | |
CN104570067A (zh) | 一种地球物理勘探中相控地震反演方法 | |
CN104516018A (zh) | 一种地球物理勘探中岩性约束下的孔隙度反演方法 | |
CN105277982A (zh) | 一种泥页岩总有机碳含量地震预测方法 | |
CN104459777A (zh) | 基于流体体积模量avo反演的流体识别的方法及系统 | |
CN104200115A (zh) | 一种基于地质统计学模拟的全地层速度建模方法 | |
Janson et al. | 3-D forward seismic model of an outcrop-based geocellular model | |
Faraji et al. | Integration of core data, well logs and seismic attributes for identification of the low reservoir quality units with unswept gas in the carbonate rocks of the world’s largest gas field | |
Oyeyemi et al. | Hydrocarbon resource evaluation using combined petrophysical analysis and seismically derived reservoir characterization, offshore Niger Delta | |
Li et al. | An Integrated quantitative modeling approach for fault-related fractures in tight sandstone reservoirs | |
CN104280767A (zh) | 一种基于柯西分布的稀疏脉冲反演方法 | |
Sanda et al. | The integrated approach to seismic attributes of lithological characterization of reservoirs: case of the F3 Block, North Sea-Dutch Sector | |
CN112965106B (zh) | 利用纵横波速度的奇异性识别天然气水合物赋存类型方法 | |
Rotimi et al. | Reservoir characterization and modeling of lateral heterogeneity using multivariate analysis | |
CN113514884A (zh) | 一种致密砂岩储层预测方法 | |
Varela et al. | On the value of 3D seismic amplitude data to reduce uncertainty in the forecast of reservoir production | |
CN113806674A (zh) | 古河道纵向尺度的量化方法、装置、电子设备及存储介质 | |
CN112147676A (zh) | 一种煤层及夹矸厚度预测方法 | |
Naseer et al. | Seismic attributes and spectral decomposition-based inverted porosity-constrained simulations for appraisal of shallow-marine lower-Cretaceous sequences of Miano gas field, Southern Pakistan | |
Naseer et al. | Delineation of stratigraphic traps within the basin floor fans of Miocene sedimentary sequences, offshore Indus, Pakistan using inverted acoustic impedance simulations | |
Lacazette et al. | Passive Seismic Methods for Unconventional Resource Development |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20150708 |