CN111598340B - 基于分数阶希尔伯特变换的薄砂体平面展布预测方法 - Google Patents
基于分数阶希尔伯特变换的薄砂体平面展布预测方法 Download PDFInfo
- Publication number
- CN111598340B CN111598340B CN202010419191.0A CN202010419191A CN111598340B CN 111598340 B CN111598340 B CN 111598340B CN 202010419191 A CN202010419191 A CN 202010419191A CN 111598340 B CN111598340 B CN 111598340B
- Authority
- CN
- China
- Prior art keywords
- sand body
- hilbert transform
- plane
- development
- interval
- 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.)
- Active
Links
- 239000004576 sand Substances 0.000 title claims abstract description 137
- 238000000034 method Methods 0.000 title claims abstract description 44
- 238000011161 development Methods 0.000 claims abstract description 53
- 238000005553 drilling Methods 0.000 claims abstract description 15
- 238000010586 diagram Methods 0.000 claims abstract description 6
- 238000004458 analytical method Methods 0.000 claims description 17
- 238000005070 sampling Methods 0.000 claims description 12
- 238000003491 array Methods 0.000 claims description 10
- 238000011160 research Methods 0.000 claims description 5
- GNFTZDOKVXKIBK-UHFFFAOYSA-N 3-(2-methoxyethoxy)benzohydrazide Chemical compound COCCOC1=CC=CC(C(=O)NN)=C1 GNFTZDOKVXKIBK-UHFFFAOYSA-N 0.000 claims description 2
- 230000008569 process Effects 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 description 8
- 238000012545 processing Methods 0.000 description 6
- 238000007796 conventional method Methods 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 3
- 230000010363 phase shift Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 230000008021 deposition Effects 0.000 description 2
- 238000003708 edge detection Methods 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- 238000004062 sedimentation Methods 0.000 description 2
- FGUUSXIOTUKUDN-IBGZPJMESA-N C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 Chemical compound C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 FGUUSXIOTUKUDN-IBGZPJMESA-N 0.000 description 1
- 208000035126 Facies Diseases 0.000 description 1
- 241000446313 Lamella Species 0.000 description 1
- 101100004031 Mus musculus Aven gene Proteins 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- 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/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/40—Transforming data representation
- G01V2210/48—Other transforms
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Strategic Management (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mathematical Physics (AREA)
- Remote Sensing (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Marine Sciences & Fisheries (AREA)
- General Health & Medical Sciences (AREA)
- Animal Husbandry (AREA)
- Agronomy & Crop Science (AREA)
- Primary Health Care (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Development Economics (AREA)
- Mining & Mineral Resources (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Algebra (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
Abstract
本发明公开了一种基于分数阶希尔伯特变换的薄砂体平面展布预测方法，包括：先计算10个阶数的三维地震资料的分数阶希尔伯特变换，再对砂体发育段在纵向上进行多个砂体发育小层段的划分，然后提取各个砂体发育小层段的属性平面图，基于实钻井在各个砂体发育小层段的发育情况优选最佳的阶数对应的振幅属性平面图进行薄砂体平面展布预测。本发明充分利用了分数阶希尔伯特变换能提高纵向分辨率和能在地震剖面上实现同相轴和声波测井曲线完美对应的优点，有效地预测薄砂体的平面展布特征，提高薄砂体平面展布预测的准确性。
Description
技术领域
本发明涉及油气田勘探技术领域，尤其涉及一种基于分数阶希尔伯特变换的薄砂体平面展布预测方法。
背景技术
薄砂体预测的重要性：常规的石油天然气勘探的储层主要是砂体，近年来随着勘探程度的不断深入，厚砂体的油层开采难度越来越大，水淹程度越来越高；薄砂体成了目前油气勘探的重点目标，其主要发育在河流—三角洲沉积末端，以三角洲分流平原末端和内前缘为主，骨架砂体以小型分流河道砂体和水下河道沉积为主。目前，学术界通常将厚度小于5米的砂体定义为薄砂体。
常规的薄砂体预测技术：目前，对于薄砂体的预测技术非常多，主要集中在地震反演技术和地震属性分析技术两大领域：一是地震反演技术，主要利用泥岩具有高波阻抗、砂岩具有低波阻抗的特征，进行波阻抗反演(参考文献：郝杰、吴鑫、孙明灯，南堡地区浅层河道砂体的识别，石油地球物理勘探，2018，53卷(增刊1)：151-157)；二是地震属性分析技术，其细分为(1)寻找砂体边界：如相干体技术、曲率分析技术、边缘检测技术等，(2)寻找砂体(岩性识别)：如瞬时吸收系数Q属性分析技术、波峰振幅属性分析技术等(参考文献：王彬、朱帅润、王量等，基于地震属性分析预测河道砂体，油气地球物理，2019，17(2)：45-49)。
目前，希尔伯特变换主要应用于4类领域：(1)基于复数道分析，提取瞬时属性(如瞬时相位、瞬时振幅、瞬时频率等)用于薄砂体预测(岩性识别)(参考文献：贺振华，地震波复数道分析的计算方法，1980，第4期，64-80)；(2)基于地震边缘检测分析，用于薄砂体的平面展布预测，如一维时间域或频率域Hilbert变换的方法、二维Hilbert变换的方法及广义希尔伯特变换(GHT)的方法等(参考文献：熊晓军、贺振华、赵明金等，一种基于GHT的裂缝检测新方法，石油地球物理勘探，2009，44(4)，442-444)；(3)对三维地震资料进行900相移处理，进行砂岩的地震沉积相带分析(参考文献：朱筱敏、曾洪流、董艳蕾，地震沉积学(译文集)，北京：石油工业出版社，2011)；(4)针对测井资料进行高分辨率处理，用于砂岩的岩性识别(参考文献：刘洪林、周磊，希尔伯特变换在测井资料数据处理中的应用，1998，22(3)：9-11)。
常规希尔伯特变换方法的缺点：利用希尔伯特变换相移900的方法，可以用极性来指示岩性，采用900相位子波处理的地震数据就可以克服零相位子波数据的不足。地震响应的最大振幅经过900移动与薄层中心相对应时，地震同相轴与地质上限定的砂岩层相一致，使解释工作较以前变得容易。但是当子波不是零相位，或当地层沉积模式是正常沉积时，再用相移900的办法已不能使地震极性和振幅很好地指示岩性。
目前，分数阶希尔伯特变换主要应用于通信信号分析领域。在油气勘探领域公开发表的文献仅有1篇(宋洪勇、别立珍，分数阶希尔伯特变换在提高纵向分辨率中的应用，2013)，其应用聚焦在提高纵向分辨率方面。尚没有基于分数阶希尔伯特变换进行薄砂体预测的研究。
发明内容
本发明的目的在于克服现有技术存在的上述问题，提供一种基于分数阶希尔伯特变换的薄砂体平面展布预测方法。本发明充分利用了分数阶希尔伯特变换能提高纵向分辨率和能在地震剖面上实现同相轴和声波测井曲线完美对应的优点，有效地预测薄砂体的平面展布特征，提高薄砂体平面展布预测的准确性。
为实现上述目的，本发明采用的技术方案如下：
一种基于分数阶希尔伯特变换的薄砂体平面展布预测方法，其特征在于：先计算10个阶数的三维地震资料的分数阶希尔伯特变换，再对砂体发育段在纵向上进行多个砂体发育小层段的划分，然后提取各个砂体发育小层段的属性平面图，基于实钻井在各个砂体发育小层段的发育情况优选最佳的阶数对应的振幅属性平面图进行薄砂体平面展布预测。
所述方法具体包括如下步骤：
a、输入叠后三维地震数据体；
b、分别计算10个阶数值的分数阶希尔伯特变换；
c、根据井震标定分析得出砂体段发育的纵向位置，对10个分数阶希尔伯特变换的结果数组在纵向上提取各个砂体发育小层段的10个三维层段数组；
d、统计研究区内的实钻井在各个砂体发育小层段内的砂体发育情况；
e、对纵向上发育的第一个砂体发育小层段，计算当前层段内的三维数组的10个均方根振幅属性，并绘制10个均方根振幅平面图；
f、根据实钻井在当前层段的砂体的发育情况，计算步骤e中10个均方根振幅平面图的砂体平面展布预测符合率，选取预测符合率最高的阶数对应的均方根振幅平面图进行薄砂体平面展布预测；
g、对纵向上的其它砂体发育的小层段，重复步骤e和f进行薄砂体平面展布预测。
所述步骤a中，输入叠后三维地震数据体w(i,j,n)，i代表Inline测线序号，总测线数记为Nline；j代表Xline测线序号，总测线数记为Nxline；n代表单道地震记录的采样点序号，总采样点数记为N。
所述步骤b中，令v＝-1.0,-0.8,-0.6,-0.4,-0.2,0.2,0.4,0.6,0.8,1，对三维数组w(i,j,n)的每个地震道记录进行分数阶希尔伯特变换，得到各个阶数的三维数组w(i,j,n)的分数阶希尔伯特变换的结果数组wv(i,j,n)，共10个三维数组。
所述步骤b具体过程如下：
b1、令v＝-1.0；
b2、令i＝1,j＝1提取当前地震道数组s(n)，s(n)为一维数组，数组的长度是N；
b3、对目标一维数组进行s(n)做DFT，得到S(k)(k＝0,1,…,N-1)，其中k＝N/2,…,N-1对应负频率；
b4、对S(k)进行分频率计算：
b5、对y(k)进行逆DFT，取其实部即得到s(n)的分数阶希尔伯特变换sv(n)，并写入三维数组wv(i,j,n)；
b6、增加i,j的数值，重复步骤b2和步骤b5，直到i＝Nline,j＝Nxline；得到当前阶数v的分数阶希尔伯特变换的三维数组wv(i,j,n)；
b7、改变v的数值，重复步骤b2和步骤b6，直到v＝1.0，得到10个阶数对应的10个三维数组。
所述步骤c中，根据井震标定分析得出砂体段发育的纵向位置，对步骤b的分数阶希尔伯特变换的结果数组wv(i,j,n)在纵向上采用等时切片方法提取各个砂体发育小层段的三维层段数组Sandl_v(i,j,nc)，l代表纵向上砂体发育小层段的序号，纵向上砂体发育小层段总数记为Nsand；nc代表当前层段内的地震记录的采样点序号，总采样点数记为Nc。
所述步骤d中，令研究区内的实钻井数目是Nwell，统计步骤c中各个层段内的单井砂体发育情况，并写入二维数组数组P(a,l)，a代表实钻井序号，l代表纵向上发育砂体层段的序号，P(a,l)＝0代表该井在该层段砂体不发育，P(a,l)＝1代表该井在该层段砂体发育。
所述步骤e中，令l＝1，计算当前砂体发育小层段内的三维数组Sandl_v(i,j,nc)(v＝-1.0,-0.8,-0.6,-0.4,-0.2,0.2,0.4,0.6,0.8,1)的均方根振幅属性，记为Amp_wv(i,j)(v＝-1.0,-0.8,-0.6,-0.4,-0.2,0.2,0.4,0.6,0.8,1)，并绘制均方根振幅平面图。
所述步骤f中，根据实钻井在当前层段l的砂体的发育情况P(a,l)，计算步骤e中10个均方根振幅平面图的砂体平面展布预测符合率，选取预测符合率最高的阶数对应的均方根振幅平面图进行薄砂体平面展布预测。
所述步骤g中，增加l的数值，重复步骤e和步骤f，直到l＝Nsand。
采用本发明的优点在于：
(1)本发明隶属于油气田勘探技术领域，但是发明人从通信信号处理分析和图像处理领域的知识，引入“分数阶希尔伯特变换能提高纵向分辨率，且能在地震剖面上实现同相轴和声波测井曲线正好对应”的优点，对叠后三维地震数据进行分数阶希尔伯特变换，再提取振幅属性进行薄砂体的平面展布预测，较常规方法直接针对叠后三维地震数据提取振幅属性进行薄砂体的平面展布预测，具有较强的创新性。
(2)为了克服阶数不能准确设定的难点，根据分数阶希尔伯特变换的周期性特征，设计10个阶数值，并根据实钻井的砂体发育情况进行阶数优选，获得针对每个砂体发育层段的最佳阶数。
(3)考虑到纵向上往往发育多套薄砂体的特征，引入等时切片方法，根据井震标定的结果，进行纵向上多个砂体发育小层段的划分，避免了砂体纵向叠置带来的预测误差，提高了砂体平面展布预测的精度。
(4)本发明通过对三维地震资料进行分数阶希尔伯特变换处理，提高地震资料的纵向分辨率，并实现薄砂体的反射同相轴与测井的砂体位置完美对应，再提取小层段内的振幅属性进行薄砂体的平面展布预测，可以有效地预测薄砂体的平面展布特征，提高薄砂体平面展布预测的准确性。
附图说明
图1为采用常规方法(基于希尔伯特变换的薄层识别方法)计算得到的过W1井的砂体预测剖面图；
图2是采用本发明方法计算得到的过W1井的砂体预测剖面图。
具体实施方式
实施例1
本发明实施例提供的一种基于分数阶希尔伯特变换的薄砂体平面展布预测方法，下面结合实施例对本发明进行详细的描述。
本发明具体步骤如下：
(1)输入叠后三维地震数据体w(i,j,n)，i代表Inline测线序号(总测线数记为Nline)，j代表Xline测线序号(总测线数记为Nxline)，n代表单道地震记录的采样点序号(总采样点数记为N)；
(2)令v＝-1.0,-0.8,-0.6,-0.4,-0.2,0.2,0.4,0.6,0.8,1，对三维数组w(i,j,n)的每个地震道记录进行分数阶希尔伯特变换，得到各个阶数的三维数组w(i,j,n)的分数阶希尔伯特变换的结果数组wv(i,j,n)，共10个三维数组，此步骤具体实施步骤如下：
(a)令v＝-1.0；
(b)令i＝1,j＝1提取当前地震道数组s(n)(一维数组，数组的长度是N)；
(c)对目标一维数组进行s(n)做DFT，得到S(k)(k＝0,1,…,N-1)，其中k＝N/2,…,N-1对应负频率；
(d)对S(k)进行分频率计算，
(e)对y(k)进行逆DFT，取其实部即得到s(n)的分数阶希尔伯特变换sv(n)，并写入三维数组wv(i,j,n)；
(f)增加i,j的数值，重复步骤b和步骤e，直到i＝Nline,j＝Nxline；得到当前阶数v的分数阶希尔伯特变换的三维数组wv(i,j,n)；
(g)改变v的数值，重复步骤b和步骤f，直到v＝1.0，得到10个阶数对应的10个三维数组；
(3)根据井震标定分析得出砂体段发育的纵向位置，对步骤2的分数阶希尔伯特变换的结果数组wv(i,j,n)在纵向上采用等时切片方法提取各个砂体发育小层段的三维层段数组Sandl_v(i,j,nc)，l代表纵向上砂体发育小层段的序号(纵向上砂体发育小层段总数记为Nsand)，nc代表当前层段内的地震记录的采样点序号(总采样点数记为Nc)；
(4)令研究区内的实钻井数目是Nwell，统计步骤3中各个层段内的单井砂体发育情况，并写入二维数组数组P(a,l)，a代表实钻井序号，l代表纵向上发育砂体层段的序号，P(a,l)＝0代表该井在该层段砂体不发育，P(a,l)＝1代表该井在该层段砂体发育；
(5)令l＝1，计算当前砂体发育小层段内的三维数组Sandl_v(i,j,nc)(v＝-1.0,-0.8,-0.6,-0.4,-0.2,0.2,0.4,0.6,0.8,1)的均方根振幅属性(共10个)，记为Amp_wv(i,j)(v＝-1.0,-0.8,-0.6,-0.4,-0.2,0.2,0.4,0.6,0.8,1)，并绘制均方根振幅平面图(共10个)；
(6)根据实钻井在当前层段(l)的砂体的发育情况P(a,l)，计算步骤5中10个均方根振幅平面图的砂体平面展布预测符合率，选取预测符合率最高的阶数对应的均方根振幅平面图进行薄砂体平面展布预测；
(7)增加l的数值，重复步骤5和步骤6，直到l＝Nsand。
实施例2
图1是采用常规方法(基于希尔伯特变换的薄层识别方法)计算得到的过W1井的砂体预测剖面，图中的黑色薄层代表砂体，W1井旁的曲线是实测的GR曲线(GR曲线低值指示砂体，GR曲线高值指示泥岩地层)。从图1可见常规方法2项缺陷：①地层顶部出现一套黑色的连续性较好的地层(地点位置：横坐标470道，纵坐标2098ms)，而该套地层对应的GR值是高值，即该套地层不是砂体，出现了预测假象；②在地层倾角变化时，图中横坐标510道位置，黑色地层断续，即横向连续性差，横向分辨率较低。
图2是采用本发明方法计算得到的过W1井的砂体预测剖面(与图1对应)，图中的黑色薄层代表砂体，W1井旁的曲线是实测的GR曲线(GR曲线低值指示砂体，GR曲线高值指示泥岩地层)。从图2可见本发明方法预测得到的黑色砂体不仅与井旁GR曲线特征一致，具有较高的纵向分辨率(图2中的2套黑色地层纵向上完全分开，分辨率优于图1)，而且横向连续性好，可连续追踪。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
Claims (9)
1.一种基于分数阶希尔伯特变换的薄砂体平面展布预测方法，其特征在于：先计算10个阶数的三维地震资料的分数阶希尔伯特变换，再对砂体发育段在纵向上进行多个砂体发育小层段的划分，然后提取各个砂体发育小层段的属性平面图，基于实钻井在各个砂体发育小层段的发育情况优选最佳的阶数对应的振幅属性平面图进行薄砂体平面展布预测；
所述方法具体包括如下步骤：
a、输入叠后三维地震数据体；
b、分别计算10个阶数值的分数阶希尔伯特变换；
c、根据井震标定分析得出砂体段发育的纵向位置，对10个分数阶希尔伯特变换的结果数组在纵向上提取各个砂体发育小层段的10个三维层段数组；
d、统计研究区内的实钻井在各个砂体发育小层段内的砂体发育情况；
e、对纵向上发育的第一个砂体发育小层段，计算当前层段内的三维数组的10个均方根振幅属性，并绘制10个均方根振幅平面图；
f、根据实钻井在当前层段的砂体的发育情况，计算步骤e中10个均方根振幅平面图的砂体平面展布预测符合率，选取预测符合率最高的阶数对应的均方根振幅平面图进行薄砂体平面展布预测；
g、对纵向上的其它砂体发育的小层段，重复步骤e和f进行薄砂体平面展布预测。
2.根据权利要求1所述的基于分数阶希尔伯特变换的薄砂体平面展布预测方法，其特征在于：所述步骤a中，输入叠后三维地震数据体w(i,j,n)，i代表Inline测线序号，总测线数记为Nline；j代表Xline测线序号，总测线数记为Nxline；n代表单道地震记录的采样点序号，总采样点数记为N。
3.根据权利要求2所述的基于分数阶希尔伯特变换的薄砂体平面展布预测方法，其特征在于：所述步骤b中，令v＝-1.0,-0.8,-0.6,-0.4,-0.2,0.2,0.4,0.6,0.8,1，对三维数组w(i,j,n)的每个地震道记录进行分数阶希尔伯特变换，得到各个阶数的三维数组w(i,j,n)的分数阶希尔伯特变换的结果数组wv(i,j,n)，共10个三维数组。
4.根据权利要求3所述的基于分数阶希尔伯特变换的薄砂体平面展布预测方法，其特征在于：所述步骤b具体过程如下：
b1、令v＝-1.0；
b2、令i＝1,j＝1提取当前地震道数组s(n)，s(n)为一维数组，数组的长度是N；
b3、对目标一维数组进行s(n)做DFT，得到S(k)(k＝0,1,…,N-1)，其中k＝N/2,…,N-1对应负频率；
b4、对S(k)进行分频率计算：
b5、对y(k)进行逆DFT，取其实部即得到s(n)的分数阶希尔伯特变换sv(n)，并写入三维数组wv(i,j,n)；
b6、增加i,j的数值，重复步骤b2和步骤b5，直到i＝Nline,j＝Nxline；得到当前阶数v的分数阶希尔伯特变换的三维数组wv(i,j,n)；
b7、改变v的数值，重复步骤b2和步骤b6，直到v＝1.0，得到10个阶数对应的10个三维数组。
5.根据权利要求4所述的基于分数阶希尔伯特变换的薄砂体平面展布预测方法，其特征在于：所述步骤c中，根据井震标定分析得出砂体段发育的纵向位置，对步骤b的分数阶希尔伯特变换的结果数组wv(i,j,n)在纵向上采用等时切片方法提取各个砂体发育小层段的三维层段数组Sandl_v(i,j,nc)，l代表纵向上砂体发育小层段的序号，纵向上砂体发育小层段总数记为Nsand；nc代表当前层段内的地震记录的采样点序号，总采样点数记为Nc。
6.根据权利要求5所述的基于分数阶希尔伯特变换的薄砂体平面展布预测方法，其特征在于：所述步骤d中，令研究区内的实钻井数目是Nwell，统计步骤c中各个层段内的单井砂体发育情况，并写入二维数组P(a,l)，a代表实钻井序号，l代表纵向上发育砂体层段的序号，P(a,l)＝0代表该井在该层段砂体不发育，P(a,l)＝1代表该井在该层段砂体发育。
7.根据权利要求6所述的基于分数阶希尔伯特变换的薄砂体平面展布预测方法，其特征在于：所述步骤e中，令l＝1，计算当前砂体发育小层段内的三维数组Sandl_v(i,j,nc)(v＝-1.0,-0.8,-0.6,-0.4,-0.2,0.2,0.4,0.6,0.8,1)的均方根振幅属性，记为Amp_wv(i,j)(v＝-1.0,-0.8,-0.6,-0.4,-0.2,0.2,0.4,0.6,0.8,1)，并绘制均方根振幅平面图。
8.根据权利要求7所述的基于分数阶希尔伯特变换的薄砂体平面展布预测方法，其特征在于：所述步骤f中，根据实钻井在当前层段l的砂体的发育情况P(a,l)，计算步骤e中10个均方根振幅平面图的砂体平面展布预测符合率，选取预测符合率最高的阶数对应的均方根振幅平面图进行薄砂体平面展布预测。
9.根据权利要求8所述的基于分数阶希尔伯特变换的薄砂体平面展布预测方法，其特征在于：所述步骤g中，增加l的数值，重复步骤e和步骤f，直到l＝Nsand。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010419191.0A CN111598340B (zh) | 2020-05-18 | 2020-05-18 | 基于分数阶希尔伯特变换的薄砂体平面展布预测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010419191.0A CN111598340B (zh) | 2020-05-18 | 2020-05-18 | 基于分数阶希尔伯特变换的薄砂体平面展布预测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111598340A CN111598340A (zh) | 2020-08-28 |
| CN111598340B true CN111598340B (zh) | 2021-07-23 |
Family
ID=72187325
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010419191.0A Active CN111598340B (zh) | 2020-05-18 | 2020-05-18 | 基于分数阶希尔伯特变换的薄砂体平面展布预测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111598340B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114527505A (zh) * | 2020-11-23 | 2022-05-24 | 中国石油天然气股份有限公司 | 薄层识别方法、装置、设备及存储介质 |
| CN117153289B (zh) * | 2023-09-14 | 2024-04-05 | 大庆油田有限责任公司 | 一种适用于窄薄砂体的油层水淹程度预测方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102707315A (zh) * | 2012-06-05 | 2012-10-03 | 中国石油天然气股份有限公司 | 利用无序度地震属性预测碳酸盐岩油气藏优质储层的方法 |
| CN108535774A (zh) * | 2018-03-20 | 2018-09-14 | 中国科学院地质与地球物理研究所 | 一种利用可控震源激发地震信号快速识别震相的方法及装置 |
| CN110389381A (zh) * | 2018-04-23 | 2019-10-29 | 中国石油化工股份有限公司 | 一种基于地震属性的砂体储层预测方法及装置 |
| US10520627B2 (en) * | 2014-11-25 | 2019-12-31 | Cgg Services Sas | Modeling an elastic stiffness tensor in a transverse isotropic subsurface medium |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE102009010751A1 (de) * | 2009-02-26 | 2010-09-02 | Ehrhardt, Wilhelm, Prof. Dr.-Ing. | Verfahren zur verbesserten Vorhersage und Bewertung von Erdbeben und/oder seismischen Erschütterungen |
| CN103643945B (zh) * | 2013-11-26 | 2016-06-22 | 辽河石油勘探局 | 薄层岩性储层识别与水平井钻井跟踪方法 |
| CN106443781A (zh) * | 2016-09-08 | 2017-02-22 | 西南石油大学 | 一种预测薄层砂体分布的方法 |
| CN106526670B (zh) * | 2016-09-21 | 2019-04-19 | 中石化石油工程技术服务有限公司 | 一种碎屑岩储层中地震属性砂体空间分布描述及评价的方法 |
| CN110632660B (zh) * | 2019-08-27 | 2021-06-01 | 中国石油天然气股份有限公司 | 基于地震数据体的薄砂体表征方法及装置 |
-
2020
- 2020-05-18 CN CN202010419191.0A patent/CN111598340B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102707315A (zh) * | 2012-06-05 | 2012-10-03 | 中国石油天然气股份有限公司 | 利用无序度地震属性预测碳酸盐岩油气藏优质储层的方法 |
| US10520627B2 (en) * | 2014-11-25 | 2019-12-31 | Cgg Services Sas | Modeling an elastic stiffness tensor in a transverse isotropic subsurface medium |
| CN108535774A (zh) * | 2018-03-20 | 2018-09-14 | 中国科学院地质与地球物理研究所 | 一种利用可控震源激发地震信号快速识别震相的方法及装置 |
| CN110389381A (zh) * | 2018-04-23 | 2019-10-29 | 中国石油化工股份有限公司 | 一种基于地震属性的砂体储层预测方法及装置 |
Non-Patent Citations (2)
| Title |
|---|
| 分数阶希尔伯特变换及其在地震沉积学中的应用;宋洪勇;《中国优秀硕士学位论文全文数据库 基础科学辑》;20120415;全文 * |
| 分数阶希尔伯特变换在提高纵向分辨率中的应用;宋洪勇 等;《华北地震科学》;20131231;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111598340A (zh) | 2020-08-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109425896B (zh) | 白云岩油气储层分布预测方法及装置 | |
| Berne et al. | Pleistocene forced regressions and tidal sand ridges in the East China Sea | |
| CN1276377C (zh) | 使用虚拟源阵列对地下岩层成像的方法 | |
| CN107783187B (zh) | 一种将测井速度和地震速度结合建立三维速度场的方法 | |
| CN110231652B (zh) | 一种基于密度的含噪声应用空间聚类的地震相提取方法 | |
| CN111598340B (zh) | 基于分数阶希尔伯特变换的薄砂体平面展布预测方法 | |
| CN109541685B (zh) | 一种河道砂体识别方法 | |
| CN107065011A (zh) | 一种应用于陆相盆地储层反演的曲线频率融合方法 | |
| CN103901465A (zh) | 全息三维地震勘探观测系统的设计方法 | |
| CN108427140B (zh) | 一种用于小尺度缝洞储集体地震识别的方法 | |
| Bruno et al. | Groundwater characterization in arid regions using seismic and gravity attributes: Al Jaww Plain, UAE | |
| CN112764103A (zh) | 一种基于稀疏编码特征dbscan聚类地震相分析方法 | |
| CN112764110A (zh) | 一种基于限制波尔兹曼机特征编码的聚类地震相分析方法 | |
| CN109283577B (zh) | 一种地震层位标定方法 | |
| CN113109875A (zh) | 一种全波形速度场约束下的盐下碳酸盐岩储层反演方法 | |
| CN115857047B (zh) | 一种地震储层综合预测方法 | |
| CN110727027B (zh) | 多期次河道砂的精细刻画方法 | |
| CN116482757A (zh) | 基于浅剖井震联合的地层结构识别分析方法 | |
| CN113219531A (zh) | 致密砂岩气水分布的识别方法及装置 | |
| Becht et al. | Analysis of multi‐offset GPR data: A case study in a coarse‐grained gravel aquifer | |
| CN115327627A (zh) | 多信息融合的致密砂岩气分流河道表征方法及装置 | |
| CN113552623A (zh) | 曲流河沉积砂体边界的确定方法和确定装置 | |
| CN111366976A (zh) | 基于尺度导向的地震属性自适应中值滤波方法 | |
| Naseer et al. | Appraisal of reservoir quality for hydrocarbon-bearing deep-water Miocene sandstones incised valley, south-east Asian offshore Indus: An application of seismic attributes and instantaneous spectral porosity quantitative reservoir simulations | |
| CN113391350B (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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |