CN106610505B - 一种基于dtw和aba联合的井震资料匹配方法 - Google Patents
一种基于dtw和aba联合的井震资料匹配方法 Download PDFInfo
- Publication number
- CN106610505B CN106610505B CN201611248347.3A CN201611248347A CN106610505B CN 106610505 B CN106610505 B CN 106610505B CN 201611248347 A CN201611248347 A CN 201611248347A CN 106610505 B CN106610505 B CN 106610505B
- Authority
- CN
- China
- Prior art keywords
- well
- data
- amplitude
- seismic
- sig
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 27
- 239000000463 material Substances 0.000 title claims abstract description 9
- 239000002131 composite material Substances 0.000 claims abstract description 31
- 238000012937 correction Methods 0.000 claims abstract description 11
- 238000000605 extraction Methods 0.000 claims abstract description 7
- 238000007788 roughening Methods 0.000 claims abstract description 4
- 238000005070 sampling Methods 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000011160 research Methods 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 5
- 238000012360 testing method Methods 0.000 claims description 5
- 238000012952 Resampling Methods 0.000 claims description 3
- 230000005856 abnormality Effects 0.000 claims description 3
- 230000003628 erosive effect Effects 0.000 claims description 3
- 238000007781 pre-processing Methods 0.000 claims description 3
- 238000007493 shaping process Methods 0.000 claims description 3
- 239000002002 slurry Substances 0.000 claims description 3
- 239000002689 soil Substances 0.000 claims description 3
- 239000004568 cement Substances 0.000 claims 1
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 238000013519 translation 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/282—Application of seismic models, synthetic seismograms
-
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/616—Data from specific type of measurement
- G01V2210/6169—Data from specific type of measurement using well-logging
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (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
本发明公开了一种基于DTW和ABA联合的井震资料匹配方法,其特征在于以下步骤:步骤S1:对测井数据进行粗化预处理,从地震资料中提取子波,利用反射系数与提取的地震子波褶积得到合成地震记录;步骤S2:采用动态时间规整算法(DTW)对测井合成记录序列和地震数据进行时间尺度上的匹配,通过拉伸和挤压达到地震数据和测井数据的时间尺度匹配;步骤S4:采用振幅均衡校正(ABA)来对地震道进行振幅匹配调整,通过振幅校正达到地震数据和测井数据的振幅大小上的匹配。本发明的有益效果是:考虑了井震数据在时间尺度和振幅大小上的特征,能够有效地完成井震匹配工作。
Description
技术领域
本发明涉及勘探地球物理领域,尤其涉及一种基于动态时间规整算法DTW和振幅均衡校正方法ABA联合的井震资料匹配方法。
背景技术
声波测井以及地面地震等多种较为先进的地球物理技术可以反映出地下相同目标的不同性质。但是由于测量仪器、测量环境以及处理技术等多种因素的不同,这些不同的地球物理技术各自得到的反映地下目标的信息会呈现出不同的特征。目前,由于油气勘探开发研究的深入,对不同地球物理技术所反映的相同信息的互相匹配,以及对声波测井和地面地震等多种信息的综合应用的研究日趋深入,因此目前研究的热点和难点,便是对地下目标进行更准确精细的描绘。井震不匹配特征具体表现如下几个方面:空间不同,地面地震资料是时间域信息,而声波测井资料却为深度域信息,对于地下相同目标,地面地震资料反映了整个三维空间信息,但声波测井资料则反映了某点的一维信息;分辨率不同,地面地震资料在纵向上分辨率较低,但在横向上却有很高的分辨能力,声波测井资料在横向上基本没有分别能力,但在纵向上却有较高的分辨率;数值不同,地面地震资料可以得到地层叠加剖面信息,而声波测井资料却可以得到地层波阻抗信息以及纵向声波速度。目前应用最广泛的井震匹配方法是闭合差法与比值法,两者均是将合成记录对应于井旁道进行整体的平移或进行局部的“拉伸”与“压缩”主要问题是存在一定的多解性,与人为干预有关很难得到最优解。刘景贤等(2012)在井、震匹配问题中引入了黏弹性岩石的杨氏谐振品质因子(Q)模型,对储层岩石进行了井震匹配与Q反演,通过谐振Q模型合理的校正声波速度,从而实现井震匹配。基于谐振模型的谐振Q理论的方法也有众多学者研究,但是由于地层品质因子等特征参数难以确定,因此其校正方法和实现的校正效果都受到很大的影响。
发明内容
本发明的目的就是针对现有技术不足,考虑测井合成资料和地震资料在时间尺度和振幅大小上存在差异,本发明提出了一种基于动态时间规整算法DTW和振幅均衡校正方法ABA联合的井震资料匹配方法,考虑了井震数据在时间尺度和振幅大小上的特征,能够有效地完成井震匹配工作。
一种基于DTW和ABA联合的井震资料匹配方法,包括如下步骤:
步骤S1:对测井数据进行预处理,针对测井数据的高采样率进行尺度粗化,从粗化后的测井数据中求取反射系数r,从井旁地震道中提取子波,利用反射系数r与理论子波或提取的地震子波褶积得到测井合成记录;
步骤S11:首先对声波测井获得的声波测试数据进行3-5次中值滤波处理,对声波测试数据序列采用滑动时窗进行初步平滑,剔除异常的数据和降低井径扩张和泥浆侵蚀的影响;然后通过Walsh变换将测井数据粗化,以便与地震数据匹配;
步骤S12:对资料做时频分析处理,考虑测井数据时频分析后的特征,确定测井数据中不同尺度的分层特征,采用最小速度差和最小厚度原理对声测井数据在深度域内进行精细分层处理,所允许的相邻两层的速度差是根据目标层位的研究需要来确定,速度差取50-200m/s;从实际地震剖面上确定时间采样间隔Δt,选取合理的分层参数ΔV,对ΔZ赋予初值Δh;计算分层后每一层的时间厚度,并找出最小的时间厚度t0min;判断t0min与Δt的大小,如果t0min等于或者第一次出现大于Δt的情况,则相应的ΔZ为最终的厚度分层参数;
步骤S13:将分层后的波阻抗P进行深时转换,并进行重采样,采样间隔与井旁道地震记录间隔一致,利用公式求出每个网格点上的反射系数,P1表示上层介质波阻抗,P2表示下层介质波阻抗;利用反射系数r与理论子波或提取的地震子波褶积得到合成地震记录sig1(t);
步骤S2:采用动态时间规整算法DTW对测井合成记录sig1(t)和井旁地震道sig2(t)进行时间尺度上的匹配,通过拉伸和挤压达到地震数据和测井数据的时间尺度匹配;两个序列的长度分别为n和m,则构建一个n×m的矩阵网格,矩阵元素d(i,j)表示和两个点的距离寻找一条通过此网格中若干格点的路径,路径通过的格点即为两个序列进行计算的对齐的点;
步骤S3:采用振幅均衡校正ABA来对测井合成记录进行振幅匹配调整,考虑到使信号的振幅或功率变化保持一致,将信号分为一个个重叠的等宽度时窗,在每个时窗内计算两个信号的平均绝对值振幅:
sig(n)表示测井合成记录与井旁地震道,n表示离散点序号,g(n)表示滑动窗函数,T表示时窗长度。
时窗内信号振幅平均值为:
振幅均衡校正过程通过把测井合成记录sig1(t)的振幅朝着井旁地震道sig2(t)的振幅的方向平衡来得到标量函数:
最后用测井合成记录sig1(t)与标量函数相乘即可完成振幅均衡校正。
附图说明
附图1为本发明的主要方法步骤。
附图2为测井数据预处理结果。
附图3为测井合成记录(未调整)和井旁道地震记录。
附图4为DTW校正路径。
附图5为DTW后测井合成记录和井旁道地震记录。
附图6为实际数据井震匹配结果,左边和右边的三道为井旁地震道,中间为DTW和ABA联合校正后的合成记录道。
具体实施方式
结合附图1-6,对本发明作进一步的描述:一种基于DTW和ABA联合的井震资料匹配方法,具体实施方案包括如下步骤:
步骤S1:对测井数据进行预处理,针对测井数据的高采样率进行尺度粗化,从地震资料中提取子波,利用反射系数r与理论子波或提取的地震子波褶积得到测井合成记录;
步骤S11:首先对声波测井获得的声波测试数据进行3-5次中值滤波处理,对声波测试数据序列采用滑动时窗进行初步平滑,剔除异常的数据和降低井径扩张和泥浆侵蚀的影响;然后通过Walsh变换将测井数据粗化,以便与地震数据匹配;
步骤S12:对资料做时频分析处理,考虑测井数据时频分析后的特征,确定测井数据中不同尺度的分层特征,采用最小速度差和最小厚度原理对声测井数据在深度域内进行精细分层处理后,所允许的相邻两层的速度差通常是根据目标层位的研究需要来确定,所允许的相邻两层的速度差是根据目标层位的研究需要来确定,速度差取50-200m/s;从实际地震剖面上确定时间采样间隔Δt,选取合理的分层参数ΔV,对ΔZ赋予初值Δh;计算分层后每一层的时间厚度,并找出最小的时间厚度t0min;判断t0min与Δt的大小,如果t0min等于或者第一次出现大于Δt的情况,则相应的ΔZ为最终的厚度分层参数;图2为测井数据预处理结果,由图可知经过处理后测井数据得到了方块化(粗化);
步骤S13:将分层后的波阻抗P进行深时转换,并进行重采样,采样间隔与井旁道地震记录间隔一致,利用公式求出每个网格点上的反射系数。利用反射系数r与理论子波或提取的地震子波褶积得到合成地震记录sig1(t);图3为测井合成记录(未调整)和井旁道地震记录,可以看出未调整的测井合成记录和井旁地震道在时间尺度和振幅大小上都存在较大差异;
步骤S2:采用动态时间规整算法DTW对测井合成记录序列sig1(t)和sig2(t)地震数据进行时间尺度上的匹配,通过拉伸和挤压达到地震数据和测井数据的时间尺度匹配;两个序列的长度分别为n和m,则构建一个n×m的矩阵网格,矩阵元素d(i,j)表示和两个点的距离寻找一条通过此网格中若干格点的路径,路径通过的格点即为两个序列进行计算的对齐的点;图4为DTW校正路径,图5为DTW后测井合成记录和井旁地震道,从图可知通过DTW处理后测井合成记录和井旁地震道在时间尺度上能够较好的对应起来。
步骤S3:采用振幅均衡校正ABA来对测井合成记录进行振幅匹配调整,考虑到使信号的振幅或功率变化保持一致,将信号分为一个个重叠的等宽度时窗,在每个时窗内计算两个信号的平均绝对值振幅:
时窗内信号振幅平均值为:
振幅均衡校正过程通过把测井合成记录sig1(t)的振幅朝着井旁地震道sig2(t)的振幅的方向平衡来得到标量函数:
最后用合成记录sig1(t)与标量函数相乘即可完成振幅均衡校正。图6所示为合成记录与井旁道记录同相轴对应准确,振幅差异较小。
通过以上流程,能够有效完成井震匹配工作。图2为测井数据预处理结果,图3为未经调整的测井合成记录和井旁道地震记录,图4为DTW校正路径,图5为DTW后测井合成记录和井旁道地震记录,从图5中可以看出处理后的测井合成结果和井旁地震道在时间尺度上能够很好地对应起来;图6是实际数据井震匹配结果,从图中可以看出通过联合处理后测井合成结果和井旁地震道的幅值也能够较好的匹配起来。应用基于DTW和ABA联合的井震资料匹配方法,考虑了井震数据在时间尺度和幅值大小上的特征,能够有效地完成井震匹配工作。
本发明的有益效果是:本发明针对井震资料在资料类型、分辨率等方面的差异,考虑了井震数据在时间尺度和振幅大小上的特征,克服了基于谐振模型需要难以确定的地层品质因子参数的缺点以及手动调整因人而异的缺点,选用了动态时间规整方法和振幅均衡校正方法能够自动有效地完成井震匹配工作,匹配后的数据在时间尺度和幅值大小能够很好地匹配起来,为后续资料的处理和解释提供强有力的技术支撑。
本发明不局限于上述具体实施方式,根据上述内容,按照本领域的普通技术知识和惯用手段,在不脱离本发明上述基本技术思想前提下,本发明还可以做出其它多种形式的等效修改、替换或变更,均属于本发明的保护范围。
Claims (2)
1.一种基于DTW和ABA联合的井震资料匹配方法,其特征在于包括如下步骤:
步骤S1:对测井数据进行预处理,针对测井数据的高采样率进行尺度粗化,从粗化后的测井数据中求取反射系数r,从井旁地震道中提取子波,利用反射系数r与理论子波或提取的地震子波褶积得到测井合成记录;
步骤S2:采用动态时间规整算法DTW对测井合成记录sig1(t)和井旁地震道sig2(t)进行时间尺度上的匹配,通过拉伸和挤压达到地震数据和测井数据的时间尺度匹配;两个序列的长度分别为n和m,则构建一个n×m的矩阵网格,矩阵元素d(i,j)表示和两个点的距离寻找一条通过此网格中若干格点的路径,路径通过的格点即为两个序列进行计算的对齐的点;
步骤S3:采用振幅均衡校正ABA来对测井合成记录和井旁地震道进行振幅匹配调整,考虑到使信号的振幅或功率变化保持一致,将信号分为一个个重叠的等宽度时窗,在每个时窗内计算两个信号的平均绝对值振幅:
时窗内信号振幅平均值为:
sig(n)表示测井合成记录与井旁地震道,n表示离散点序号,g(n)表示滑动窗函数,T表示时窗长度;
振幅均衡校正过程通过把测井合成记录sig1(t)的振幅朝着井旁地震道sig2(t)的振幅的方向平衡来得到标量函数:
sig1(n)表示测井合成记录,sig2(n)表示井旁地震道,n表示离散点序号,g(n)表示滑动窗函数,T表示时窗长度,τ表示窗函数中心点;
最后用测井合成记录sig1(t)与标量函数相乘即可完成振幅均衡校正。
2.根据权利要求1所述的一种基于DTW和ABA联合的井震资料匹配方法,其特征在于步骤S1包括如下步骤:
步骤S11:首先对声波测井获得的声波测试数据进行3-5次中值滤波处理,对声波测井数据采用滑动时窗进行初步平滑,剔除异常的数据和降低井径扩张和泥浆侵蚀的影响;然后通过Walsh变换将测井数据粗化,以便与地震数据匹配;
步骤S12:对资料做时频分析处理,考虑声波测井数据时频分析后的特征,确定声波测井数据中不同尺度的分层特征,采用最小速度差和最小厚度原理对声波测井数据在深度域内进行精细分层处理,所允许的相邻两层的速度差是根据目标层位的研究需要来确定,速度差取50-200m/s;从实际地震剖面上确定时间采样间隔Δt,选取合理的分层参数ΔV,对ΔZ赋予初值Δh;计算分层后每一层的时间厚度,并找出最小的时间厚度t0min;判断t0min与Δt的大小,如果t0min等于或者第一次出现大于Δt的情况,则相应的ΔZ为最终的厚度分层参数;
步骤S13:将分层后的波阻抗P进行深时转换,并进行重采样,采样间隔与井旁地震道间隔一致,利用公式求出每个网格点上的反射系数,P1表示上层介质波阻抗,P2表示下层介质波阻抗;利用反射系数r与理论子波或提取的地震子波褶积得到测井合成记录sig1(t)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611248347.3A CN106610505B (zh) | 2016-12-29 | 2016-12-29 | 一种基于dtw和aba联合的井震资料匹配方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611248347.3A CN106610505B (zh) | 2016-12-29 | 2016-12-29 | 一种基于dtw和aba联合的井震资料匹配方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106610505A CN106610505A (zh) | 2017-05-03 |
CN106610505B true CN106610505B (zh) | 2019-03-22 |
Family
ID=58636507
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611248347.3A Expired - Fee Related CN106610505B (zh) | 2016-12-29 | 2016-12-29 | 一种基于dtw和aba联合的井震资料匹配方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106610505B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20210132253A1 (en) * | 2019-11-01 | 2021-05-06 | Saudi Arabian Oil Company | Automatic geological formations tops picking using dynamic time warping (dtw) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107831542B (zh) * | 2017-07-24 | 2019-05-31 | 中国石油化工股份有限公司 | Ddw高精度深度域井震匹配方法 |
CN107329170B (zh) * | 2017-08-07 | 2019-02-15 | 中国石油天然气集团公司 | 一种地震波形匹配方法及装置 |
CN107563297B (zh) * | 2017-08-07 | 2020-06-09 | 中国石油天然气集团公司 | 一种波形匹配方法及装置 |
CN107356970B (zh) * | 2017-09-11 | 2020-05-01 | 中国石油大学(华东) | 一种高精度的井震数据匹配方法 |
CN107843919B (zh) * | 2017-10-27 | 2020-02-07 | 中国石油天然气集团有限公司 | 一种微地震监测数据振幅均衡方法 |
CN112824940B (zh) * | 2019-11-20 | 2024-06-25 | 中国石油天然气集团有限公司 | 反射纵波与反射转换横波的振幅匹配方法及装置 |
CN117333575B (zh) * | 2023-06-21 | 2024-05-31 | 广州海洋地质调查局 | 一种井筒分布式光纤声振动监测数据成像方法及处理终端 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102062873A (zh) * | 2009-11-13 | 2011-05-18 | 中国石油化工股份有限公司 | 一种纵横波匹配方法 |
CN103675904A (zh) * | 2012-09-10 | 2014-03-26 | 中国石油集团东方地球物理勘探有限责任公司 | 一种井震匹配解释性目标处理方法及装置 |
CN104237946A (zh) * | 2014-09-19 | 2014-12-24 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 基于井控的单层反射纵波和反射转换横波的振幅匹配方法 |
CN104237938A (zh) * | 2014-09-04 | 2014-12-24 | 中国石油天然气集团公司 | 利用动态时间调整进行纵波和转换波数据匹配的方法 |
CN104536042A (zh) * | 2014-12-24 | 2015-04-22 | 中国石油天然气股份有限公司 | 一种二维叠后地震资料振幅补偿方法及装置 |
CN105044773A (zh) * | 2015-08-24 | 2015-11-11 | 中国石油集团东方地球物理勘探有限责任公司 | 转换波井震标定方法及装置 |
CN105445802A (zh) * | 2014-09-04 | 2016-03-30 | 中国石油化工股份有限公司 | 一种基于时频谱约束的高精度时深校正方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060256657A1 (en) * | 2005-05-11 | 2006-11-16 | Prism Seismic, Inc. | Method for improving the time-depth tie of well log data and seismic data |
WO2016115437A1 (en) * | 2015-01-15 | 2016-07-21 | Schlumberger Canada Limited | Approaches to tying well logs to seismic data |
-
2016
- 2016-12-29 CN CN201611248347.3A patent/CN106610505B/zh not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102062873A (zh) * | 2009-11-13 | 2011-05-18 | 中国石油化工股份有限公司 | 一种纵横波匹配方法 |
CN103675904A (zh) * | 2012-09-10 | 2014-03-26 | 中国石油集团东方地球物理勘探有限责任公司 | 一种井震匹配解释性目标处理方法及装置 |
CN104237938A (zh) * | 2014-09-04 | 2014-12-24 | 中国石油天然气集团公司 | 利用动态时间调整进行纵波和转换波数据匹配的方法 |
CN105445802A (zh) * | 2014-09-04 | 2016-03-30 | 中国石油化工股份有限公司 | 一种基于时频谱约束的高精度时深校正方法 |
CN104237946A (zh) * | 2014-09-19 | 2014-12-24 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 基于井控的单层反射纵波和反射转换横波的振幅匹配方法 |
CN104536042A (zh) * | 2014-12-24 | 2015-04-22 | 中国石油天然气股份有限公司 | 一种二维叠后地震资料振幅补偿方法及装置 |
CN105044773A (zh) * | 2015-08-24 | 2015-11-11 | 中国石油集团东方地球物理勘探有限责任公司 | 转换波井震标定方法及装置 |
Non-Patent Citations (3)
Title |
---|
A semiautomatic method to tie well logs to seismic data;Roberto H.Herrera 等;《GEOPHYSICS》;20140630;第79卷(第3期);第V47-V54页 |
Automatic approaches for seismic to well tying;Roberto H.Herrera1 等;《Interpretation》;20140531;第2卷(第2期);第SD101-SD109页 |
井震资料尺度匹配过程中声测井数据的精细分层方法研究;周建科 等;《物探化探计算技术》;20150331;第37卷(第2期);第242-246页 |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20210132253A1 (en) * | 2019-11-01 | 2021-05-06 | Saudi Arabian Oil Company | Automatic geological formations tops picking using dynamic time warping (dtw) |
Also Published As
Publication number | Publication date |
---|---|
CN106610505A (zh) | 2017-05-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106610505B (zh) | 一种基于dtw和aba联合的井震资料匹配方法 | |
CN106855636B (zh) | 基于碳酸盐岩储层露头的原型地质模型地震正演方法 | |
CN106873033B (zh) | 一种基于神经网络的多属性地震信息融合裂缝预测的方法 | |
CN104614765B (zh) | 一种增强地震波激发照明的设计方法 | |
CN106443768B (zh) | 一种叠前深度域合成地震记录制作方法 | |
CN103954995B (zh) | 一种砂岩型铀矿勘查中的砂体识别方法 | |
CN104820243B (zh) | 一种三维叠后地震资料振幅补偿方法及装置 | |
CN109343120B (zh) | 融入约束稀疏脉冲反演低频补偿的声波曲线重构方法 | |
CN106597545B (zh) | 一种水平裂缝地震叠前反演方法和装置 | |
CN107065013B (zh) | 一种地震尺度下的层速度确定方法及装置 | |
CN110954943A (zh) | 被动源地震频率谐振勘探方法 | |
CN108398720A (zh) | 一种基于杨氏模量、泊松比的两项式地震叠前反演方法 | |
CN110007340A (zh) | 基于角度域直接包络反演的盐丘速度密度估计方法 | |
CN109425900A (zh) | 一种地震储层预测方法 | |
CN110261898A (zh) | 基于地震岩石物理实验分析的测井与地震速度匹配方法 | |
CN106842299B (zh) | 一种基于地震属性的裂缝定量化预测的方法 | |
CN108594301B (zh) | 一种具有差异特征的地震数据融合的方法及处理终端 | |
Holliger et al. | Seismic structure of gneissic/granitic upper crust: geological and petrophysical evidence from the Strona-Ceneri Zone (northern Italy) and implications for crustal seismic exploration | |
Cai et al. | Intelligent calculation method of relative sonic attenuation and its application to fracture evaluation in tight sandstone reservoir | |
CN102914797A (zh) | 一种获得地层各向异性系数的方法及装置 | |
CN109324344A (zh) | 基于纯纵波和拟声波反演的页岩厚度预测方法及系统 | |
CN110244383A (zh) | 基于近地表数据的地质岩性综合模型创建方法 | |
CN104570104B (zh) | 一种基于两步法avf的纵横波地震品质因子提取方法 | |
CN107561581B (zh) | 基于相关系数法建立多井模型的方法 | |
Pratt et al. | Crosswell waveform tomography in fine-layered sediments-Meeting the challenges of anisotropy |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190322 Termination date: 20191229 |