CN112346121A - 一种基于全波形的储集体地层分离处理方法 - Google Patents
一种基于全波形的储集体地层分离处理方法 Download PDFInfo
- Publication number
- CN112346121A CN112346121A CN202011215112.0A CN202011215112A CN112346121A CN 112346121 A CN112346121 A CN 112346121A CN 202011215112 A CN202011215112 A CN 202011215112A CN 112346121 A CN112346121 A CN 112346121A
- Authority
- CN
- China
- Prior art keywords
- data
- energy
- reservoir
- stratum
- counting
- 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
- 238000000926 separation method Methods 0.000 title claims abstract description 19
- 238000003672 processing method Methods 0.000 title description 4
- 238000000034 method Methods 0.000 claims abstract description 47
- 238000001228 spectrum Methods 0.000 claims abstract description 44
- 238000012545 processing Methods 0.000 claims abstract description 30
- 238000011160 research Methods 0.000 claims abstract description 20
- 238000010606 normalization Methods 0.000 claims abstract description 7
- 238000005070 sampling Methods 0.000 claims description 42
- 230000015572 biosynthetic process Effects 0.000 claims description 20
- 230000002159 abnormal effect Effects 0.000 claims description 16
- 230000003595 spectral effect Effects 0.000 claims description 12
- 230000002087 whitening effect Effects 0.000 claims description 11
- 238000004458 analytical method Methods 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 7
- 230000004044 response Effects 0.000 claims description 5
- 238000009825 accumulation Methods 0.000 claims description 2
- 239000011800 void material Substances 0.000 claims description 2
- 238000004519 manufacturing process Methods 0.000 abstract description 10
- 239000011148 porous material Substances 0.000 abstract description 9
- 238000005516 engineering process Methods 0.000 description 17
- 235000011850 desserts Nutrition 0.000 description 7
- 239000012530 fluid Substances 0.000 description 7
- 230000008859 change Effects 0.000 description 6
- 235000021185 dessert Nutrition 0.000 description 6
- 238000011161 development Methods 0.000 description 6
- 235000009508 confectionery Nutrition 0.000 description 5
- 238000013508 migration Methods 0.000 description 4
- 230000005012 migration Effects 0.000 description 4
- 239000003208 petroleum Substances 0.000 description 4
- 238000003384 imaging method Methods 0.000 description 3
- 238000012937 correction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 238000005553 drilling Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000000945 filler Substances 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 235000015055 Talinum crassifolium Nutrition 0.000 description 1
- 244000010375 Talinum crassifolium Species 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000000712 assembly Effects 0.000 description 1
- 238000000429 assembly Methods 0.000 description 1
- 239000011324 bead Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005282 brightening Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 238000011031 large-scale manufacturing process Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001502 supplementing effect Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
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
- G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
-
- 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/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/30—Noise handling
- G01V2210/32—Noise reduction
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/53—Statics correction, e.g. weathering layer or transformation to a datum
-
- 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/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
-
- 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/63—Seismic attributes, e.g. amplitude, polarity, instant phase
-
- 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/64—Geostructures, e.g. in 3D data cubes
- G01V2210/642—Faults
-
- 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/64—Geostructures, e.g. in 3D data cubes
- G01V2210/646—Fractures
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
本发明公开了一种基于全波形的储集体地层分离处理方法，包括以下步骤：第一步，完成研究区内地震资料的能量归一化处理；第二步，选取研究区域的核心地层数据，并把核心地层数据分为四个频率段，第三步，选择第二步中0‑10Hz的第一分频数据体，统计单个椭圆频谱的能量，第四步，用第三步的方法，分别其他分频数据体能量统计；第五步，进行相似系数的计算，用椭圆频谱间的波形进行处理，并统计出真实地层背景；第六步，细分为多种不同尺寸规模的储集体；本发明能够通过精确统计地层信息，并从地震资料中去掉，留下的信息为大的裂缝、孔隙的信息，提高有效油气藏的准确性，提高生产效率。
Description
技术领域
本发明涉及石油物理勘探技术领域，具体为一种基于全波形的储集体地层分离处理方法。
背景技术
在石油地球物理勘探开发中发现，油气的储存空间和有效的油气运移通道，简称“储集体”(储层集合体)，是判断油气的基础资料；储集体的寻找过程，经历了递进及平行分支的亮点、平点、甜点及串珠再到断溶体、深部油源断裂等“储集体”的关键技术，在油气勘探开发的进程中，这些技术发挥了重要的作用。
以上所有技术都是基于资料的相对振幅的变化，因其在地震剖面上相对能量较强(发亮)而得名，也就是反射波振幅比较强，并且证明和油气的相关性比较大，最大的优点是可以在资料上可以直观看到，后因技术的停滞而受阻；之后的技术更多的通过各种反演、地质分析、测井分析等技术手段，间接成图进行储层预测等，但也同样是用到了资料的强振幅信息。目前最大的问题是，当地层的能量接近或大于亮点的能量时，其误差就开始变大，人们也尝试用衰减地层能量的方法应用于生产，都是因为地层的存在而效果不佳。
目前有利于发现大型油田或大规模提高产量的有以下两种主要技术：①断溶体储层技术。发表在石油地球物理勘探上的《塔里木盆地顺北地区断溶体地震反射特征与识别》一文中对断溶体描述和识别(见图1和图2)，是亮点等技术的升级，是地球物理最本质的方法，是通过多个亮点的组合，并结合断层确定油气的存储、运移及油源通道来描述评价优质储集体的理论，未来的发展空间较大，部分资料信息可以从资料上直接看到。
②甜点储层技术。甜点是有大段含油地层中，油气的分布在横向上和纵向上都不连片，成点状分布，被称为sweetheart。在西得克萨斯州特拉华和米德兰盆地致密油藏中得到应用。通过收集整理134口井的井位信息、地震和测井、储层状况等资料，建立人工神经网络模型，最终确定了未来勘探的潜力区(见图3)，明确、客地指出了未勘探地区的甜点储集体分布，避免了地震解释造成的主观判断错误，及时发现成熟探区遗漏的储集体；可直接排除低产区块，大幅提高工作效率。
其中，①和早期的亮点、平点一样，部分资料可以直接看到由亮点所组成的断溶体，因为资料处理技术还没有突破，不能去掉资料背景的影响，信噪比高的区域，无法直接识别断溶体的存在，只能用推理的方法，最重要的是，深部大断裂是通过推断而得到，其断裂是否为通道也很难证明。如果能够通过技术减弱或剥离掉地层，不但亮点可以分离出来，则由亮点所组成的断溶体也可以清楚看到，并且可以识别大断裂是否被流体所充填，进一步识别是否能成为通道；方法②中甜点成像识别，除美国之外，几乎没有成功的案例，原因是高质量的地震资料使地层的能量更强，或甜点能量本身相对较弱，因为地层的存在，常规的方法基本无法获得甜点的信息，而通过测井、钻进信息进行地质分析后再识别，则生产效率太降低。
对于深部探井的设计开发，因为没有相应的井资料信息，目前的多数方法很难预测准确，增加了深井钻探的风险。
对于储层的分析，是地质解释人员从最终的偏移资料中去寻找。如果处理资料时，在去噪时就已经把储集体去掉了，那么一个工期1-3年的时间大家都在做无用功，并且容易得到错误的结论。如果用储集体成像技术进行监控，则可能会避免出现这个问题。
目前所有储层预测，弹性波反演，流体检测等技术，都是以能量的异常为依据，但所有的计算，都把地层的能量计算进去再找差异，无法实现上下强弱地层间的储层关系对比，从而精度降低。同时在多次进行预测工作时，储层的位置在变，储层区域的位置和大小也在变化，油气的运移的通道，特别是纵向和侧向的通道等空间关联，基本没有描述。
发明内容
本发明的目的在于提供一种基于全波形的储集体地层分离处理方法，通过精确统计地层信息，并从地震资料中去掉，留下的信息为大的裂缝、孔隙的信息，其中的亮点、甜点等和油气相关的信息等可以直接看到，可以识别断裂与甜点的关系，进一步提高有效油气藏的准确性，提高生产效率。
为实现上述目的，本发明提供如下技术方案：一种基于全波形的储集体地层分离处理方法，包括以下步骤：
第一步，完成研究区内地震资料的能量归一化处理；
从研究区内地震资料中选取时间偏移或深度偏移的道集或叠加资料，进行谱白化处理；
第二步，选取研究区的核心地层数据，并把核心地层数据分为四个频率段，选取0-10Hz的核心地层数据作为第一分频数据体，11-20Hz的核心地层数据作为第二分频数据体，21-30Hz的核心地层数据作为第三分频数据体，31-50Hz的核心地层数据作为第四分频数据体；
第三步，选择第二步中的第一分频数据体，分析第一分频数据体中研究区内的地层角度，选择第一分频数据体中研究区内地层角度最大的地层，从水平方向开始进行采样，首先统计一个采样点平均椭圆频谱能量，然后对下一个采样点进行统计，最后完成研究区每个采样点的统计和记录；
第四步，用第三步的方法，分别完成对11-20Hz的第二分频数据体，21-30Hz的第三分频数据体，31-50Hz的第四分频数据体中每个采样点的平均椭圆频谱能量的统计和记录；
第五步，利用第三步和第四步的得到的数据，进行相似系数的计算，利用每个采样点平均椭圆频谱能量间的波形进行处理，利用椭圆频谱能量统计异常地质体在某个方向延伸，并统计出真实地层背景；
第六步，细分为多种不同尺寸规模的储集体；
根据不同的构造特征和储集体的形状，先对真实地层背景的不同地区进行地质分析，并根据真实地质体的固有频率来从第五步得到的数据中寻找最接近真实地质体的最强地震响应，优化和储层不吻合的信息，最终得到多种不同尺寸规模的储集体。
在上述任一方案中优选的是，第一步中，谱白化处理用于均衡非全频带处理的资料，使不同级别、不同频率的储集体所对应的强反射信息，以一定的能量存在，使程序能进行捕捉和识别。
在上述任一方案中优选的是，在第三步中，统计一个采样点平均椭圆频谱能量的步骤如下：以采样点为椭圆圆心，并计算出该采样点距离W个相邻道中最高点和最低点之间的时差，作为椭圆的纵轴长度，将相邻道的数量W作为椭圆的横轴长度，统计以该采样点单个椭圆频谱的能量；保持椭圆的圆心位置不变，每旋转1度后，再次统计旋转后椭圆的频谱能量，统计完该采样点的1-360度的频谱能量后，取1-360度的频谱能量平均值作为该采样点的平均椭圆频谱能量。
在上述任一方案中优选的是，在第五步中，计算相似系数的公式如下：
其中，r为相似系数，x(n)、y(n)是相邻两个椭圆的信号，n为采样点的数量，j为沿长轴方向的采样点距离，其中j＝0，1，...。
在上述任一方案中优选的是，对于已知的信号x(n)、y(n)，它们的能量是确定的，则采用下述公式描述：
在上述任一方案中优选的是，对于椭圆频谱的进行波形相似统计，即对同一频率段内的能量做归一化处理，归一化处理公式如下：
其中，x(n)，y(n)是相邻两个椭圆的信号，ρxy是x(n)、y(n)的互相关系数，ρxy范围是[-1～+1]。
在上述任一方案中优选的是，根据互相关系数的统计结果，反映相邻椭圆频谱波形一致性的级别，是其中0代表不相关，-1代表相位相反，1代表相关性最好；计算中互相关系数大于0.5的为相关数据，用来统计背景能量，小于0.5的为不相关数据，用来统计异常体能量，不参与地层的能量统计；相关部分的平均能量和不相关部分的平均能量倍数m，其计算公式为：
m＝Aρ0/Aρ1；
其中Aρ0为不相关数据能量平均值，Aρ1为相关数据的能量平均值。
在上述任一方案中优选的是，从地震数据中统计相关数据能量，初步识别地层信息并剥离出去，非地层异常体的能量相对增加，多次进行迭代后，则充填低密度的空隙体组合的真实能量就相对加强，当m>1时，背景能量就弱于异常体能量；在处理的过程中，对所选择的椭圆频谱，进行调整和试验，使其大小要比储集体大，防止一个地质体分为两个或多个地质体的假象。
在上述任一方案中优选的是，在第六步中，将多种不同尺寸规模的储集体数据细分为3个等级，大尺寸储集体数据、中尺寸储集体数据、小尺寸储集体数据。
在上述任一方案中优选的是，同时对不同固有频率响应的资料进行展示时，对相同等级的储集体数据进行数据叠合或进行图形累加，以及和储层数据的叠加。
与现有技术相比，本发明的有益效果是：
1、本发明能够通过精确统计地层信息，并从地震资料中去掉，留下的信息为大的裂缝、孔隙的信息，其中的亮点、甜点等和油气相关的信息等可以直接看到，可以识别断裂与甜点的关系，进一步提高有效油气藏的准确性，提高生产效率，是目前在深层、老地层中寻找油气的较好方法。
2、本发明从资料处理开始的第一步就可以获得可能是储层的信息，其中每做一步，每选择一个参数，就可以看到初始储集体的变化，如深部大断裂等，从而实现储层的保幅处理。还可以与球面补偿、去面波、去随机噪音、反褶积、拼接处理、剩余静校正、速度分析、偏移成像技术等技术进行结合应用，随时检查储集体的变化。
3、本发明提供一套解释性处理方案，通过解释性处理，把地震资料中的地层剥离出去，留下非地层信息的资料，其中含有孔洞、孔隙、裂缝、断层、破碎带等地质体，以及地质体内所含充填物的信息，相对地层来说就是异常体；处理结果不只是得到地层中的甜点，同时可以获得隐藏在地层周围的甜点，可以更好的识别出储集体的连接通道和空间关系信息等，是有效储层空间的扩展。同时可以区分断裂是否封闭还是被流体所充填，最终实现工区内储集体的可视化数据，为油田开发，新区开发，寻找优质储层和井位部署提供更加可靠的依据。
4、本发明采用椭圆频谱的统计方法，其计算精度高，能够精确的找到异常体，提高工作效率。
5、本发明采用相似系数的计算，利用每个采样点平均椭圆频谱能量间的波形进行处理，利用椭圆频谱能量统计异常地质体在某个方向延伸，能够精确计算出真实地层背景，并快速去掉地层资料，其工作效率非常高，且准确度也很高。
本发明附加的方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明的实践了解到。
附图说明
图1为低信噪比资料断溶体剖面；
图2为高信噪比资料断溶体剖面；
图3为未来勘探潜力区预测的甜点图；
图4(a)为一个井区地层资料的原始资料图；
图4(b)为本发明对图4(a)进行第一至第四步以后的储集体的资料图；
图5(a)为本发明对图4(b)进行第五步后的储集体的资料图；
图5(b)为本发明对图5(a)进行第六步后的储集体的资料图；
图6(a)为一工区储集体通道的资料图；
图6(b)为本发明对图6(a)进行处理后的储集体通道的资料图；
图6(c)为实际井的产量资料图；
图7为本发明的流程图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
在本发明中，除非另有明确的规定和限定，术语“安装”、“相连”、“连接”、“固定”等术语应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以根据具体情况理解上述术语在本发明中的具体含义。
为了解决背景技术中所提到的技术问题，本发明提供一种基于全波形的储集体地层分离处理方法，主要是利用地震资料中充满低密度流体，大于一个面元的孔洞体、裂缝体、断裂带等地质体，分析对比从低频到高频的波形与地层信息的波形的差异，主要依据地层的频率从低频到高频，其能量相对稳定，通过计算每个频率段地层的平均能量，并从数据中减去，则剩下的就是要得到的非地层信息。所得数据不包括地层(也可以人为包含部分地层)信息，以储层信息为主，包括隐藏在地层中的亮点、平点、甜点等优质储层的储集体。
如图7所示，本发明提供一种基于全波形的储集体地层分离处理方法，包括以下步骤：
第一步，完成研究区内地震资料的能量归一化处理；
从研究区内地震资料中选取时间偏移或深度偏移的道集或叠加资料，进行谱白化处理；用于均衡非全频带处理的资料，使不同级别、不同频率的储集体所对应的强反射信息，以一定的能量存在，使程序能进行捕捉和识别。
谱白化处理(Spectral Whitening)是地震勘探数据处理的一种常用方法,可以展宽信号的频谱提高信号分辨能力,是一种“纯振幅”的滤波过程。
谱白化处理是在假设反射系数序列的谱是白的基础上的，即在全频带内其谱值都是1。根据这个假设，可以将反摺积以后(这里是指偏移或叠加)的地震道分成许多频带比较窄的地震道，然后，分别对分解后的各道能量进行均衡处理，将均衡处理后的各道再组合就得到谱白化后的地震道。
虽然谱白化处理的技术手段是一种常用的方法，但是在本技术领域中，不做谱白化的处理，更多较弱的信息就无法识别，因为资料处理技术还没有突破，不能去掉资料背景的影响，特别是信噪比高的区域，对于识别断溶体正常的方法是用推理的方法，谱白化后的资料，识别断溶体就更加困难，并且当地层的能量接近或大于亮点的能量时，其误差就开始变大，人们也尝试用衰减地层能量的方法应用于生产，都是因为地层的存在而效果不佳。
基于此，本方案采用下面的步骤来去除地层的能量，不但亮点可以分离出来，则由亮点所组成的断溶体也可以清楚看到，并且可以识别大断裂是否被流体所充填，进一步识别是否能成为通道。
第二步，选取研究区的核心地层数据，并把核心地层数据分为四个频率段，选取0-10Hz的核心地层数据作为第一分频数据体，11-20Hz的核心地层数据作为第二分频数据体，21-30Hz的核心地层数据作为第三分频数据体，31-50Hz的核心地层数据作为第四分频数据体；
第三步，选择第二步中的第一分频数据体，分析第一分频数据体中研究区内的地层角度，选择第一分频数据体中研究区内地层角度最大的地层，从水平方向开始进行采样。
首先统计一个采样点平均椭圆频谱能量，以采样点为椭圆圆心，并计算出该采样点距离W个相邻道中最高点和最低点之间的时差(相邻道的数据可以根据实际情况来选择，对于一般的地震资料而言，选取6个相邻道最佳，此外，如果不是水平地层，则用全波形的时差，如果是水平地层，则用半个波形的时差就行)，作为椭圆的纵轴长度，将相邻道的数量W作为椭圆的横轴长度，统计以该采样点单个椭圆频谱的能量；保持椭圆的圆心位置不变，每旋转1度后，再次统计旋转后椭圆的频谱能量，统计完该采样点的1-360度的频谱能量后，取1-360度的频谱能量平均值作为该采样点的平均椭圆频谱能量。
然后对下一个采样点进行统计，最后完成研究区每个采样点的统计和记录。
第四步，用第三步的方法，分别完成对11-20Hz的第二分频数据体，21-30Hz的第三分频数据体，31-50Hz的第四分频数据体中每个采样点的平均椭圆频谱能量的统计和记录；
如图4(a)-图4(b)所示，通过对比可知，图4(a)中的井区地层资料经过第一步至第四步处理后，可以将地震资料中的地层剥离出去，留下非地层信息的资料，使得含有孔洞、孔隙、裂缝、断层、破碎带等地质体，以及地质体内所含充填物的信息，相对地层来说就是异常体，能够清楚的显现出来如图4(b)所示。
第五步，利用第三步和第四步的得到的数据，进行相似系数的计算，用椭圆频谱间的波形进行处理，利用椭圆统计异常地质体在某个方向延伸，并统计出真实地层背景；
计算相似系数的公式如下：
其中，r为相似系数，x(n)、y(n)是相邻两个椭圆的信号，n为采样点的数量，j为沿长轴方向的采样点距离，其中j＝0，1，...。
对于已知的信号x(n)、y(n)，它们的能量是确定的，则采用下述公式描述：
对于椭圆频谱的进行波形相似统计，即对同一频率段内的能量做归一化处理，归一化处理公式如下：
其中，x(n)，y(n)是相邻两个椭圆的信号，ρxy是x(n)、y(n)的互相关系数，ρxy范围是[-1～+1]。
根据互相关系数的统计结果，反映相邻椭圆频谱波形一致性的级别，是其中0代表不相关，-1代表相位相反，1代表相关性最好；计算中互相关系数大于0.5的为相关数据，用来统计背景能量，小于0.5的为不相关数据，用来统计异常体能量，不参与地层的能量统计；相关部分的平均能量和不相关部分的平均能量倍数m，其计算公式为：
m＝Aρ0/Aρ1；
其中Aρ0为不相关数据能量平均值，Aρ1为相关数据的能量平均值。
从地震资料数据中减去背景能量，即减去不相关数据能量，初步识别地层信息并剥离出去，非地层异常体的能量相对增加，多次进行迭代后，则充填低密度的空隙体组合的真实能量就相对加强，当m>1时，背景能量就弱于异常体能量；在处理的过程中，对所选择的椭圆频谱，进行调整和试验，使其大小要比储集体大，防止一个地质体分为两个或多个地质体的假象。
第六步，细分为多种不同尺寸规模的储集体；
根据不同的构造特征和储集体的形状，先对真实地层背景的不同地区进行地质分析，并根据真实地质体的固有频率来从第五步得到的数据中寻找最接近真实地质体的最强地震响应(这点是基于下面的原理来实现的，当外界的力作用到某一地质体时，外力的频率和地层的固有频率一致时，才能产生共振，产生最大的振幅)，优化和储层不吻合的信息，最终得到多种不同尺寸规模的储集体。
将多种不同尺寸规模的储集体数据细分为多个等级，对于一般资料而言，可以分为3个等级，细分为大尺寸储集体数据、中尺寸储集体数据、小尺寸储集体数据，对于特殊的资料而言可以根据特殊资料来分级。
如图5(a)-图5(b)，将图4(b)进行相似系数计算以及储集体细分以后，能够将剩下的孔洞、裂缝等充填低密度流体的孔隙的组合体，直接显示能够成为亮点、平点、甜点、断溶体等优质储集体。帮助人们对油气的识别，可以不用考虑地层能量的影响，提高生产效率和准确率。
通过图6(a)-图6(c)对比可知，采用本发明的基于全波形的储集体地层分离处理方法处理后的储集体通道余实际井的产量是完全吻合的。本发明从资料处理开始的第一步就可以获得可能是储层的信息，其中每做一步，每选择一个参数，就可以看到初始储集体的变化，如深部大断裂等，从而实现储层的保幅处理。其中在监控处理过程中的球面补偿、去面波、去随机噪音、反褶积、拼接处理、剩余静校正、速度分析、偏移成像技术应用时，随时检查储集体的变化，可以将从地震资料中剥离地层信息，识别较大孔隙中是否有甜点等的存在，还能识别深部大断裂的流通性。
本发明属于石油地球物理勘探领域，具体涉及从地层中分离优质储层的方法，是对目前储层预测、寻找有利储层体的有效补充和进一步完善；能够把地层从地震资料中剥离出去或减弱，然后就只剩下孔洞、裂缝等充填低密度流体的孔隙的组合体，直接显示能够成为亮点、平点、甜点、断溶体等优质储集体。帮助人们对油气的识别，可以不用考虑地层能量的影响，提高生产效率和准确率。对于深部油气和浅层油气的技术应用，该方法的应用没有任何区别，准确性是一样的。
在本发明的描述中，需要理解的是，指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
尽管已经示出和描述了本发明的实施例，对于本领域的普通技术人员而言，可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型，本发明的范围由所附权利要求及其等同物限定。
Claims (10)
1.一种基于全波形的储集体地层分离处理方法，其特征在于：包括以下步骤：
第一步，完成研究区内地震资料的能量归一化处理；
从研究区内地震资料中选取时间偏移或深度偏移的道集或叠加资料，进行谱白化处理；
第二步，选取研究区的核心地层数据，并把核心地层数据分为四个频率段，选取0-10Hz的核心地层数据作为第一分频数据体，11-20Hz的核心地层数据作为第二分频数据体，21-30Hz的核心地层数据作为第三分频数据体，31-50Hz的核心地层数据作为第四分频数据体；
第三步，选择第二步中的第一分频数据体，分析第一分频数据体中研究区内的地层角度，选择第一分频数据体中研究区内地层角度最大的地层，从水平方向开始进行采样，首先统计一个采样点平均椭圆频谱能量，然后对下一个采样点进行统计，最后完成研究区每个采样点的统计和记录；
第四步，用第三步的方法，分别完成对11-20Hz的第二分频数据体，21-30Hz的第三分频数据体，31-50Hz的第四分频数据体中每个采样点的平均椭圆频谱能量的统计和记录；
第五步，利用第三步和第四步的得到的数据，进行相似系数的计算，利用每个采样点平均椭圆频谱能量间的波形进行处理，利用椭圆频谱能量统计异常地质体在某个方向延伸，并统计出真实地层背景；
第六步，细分为多种不同尺寸规模的储集体；
根据不同的构造特征和储集体的形状，先对真实地层背景的不同地区进行地质分析，并根据真实地质体的固有频率来从第五步得到的数据中寻找最接近真实地质体的最强地震响应，优化和储层不吻合的信息，最终得到多种不同尺寸规模的储集体。
2.根据权利要求1所述的基于全波形的储集体地层分离处理方法，其特征在于：第一步中，谱白化处理用于均衡非全频带处理的资料，使不同级别、不同频率的储集体所对应的强反射信息，以一定的能量存在，使程序能进行捕捉和识别。
3.根据权利要求1所述的基于全波形的储集体地层分离处理方法，其特征在于：在第三步中，统计一个采样点平均椭圆频谱能量的步骤如下：以采样点为椭圆圆心，并计算出该采样点距离W个相邻道中最高点和最低点之间的时差，作为椭圆的纵轴长度，将相邻道的数量W作为椭圆的横轴长度，统计以该采样点单个椭圆频谱的能量；保持椭圆的圆心位置不变，每旋转1度后，再次统计旋转后椭圆的频谱能量，统计完该采样点的1-360度的频谱能量后，取1-360度的频谱能量平均值作为该采样点的平均椭圆频谱能量。
4.根据权利要求1所述的基于全波形的储集体地层分离处理方法，其特征在于：在第五步中，计算相似系数的公式如下：
其中，r为相似系数，x(n)、y(n)是相邻两个椭圆的信号，n为采样点的数量，j为沿长轴方向的采样点距离，其中j＝0，1，...。
5.根据权利要求4所述的基于全波形的储集体地层分离处理方法，其特征在于：对于已知的信号x(n)、y(n)，它们的能量是确定的，则采用下述公式描述：
6.根据权利要求4所述的基于全波形的储集体地层分离处理方法，其特征在于：对于椭圆频谱的进行波形相似统计，即对同一频率段内的能量做归一化处理，归一化处理公式如下：
其中，x(n)，y(n)是相邻两个椭圆的信号，ρxy是x(n)、y(n)的互相关系数，ρxy范围是[-1～+1]。
7.根据权利要求6所述的基于全波形的储集体地层分离处理方法，其特征在于：根据互相关系数的统计结果，反映相邻椭圆频谱波形一致性的级别，是其中0代表不相关，-1代表相位相反，1代表相关性最好；计算中互相关系数大于0.5的为相关数据，用来统计背景能量，小于0.5的为不相关数据，用来统计异常体能量，不参与地层的能量统计；相关部分的平均能量和不相关部分的平均能量倍数m，其计算公式为：
m＝Aρ0/Aρ1；
其中Aρ0为不相关数据能量平均值，Aρ1为相关数据的能量平均值。
8.根据权利要求7所述的基于全波形的储集体地层分离处理方法，其特征在于：从地震数据中统计相关数据能量，初步识别地层信息并剥离出去，非地层异常体的能量相对增加，多次进行迭代后，则充填低密度的空隙体组合的真实能量就相对加强，当m>1时，背景能量就弱于异常体能量；在处理的过程中，对所选择的椭圆频谱，进行调整和试验，使其大小要比储集体大，防止一个地质体分为两个或多个地质体的假象。
9.根据权利要求1所述的基于全波形的储集体地层分离处理方法，其特征在于：在第六步中，将多种不同尺寸规模的储集体数据细分为3个等级，大尺寸储集体数据、中尺寸储集体数据、小尺寸储集体数据。
10.根据权利要求9所述的基于全波形的储集体地层分离处理方法，其特征在于：同时对不同固有频率响应的资料进行展示时，对相同等级的储集体数据进行数据叠合或进行图形累加，以及和储层数据的叠加。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011215112.0A CN112346121B (zh) | 2020-11-04 | 2020-11-04 | 一种基于全波形的储集体地层分离处理方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011215112.0A CN112346121B (zh) | 2020-11-04 | 2020-11-04 | 一种基于全波形的储集体地层分离处理方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112346121A true CN112346121A (zh) | 2021-02-09 |
| CN112346121B CN112346121B (zh) | 2023-09-12 |
Family
ID=74356433
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011215112.0A Active CN112346121B (zh) | 2020-11-04 | 2020-11-04 | 一种基于全波形的储集体地层分离处理方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112346121B (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113093277A (zh) * | 2021-04-08 | 2021-07-09 | 中石化石油工程技术服务有限公司 | 基于时效的隐蔽油气藏目标增强解释性处理方法 |
| CN117908128A (zh) * | 2024-03-08 | 2024-04-19 | 成都星辉科瑞科技有限责任公司 | 一种波数域优势频带识别岩性体断溶体边界的检测方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2713185A1 (en) * | 2012-09-28 | 2014-04-02 | Dal Moro Giancarlo Ditta Individuale | Method and apparatus to detect and analyze seismic signals |
| CN105093313A (zh) * | 2015-07-07 | 2015-11-25 | 中国地质大学（北京） | 一种岩溶型油藏单井油气产能预测方法及装置 |
-
2020
- 2020-11-04 CN CN202011215112.0A patent/CN112346121B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2713185A1 (en) * | 2012-09-28 | 2014-04-02 | Dal Moro Giancarlo Ditta Individuale | Method and apparatus to detect and analyze seismic signals |
| CN105093313A (zh) * | 2015-07-07 | 2015-11-25 | 中国地质大学（北京） | 一种岩溶型油藏单井油气产能预测方法及装置 |
Non-Patent Citations (1)
| Title |
|---|
| 刘浩杰;王延光;孟宪军;夏吉庄;: "高勘探程度区储层描述技术的探索与思考", 油气地球物理, no. 01 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113093277A (zh) * | 2021-04-08 | 2021-07-09 | 中石化石油工程技术服务有限公司 | 基于时效的隐蔽油气藏目标增强解释性处理方法 |
| CN113093277B (zh) * | 2021-04-08 | 2022-06-24 | 中石化石油工程技术服务有限公司 | 基于时效的隐蔽油气藏目标增强解释性处理方法 |
| CN117908128A (zh) * | 2024-03-08 | 2024-04-19 | 成都星辉科瑞科技有限责任公司 | 一种波数域优势频带识别岩性体断溶体边界的检测方法 |
| CN117908128B (zh) * | 2024-03-08 | 2024-07-30 | 成都星辉科瑞科技有限责任公司 | 一种波数域优势频带识别岩性体断溶体边界的检测方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112346121B (zh) | 2023-09-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101738639B (zh) | 提高岩石裂缝参数精度的方法 | |
| CN102466815B (zh) | 三叠系碎屑岩油气藏识别方法 | |
| CN103728659B (zh) | 一种提高地下岩溶探测精度的方法 | |
| CN109655896B (zh) | 岩溶古河道识别及刻画的方法及系统 | |
| CN111158044B (zh) | 潜山裂缝储集体油藏预测方法 | |
| CN107356965B (zh) | 基于加权叠加噪音压制策略的反射系数反演储层预测方法 | |
| CN104360382A (zh) | 一种利用叠后地震数据进行油气检测的方法 | |
| CN114114459B (zh) | 一种相控约束下的深层-超深层碳酸盐岩薄储层预测方法 | |
| CN110658556A (zh) | 一种碳酸盐岩走滑断裂破碎带判识与评价的地震技术组合方法 | |
| CN112346121B (zh) | 一种基于全波形的储集体地层分离处理方法 | |
| CN105240006B (zh) | 一种适用于火山岩储层的油水层识别方法 | |
| CN108254783A (zh) | 一种基于时频分析的叠后地震流体识别方法 | |
| CN111856566B (zh) | 湖相滩坝砂体中薄储层预测方法及装置 | |
| CN112835124B (zh) | 基于成像测井和阵列声波测井资料的裂缝有效性评价方法 | |
| CN113279746A (zh) | 一种套管变形风险区域的确定方法及应用 | |
| CN104142517A (zh) | 一种利用地震数据动态频谱属性的油气检测方法 | |
| CN106353806A (zh) | 一种检测生物礁储层连通性的方法 | |
| CN109425889B (zh) | 一种用于刻画古岩溶暗河的方法 | |
| CN115857047B (zh) | 一种地震储层综合预测方法 | |
| CN113514884A (zh) | 一种致密砂岩储层预测方法 | |
| CN112698397A (zh) | 一种台盆区走滑断裂空腔储层刻画方法 | |
| CN113093275B (zh) | 曲流河及曲流河三角洲油田钻井成功率提高方法及装置 | |
| CN111352154B (zh) | 一种基于宽方位地震的储层预测方法 | |
| Thompson | Induced fracture detection in the Barnett Shale, Ft. Worth Basin, Texas | |
| CN111812743B (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 |