CN110954953A - 一种基于柯西分布的叠后波阻抗反演方法 - Google Patents
一种基于柯西分布的叠后波阻抗反演方法 Download PDFInfo
- Publication number
- CN110954953A CN110954953A CN201911290541.1A CN201911290541A CN110954953A CN 110954953 A CN110954953 A CN 110954953A CN 201911290541 A CN201911290541 A CN 201911290541A CN 110954953 A CN110954953 A CN 110954953A
- Authority
- CN
- China
- Prior art keywords
- data
- inversion
- wave impedance
- reflection coefficient
- curve
- 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 45
- 238000001228 spectrum Methods 0.000 claims abstract description 22
- 238000012545 processing Methods 0.000 claims abstract description 19
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 7
- 238000005070 sampling Methods 0.000 claims description 4
- 101100356682 Caenorhabditis elegans rho-1 gene Proteins 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 239000010410 layer Substances 0.000 description 22
- 239000004576 sand Substances 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 230000000994 depressogenic effect Effects 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 239000011435 rock Substances 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 239000011229 interlayer Substances 0.000 description 1
- 238000004519 manufacturing process Methods 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
- 238000003909 pattern recognition Methods 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 238000006467 substitution reaction Methods 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
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- 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
- G01V1/30—Analysis
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
本发明公开了一种基于柯西分布的叠后波阻抗反演方法,采用常规的地震勘探方法采集地震数据,对地震数据进行常规处理之后得到叠后地震数据,然后对叠后地震数据进行进行处理得到层位数据,本发明不仅具有一般递推反演方法的特性,而且相对于一般递推反演方法得到的结果较为准确,不易出现假象,首先利用柯西分布对其数据进行处理,其次通过奇偶分解原理进行数据处理,将谱反演和波阻抗递推反演相结合,然后利用谱反演准确求取反射系数,此外在有井的情况下,以井为控制,在地层等厚且产状平缓时,进行简单的线性内插确定该点的低频分量,在无井的情况下,将反褶积得到的反射系数转为波阻抗,从而得到该点精准的低频分量,使得反演结果更加的精确。
Description
技术领域
本发明公开了一种基于柯西分布的叠后波阻抗反演方法,属于波阻抗反演技术领域。
背景技术
随着科学技术的进步,地震勘探技术逐步发展,地震反演技术正是这一发展过程的产物,地震反演的基本日的是利用地震波在地下介质中的传播规律,通过数据采集、处理与解释等流程,来推测地下岩层结构和物性参数的空间分布,为勘探开发提供重要依据,而常规的地表波阻抗反演就是指利用地震贤料反演地层/岩层波阻抗的地震特殊处理技术,与地震多参数模式识别预测储层油气、神经网络预测地层参数、振幅拟合预测储层厚度等统计性方法相比,波阻抗反演具有明确的物理总义,它是储层岩性预测、汕藏特征描述的确定性方法,在实际应用中取得了显著的地质效果。
但是目前的递推反演方法由于子波的混合相位因素,会使反射剖面的极性发生变化,而地震反射波的极性是正还是负直接影响到反演波阻抗后速度变高还是变低,同时地震反演中对反射系数的标定,通常是根据井中反射系数来标定反褶积后的振幅值的,但是,求波阻抗是一个积分的过程,反褶积后的地震道振幅实际上还不是反射系数,而是相当于反射系数再褶积一个剩余子波,这个剩余子波一般在浅层主频高些,深层主频低些,频率低的波积分后数值偏大,会使深层产生偏大的波阻抗值,此外在有井的情况下,以井为控制,能够得到该点的低频分量,但是井与井之间低频分量的内插又是一个难题,简单的线性内插无法得到该点精准的低频分量,在无井区,作波阻抗反演往往要从叠加速度谱中提取低频分量,由于速度谱的质量和分辨率问题,这些问题解决的好坏直接影响着地震反演结果的可靠性。
发明内容
本发明要解决的技术问题是克服传统方法的技术缺陷,提供一种基于柯西分布的叠后波阻抗反演方法,从而解决上述问题。
为实现上述目的,本发明提供如下技术方案:一种基于柯西分布的叠后波阻抗反演方法,步骤如下:
(1)采用常规的地震勘探方法采集地震数据,对地震数据进行常规处理之后得到叠后地震数据,然后对叠后地震数据进行进行处理得到层位数据;
(2)采用常规的测井方法得到测井数据,将测井数据进行处理得到测井声波时差曲线和密度曲线,然后通过分析测井声波时差曲线和密度曲线得到叠加速度谱质量曲线和分辨率曲线。
(3)根据叠后地震数据、层位数据以及测井数据实现测井分层到地震层位的标定和映射,由此得到时深关系曲线,然后根据测井曲线与时深关系曲线得到时间域曲线,接着根据时间域波抗曲线,生成初始波阻抗模型。
(4)把所述步骤(1)中的叠后地震数据、步骤(3)中生成的初始波阻抗模型一道数据和步骤(2)中得到的声波密度数据,代入以下公式,然后根据波的传播理论,当平面弹性波垂直入射到两种介质的平面分界面上时,通过反演得到该道的反射系数:
其中,ρ1和ρ2是界面两侧介质的密度;V1和V2是界面两侧介质的波速;Z=ρV是界面两侧介质的波阻抗;γ是界面反射系数;
当地面下的介质为层状介质,并存在着一系列互相平行的N个反射界面时,对于其中第i个界面,当波垂直入射时,其反射系数为:
当r、Z1已知时,可以由(1式)导出求Z2的公式
当r、Z2已知时,则:
(5)然后将所述步骤(4)中得到的反射系数数据,根据奇偶分解原理,将反射系数分解为奇偶分量,即:
r(t)=r0(t)+re(t) (5式)
其中,r0(t)=【r(t)-r(-t)】/2为奇分量;
re(t)=【r(t)+r(-t)】/2为偶分量;
然后建立频率域反射系数反演目标函数,通过最小平方法得到各个采样点的反射系数,最后利用谱反演准确求取反射系数。
(6)在所述步骤(5)中得到的数据基础上,当地下存在一系列互相平行的反射界面时,已知第N层界面以上N-1层界面的反射系数ri(i=1,2,...N-1)和第一层介质的波阻抗Z1,由(3式)可导出任意第n层介质的波阻抗Zn,得到:
(7)而已知第n层介质波阻抗Zn时,由(4式)又可导出任意第m层(m<n)介质的Zm,即:
(8)利用(6式)或(7式),可以从声波时差曲线及密度曲线上(没有密度测井时可利用Gardnar公式,ρ=0.31V0.25换算)选择标准层波阻抗作为基准波阻抗,将反褶积得到的反射系数转换为波阻抗。
(9)对所有的地震道重复所述步骤(6)至步骤(8)过程得到所有道的反演结果。
本发明具有如下特点,主要表现为:
(1)对于子波的混合相位进行分部逐一通过服从柯西分布进行处理,单次进行,避免由于子波的混合相位因素造成反射剖面的极性发生变化,从而避免直接影响到反演波阻抗后速度变高还是变低。
(2)通过常规的地震勘探方法多次采集地震数据,然后对于采集到的地震数据进行反复多次处理,通过公式计算出反射系数,然后将多次反射系数根据奇偶分解原理进行处理,同时除了考虑时变的振幅因素外,还要考虑时变的主频变化,将谱反演和波阻抗递推反演相结合,然后利用谱反演准确求取反射系数。
(3)在有井的情况下,以井为控制,在地层等厚且产状平缓时,进行简单的线性内插,从而保证能够得到该点精准的低频分量,在无井的情况下,需保证速度谱的质量和分辨率,然后才能从叠加速度谱中提取低频分量,从而更好的保证作波阻抗反演中低频分量的准确无误。
与现有技术相比,本发明的有益效果如下:
本发明与现有技术相比,不仅具有一般递推反演方法的特性,即可以得到地层的绝对波阻抗,分辨率也比道积分剖面高,而且相对于一般递推反演方法得到的结果较为准确,不易出现假象,首先在子波的混合相位中利用柯西分布对其数据进行处理,避免由于子波的混合相位因素造成反射剖面的极性发生变化,其次通过奇偶分解原理进行数据处理,将谱反演和波阻抗递推反演相结合,然后利用谱反演准确求取反射系数,此外在有井的情况下,以井为控制,在地层等厚且产状平缓时,进行简单的线性内插确定该点的低频分量,同时,在无井的情况下,通过在速度谱的质量曲线与速度谱的分辨率曲线上选择标准层波阻抗作为基准波阻抗,将反褶积得到的反射系数转为波阻抗,从而得到该点精准的低频分量,使得反演结果更加的精确。
具体实施方式
下面对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
一种基于柯西分布的叠后波阻抗反演方法,步骤如下:
(1)采用常规的地震勘探方法采集地震数据,对地震数据进行常规处理之后得到叠后地震数据,然后对叠后地震数据进行进行处理得到层位数据;
(2)采用常规的测井方法得到测井数据,将测井数据进行处理得到测井声波时差曲线和密度曲线,然后通过分析测井声波时差曲线和密度曲线得到叠加速度谱质量曲线和分辨率曲线,在对于子波的混合相位进行处理时,进行分部逐一通过服从柯西分布进行处理,单次进行,避免由于子波的混合相位因素造成反射剖面的极性发生变化,从而避免直接影响到反演波阻抗后速度变高还是变低。
(3)根据叠后地震数据、层位数据以及测井数据实现测井分层到地震层位的标定和映射,由此得到时深关系曲线,然后根据测井曲线与时深关系曲线得到时间域曲线,接着根据时间域波抗曲线,生成初始波阻抗模型。
(4)把步骤(1)中的叠后地震数据、步骤(3)中生成的初始波阻抗模型一道数据和步骤(2)中得到的声波密度数据,代入以下公式,然后根据波的传播理论,当平面弹性波垂直入射到两种介质的平面分界面上时,通过反演得到该道的反射系数:
其中,ρ1和ρ2是界面两侧介质的密度;V1和V2是界面两侧介质的波速;Z=ρV是界面两侧介质的波阻抗;γ是界面反射系数;
当地面下的介质为层状介质,并存在着一系列互相平行的N个反射界面时,对于其中第i个界面,当波垂直入射时,其反射系数为:
当r、Z1已知时,可以由(1式)导出求Z2的公式
当r、Z2已知时,则:
(5)然后将步骤(4)中得到的反射系数数据,根据奇偶分解原理,将反射系数分解为奇偶分量,即:
r(t)=r0(t)+re(t) (5式)
其中,r0(t)=【r(t)-r(-t)】/2为奇分量;
re(t)=【r(t)+r(-t)】/2为偶分量;
然后建立频率域反射系数反演目标函数,通过最小平方法得到各个采样点的反射系数,同时除了考虑时变的振幅因素外,还要考虑时变的主频变化,将谱反演和波阻抗递推反演相结合,最后利用谱反演准确求取反射系数。
(6)在步骤(5)中得到的数据基础上,当地下存在一系列互相平行的反射界面时,已知第N层界面以上N-1层界面的反射系数ri(i=1,2,...N-1)和第一层介质的波阻抗Z1,由(3式)可导出任意第n层介质的波阻抗Zn,得到:
(7)而已知第n层介质波阻抗Zn时,由(4式)又可导出任意第m层(m<n)介质的Zm,即:
(8)利用(6式)或(7式),可以从声波时差曲线及密度曲线上(没有密度测井时可利用Gardnar公式,ρ=0.31V0.25换算)选择标准层波阻抗作为基准波阻抗,将反褶积得到的反射系数转换为波阻抗,同时,在有井的情况下,以井为控制,在地层等厚且产状平缓时,进行简单的线性内插,从而保证能够得到该点精准的低频分量,在无井的情况下,需保证速度谱的质量和分辨率,然后才能从叠加速度谱中提取低频分量,从而更好的保证作波阻抗反演中低频分量的准确无误。
(9)对所有的地震道重复步骤(6)至步骤(8)过程得到所有道的反演结果。
以下通过具体实施例说明本发明的有益效果:
发明应用实施例1:
使用本专利方法对济阳坳陷地区进行反演,河流相储集体主要发育在上第三系,地层较为平缓,断裂***不发育,存在的主要问题是“迷宫状”的河道无法准确的建立三维初始模型,同时由于地理因素以及子波的混合相位因素的影响可能会使得反射剖面的极性发生变化,从而直接影响到反演波阻抗后速度变高还是变低,而一般的多采用对模型精度要求低的反演方法,进而使得得到的反演结果拥有较大的误差,容易出现假象,但是依据前面的步骤(1)至(9)对济阳坳陷地区进行反演,通过对于子波的混合相位进行分部逐一服从柯西分布进行处理,单次进行,避免由于子波的混合相位因素造成反射剖面的极性发生变化,从而避免直接影响到反演波阻抗后速度变高还是变低,极大的增加反演结构的准确性。
发明应用实施例2:
使用本专利方法对洼陷带滨浅湖滩坝砂地区进行反演,洼陷带滨浅湖滩坝砂一般构造平缓,砂体多呈席状,存在的主要问题是砂泥岩薄互层的问题,般砂层和隔层厚度都小于10m,地震无法识别,同时在此地区的地震反演中对反射系数的标定,通常是根据井中反射系数来标定反褶积后的振幅值的,但是,求波阻抗是一个积分的过程,反褶积后的地震道振幅实际上还不是反射系数,而是相当于反射系数再褶积一个剩余子波,这个剩余子波一般在浅层主频高些,深层主频低些,频率低的波积分后数值偏大,会使深层产生偏大的波阻抗值,而一般采取的反演方式是采用稀疏脉冲反演,并根据砂体的连续性适当进行高频补偿提高分辨率可以取得一定的反演效果,在井资料允许的情况下,采用随机反演等纵向分辨率高的反演方法,但是依据前面的步骤(1)至(9)对洼陷带滨浅湖滩坝砂地区进行反演,通过常规的地震勘探方法多次采集地震数据,然后对于采集到的地震数据进行反复多次处理,通过公式计算出反射系数,然后将多次反射系数根据奇偶分解原理进行处理,同时除了考虑时变的振幅因素外,还要考虑时变的主频变化,将谱反演和波阻抗递推反演相结合,然后建立频率域反射系数反演目标函数,通过最小平方法得到各个采样点的反射系数,最后利用谱反演准确求取反射系数,从而能够更好的得到对该地区的反演结果。
尽管已经示出和描述了本发明的实施例,对于本领域的普通技术人员而言,可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型,本发明的范围由所附权利要求及其等同物限定。
Claims (1)
1.一种基于柯西分布的叠后波阻抗反演方法,步骤如下:
(1)采用常规的地震勘探方法采集地震数据,对地震数据进行常规处理之后得到叠后地震数据,然后对叠后地震数据进行进行处理得到层位数据;
(2)采用常规的测井方法得到测井数据,将测井数据进行处理得到测井声波时差曲线和密度曲线,然后通过分析测井声波时差曲线和密度曲线得到叠加速度谱质量曲线和分辨率曲线。
(3)根据叠后地震数据、层位数据以及测井数据实现测井分层到地震层位的标定和映射,由此得到时深关系曲线,然后根据测井曲线与时深关系曲线得到时间域曲线,接着根据时间域波抗曲线,生成初始波阻抗模型。
(4)把所述步骤(1)中的叠后地震数据、步骤(3)中生成的初始波阻抗模型一道数据和步骤(2)中得到的声波密度数据,代入以下公式,然后根据波的传播理论,当平面弹性波垂直入射到两种介质的平面分界面上时,通过反演得到该道的反射系数:
其中,ρ1和ρ2是界面两侧介质的密度;V1和V2是界面两侧介质的波速;Z=ρV是界面两侧介质的波阻抗;γ是界面反射系数;
当地面下的介质为层状介质,并存在着一系列互相平行的N个反射界面时,对于其中第i个界面,当波垂直入射时,其反射系数为:
当r、Z1已知时,可以由(1式)导出求Z2的公式
当r、Z2已知时,则:
(5)然后将所述步骤(4)中得到的反射系数数据,根据奇偶分解原理,将反射系数分解为奇偶分量,即:
r(t)=r0(t)+re(t) (5式)
其中,r0(t)=【r(t)-r(-t)】/2为奇分量;
re(t)=【r(t)+r(-t)】/2为偶分量;
然后建立频率域反射系数反演目标函数,通过最小平方法得到各个采样点的反射系数,最后利用谱反演准确求取反射系数。
(6)在所述步骤(5)中得到的数据基础上,当地下存在一系列互相平行的反射界面时,已知第N层界面以上N-1层界面的反射系数ri(i=1,2,...N-1)和第一层介质的波阻抗Z1,由(3式)可导出任意第n层介质的波阻抗Zn,得到:
(7)而已知第n层介质波阻抗Zn时,由(4式)又可导出任意第m层(m<n)介质的Zm,即:
(8)利用(6式)或(7式),可以从声波时差曲线及密度曲线上(没有密度测井时可利用Gardnar公式,ρ=0.31V0.25换算)选择标准层波阻抗作为基准波阻抗,将反褶积得到的反射系数转换为波阻抗。
(9)对所有的地震道重复所述步骤(6)至步骤(8)过程得到所有道的反演结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911290541.1A CN110954953A (zh) | 2019-12-16 | 2019-12-16 | 一种基于柯西分布的叠后波阻抗反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911290541.1A CN110954953A (zh) | 2019-12-16 | 2019-12-16 | 一种基于柯西分布的叠后波阻抗反演方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110954953A true CN110954953A (zh) | 2020-04-03 |
Family
ID=69981800
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911290541.1A Pending CN110954953A (zh) | 2019-12-16 | 2019-12-16 | 一种基于柯西分布的叠后波阻抗反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110954953A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112363221A (zh) * | 2020-10-27 | 2021-02-12 | 中国石油天然气集团有限公司 | Walkaway VSP测线的布设方法及装置 |
CN112485826A (zh) * | 2020-11-12 | 2021-03-12 | 中国地质大学(武汉) | 绝对波阻抗反演成像方法、装置、设备及存储介质 |
CN113311482A (zh) * | 2020-06-09 | 2021-08-27 | 中国海洋石油集团有限公司 | 一种基于叠前谱反演优化的高分辨率中深层储层预测方法 |
CN117310802A (zh) * | 2023-09-13 | 2023-12-29 | 成都捷科思石油天然气技术发展有限公司 | 一种深度域储层地震反演方法 |
WO2024002027A1 (zh) * | 2022-06-29 | 2024-01-04 | 深圳市中兴微电子技术有限公司 | 互连阻抗在线监控方法及*** |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104769458A (zh) * | 2014-07-15 | 2015-07-08 | 杨顺伟 | 一种基于柯西分布的叠后波阻抗反演方法 |
-
2019
- 2019-12-16 CN CN201911290541.1A patent/CN110954953A/zh active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104769458A (zh) * | 2014-07-15 | 2015-07-08 | 杨顺伟 | 一种基于柯西分布的叠后波阻抗反演方法 |
Non-Patent Citations (1)
Title |
---|
刘洁 等: "基于地震谱反演的地层速度估算方法及应用", 《石油地球物理勘探》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113311482A (zh) * | 2020-06-09 | 2021-08-27 | 中国海洋石油集团有限公司 | 一种基于叠前谱反演优化的高分辨率中深层储层预测方法 |
CN113311482B (zh) * | 2020-06-09 | 2022-12-23 | 中国海洋石油集团有限公司 | 一种基于叠前谱反演优化的高分辨率中深层储层预测方法 |
CN112363221A (zh) * | 2020-10-27 | 2021-02-12 | 中国石油天然气集团有限公司 | Walkaway VSP测线的布设方法及装置 |
CN112363221B (zh) * | 2020-10-27 | 2024-05-28 | 中国石油天然气集团有限公司 | Walkaway VSP测线的布设方法及装置 |
CN112485826A (zh) * | 2020-11-12 | 2021-03-12 | 中国地质大学(武汉) | 绝对波阻抗反演成像方法、装置、设备及存储介质 |
WO2024002027A1 (zh) * | 2022-06-29 | 2024-01-04 | 深圳市中兴微电子技术有限公司 | 互连阻抗在线监控方法及*** |
CN117310802A (zh) * | 2023-09-13 | 2023-12-29 | 成都捷科思石油天然气技术发展有限公司 | 一种深度域储层地震反演方法 |
CN117310802B (zh) * | 2023-09-13 | 2024-06-07 | 成都捷科思石油天然气技术发展有限公司 | 一种深度域储层地震反演方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110954953A (zh) | 一种基于柯西分布的叠后波阻抗反演方法 | |
CN108802812B (zh) | 一种井震融合的地层岩性反演方法 | |
WU et al. | S wave velocity structure beneath digital seismic stations of Yunnan province inferred from teleseismic receiver function modeling | |
WO2016008105A1 (zh) | 一种基于柯西分布的叠后波阻抗反演方法 | |
CN108897041B (zh) | 一种铀矿富集区的预测方法和装置 | |
CN109521474B (zh) | 一种三维双控下的叠前地质统计学反演方法 | |
CN104570067A (zh) | 一种地球物理勘探中相控地震反演方法 | |
CN105242318A (zh) | 一种确定砂体连通关系的方法及装置 | |
CN108957532B (zh) | 储层预测方法及装置 | |
CN113109875B (zh) | 一种全波形速度场约束下的盐下碳酸盐岩储层反演方法 | |
CN117406272A (zh) | 一种快速多元信息约束的反褶积宽频处理方法及装置 | |
Singh et al. | Neural networks and their applications in lithostratigraphic interpretation of seismic data for reservoir characterization | |
CN109188522B (zh) | 速度场构建方法及装置 | |
Wang et al. | Joint rough sets and Karhunen-Loève transform approach to seismic attribute selection for porosity prediction in a Chinese sandstone reservoir | |
CN109143335B (zh) | 一种合成地震记录的制作方法、***、介质及设备 | |
CN111830560B (zh) | 一种基于降秩算法的地震数据重建方法 | |
CN112213776B (zh) | 叠前道集和vsp资料联合层控q模型建立方法 | |
Butorin | Application of wavelet spectral decomposition for geological interpretation of seismic data | |
CN104375171B (zh) | 一种高分辨率地震反演方法 | |
Corchete | Shear-wave velocity structure of Antarctica from Rayleigh-wave analysis | |
Chakraborty et al. | Multi-trace Inversion with Selective Wedge to Increase Thin Layer Resolution in Horizontal Layer | |
Zhang* et al. | Compressed sensing Inversion for thin-bed resolution in depth domain | |
Lv et al. | Geological Modeling of Semideep Reservoir Based on Seismic Frequency Division Processing and Sedimentary Numerical Simulation | |
Dao et al. | Simultaneous source separation by shot collocation and strength variation | |
Fu et al. | USING THE WAVELET TRANSFORM FOR SEISMIC WAVE IMPEDANCE INVERSION |
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: 20200403 |
|
RJ01 | Rejection of invention patent application after publication |