CN1073705C - 地震-测井多信息储层数据采集处理方法 - Google Patents
地震-测井多信息储层数据采集处理方法 Download PDFInfo
- Publication number
- CN1073705C CN1073705C CN 96118790 CN96118790A CN1073705C CN 1073705 C CN1073705 C CN 1073705C CN 96118790 CN96118790 CN 96118790 CN 96118790 A CN96118790 A CN 96118790A CN 1073705 C CN1073705 C CN 1073705C
- Authority
- CN
- China
- Prior art keywords
- centerdot
- parameter
- well logging
- information
- formula
- 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
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种预测地层结构中地下地质体的储存信息的处理方法。其特征在于:(1)将多个传感器和弹药埋于地层内并引爆弹药或用重锤冲击地层,由埋入地层内的传感器采集到相关数据,(2)将由传感器采集到的数据的地震、测井、地质资料数据输入,(3)计算地震剖面上空间距离为△X的数值差的期望值E丨BH(X+△X)-BH(X)丨2;(4)确定尺度极限参数,计算丨△X丨min、丨△X丨max;画出分维图,即logE丨BH(X+△X)-BH(X)丨2相对log丨△X丨的曲线;(5)计算参数H和地震数据正态分布的标准差δ;(6)得到各种测井曲线剖面。
Description
本发明涉及一种预测地层结构中地下地质体的储存信息的处理方法。
目前石油或天然气的地质勘探是通过地震或重锤法来采集地层下方的多种信息,它们都将多个传感器和弹药埋于地层内并引爆弹药或用重锤冲击地层,由埋入地层内的传感器采集到相关数据,然后对采集到的数据进行处理(也有称作反演处理法),我国在这方面进行了大量的研究工作,并取得了一定的成效。近年来引进和研制的slim、Delog,这些方法在处理过程中没有充分利用测井资料的高频信息,因此分辨率较低。所以叠后地震道反演技术存在的问题有:纵、横向分辨率低;反演结果的多解性严重;横向外推精度不高;多井处理难于闭合;三维反演实现困难;多参数综合反演方法尚未成熟。
本发明的目的是提供一种专门针对目前叠后地震道反演技术存在的问题,采用“局部最优化算法”和WFBR方法,充分利用地震、测井和地质概念进行综合反演的地震-测井多信息储层数据采集处理方法。
本发明是最优化算法和WFBR方法的结合,首次提出地震波场是地下地质体微观波动的结果总和。地震波场是对地下地质体特征的总体反映,测井是对地下地质体特征的局部反映,地震和测井都包含地下地质体的储层信息。其特点是利用地震、测井、地质资料,采用信息优化和WFBR方法,处理出各种测井曲线剖面,然后将多种信息综合,得到储层参数剖面。有别于传统处理方法。传统方法是:利用地震、测井、地质资料以褶积模型为基础,采用优化相关外推方法处理出波阻抗剖面,然后利用经验公式或量板,计算出储层参数剖面。本发明的实现过程是采用滑动小时窗,首先计算地震剖面上空间距离为ΔX的数值差的期望值E|BH(X+ΔX)-BH(X)|2;然后确定尺度极限参数|ΔX|min,|ΔX|max;最后计算参数H和地震数据正态分布的标准差б。分数Brown运动BH(t)是一非平稳的具有均值为零的Gauss随机函数,其定义如下:
式中BH(t)是通常的Brown运动,0<H<1,当
时,FBM就是Brown运动,Pentland把FBM扩展到高维情况,给出了分数布郎随机场定义,设:X,ΔX∈R2,0≤H<1,F(y)是均值为零的高斯随机函数,Pr(·)表示概率密度,‖·‖表示范数,若随机场BH(X)满足:
则,BH(X)为分数布郎随机场(FBR),H参数反映了地震剖面的粗糙度,由H参数值可得地震剖面的分形维数D为
D=DT+1-H (3)
式中,DT为地震剖面的拓扑维数。
BH(X)具有如下性质:
E|BH(X+ΔX)-BH(X)|2=E|BH(X+1)-BH(X)|2·‖ΔX‖2H (4)
利用公式4即可方便地计算H参数。确定尺度极限参数|ΔX|min,|ΔX|max;如果地震剖面是完全理想分形的,则在任何尺度下分维均保持为常数,但实际地震剖面并不是完全理想分形的,所以要确定一个尺度范围,在此范围分维保持常数;确定方法为:画出分维图,即:logE|BH(X+ΔX〕-BH(X)|2相对log|ΔX|的曲线,由图中可见在一段范围内保持为直线,该范围的上下限,即为|ΔX|min,|ΔX|max。计算参数H和地震数据正态分布的标准差σ,根据公式(4),可以得到关系:
logE|BH(X+ΔX)-BH(X)|2-2Hlog|ΔX|=logб2 (5)
式中:б2=E|BH(X+1)-BH(X)|2,采用最小二乘法解上述方程,即可计算出参数H和б。这种类似二维分形内插,实质上是一种递归中点位移的过程,其递推公式按如下方式进行,对点(i,j),假定当i,j均为奇数时,它的灰度值BH已经确定,则对i,j均为偶数时:
而当i,j中有且仅有一个偶数时,
其中G是Gauss随机变量,服从N(0,1)分布,‖ΔX‖是样本的间距。由此可见,插值点的亮度完全由描述原始数据的分数布郎函数的H和б决定。上述步骤可以重复迭代直到所需的空间分辨率为止。
本发明利用地下唯一的地质体,把地震与测井有机的结合起来,从而使测井曲线延拓到二维,并且采用“局部最优化算法”和WFBR结合的方法,将地震信息和测井信息进行分解、提取、合成、重建等手段,处理出多种质量高,实用性强的储层信息,为油气预测找到了一条新的途径。而且分辨率高,计算精度好,能实现多井闭合,可实现三维处理,更重要的是建立各种测井曲线剖面,用于划分地层岩性、油水界面等。使测井的一孔之见扩展到一个面上,从而为下一步的工作,如提取砂泥岩百分含量、孔隙度、渗透率、饱和度等空间分布提供了依据,也为进一步综合地质解释和储层横向预测工作创造了良好的条件。使解释专家从繁重的计算、推导、绘图中解脱出来,更好的将精力投入到综合分析中去。
附图给出了本发明实施例的有关图示。
图1是地震剖面;
图2是测线联井剖面;
图3是SEIMPAR深度域波阻抗剖面;
图4是SEIMPAR深度域自然电位剖面;
图5是SEIMPAR深度域电阻率剖面;
图6是SEIMPAR自然伽码剖面;
图7是SEIMPAR补偿中子剖面;
图8是A测线地震测面;
图9是A测线SEIMPAR处理波阻抗剖面;
图10是B测线地震剖面;
图11是B测线SEIMPAR处理波阻抗剖面;
图12是本发明处理流程框图。
下面结合实施例与附图,对本发明SEISMPAR作进一步说明。
图12描述了本发明在计算机中实现的流程:
(1)将多个传感器和弹药埋于地层内并引爆弹药或用重锤冲击地层,由埋入地层内的传感器采集到相关数据,(2)将由传感器采集到的数据的地震、测井、地质资料数据输入,(3)计算地震剖面上空间距离为ΔX的数值差的期望值E|BH(X+ΔX)-BH(X)|2;(4)确定尺度极限参数,计算|ΔX|min、|ΔX|max;画出分维图,即logE|BH(X+ΔX)-BH(X)|2相对log|ΔX|的曲线;(5)计算参数H和地震数据正态分布的标准差δ;(6)得到各种测井曲线剖面。
实施例1:该方法首次在新疆吐鲁番-哈密地区某油田实际资料上应用,取得了可喜的成果。图1是该油田一条地震测线,由于受野外采集条件的限制,该资料在信噪比和分辨率方面均不理想,若采用传统处理方法进行储层预测,其前提条件都不能满足。而本技术是以“信息优化论”为基础,能够从复杂的地震资料中提取其共性参数。图2是利用该测线上的九口实际测井资料所得的连井剖面,它反映了地下地层中的砂体分布情况,它是检验我们该技术处理结果正确性的标准,其砂层数目和砂层厚度将做为主要的考查对象。图3是采用该技术处理的波阻抗剖面,其分辨率高,可用于储层评价阶段是储层岩性参数处理的基础,仅靠波阻抗处理还不能准确地预测储层,理由就是波阻抗不能唯一地确定岩性。图4是自然电位剖面,该剖面是在波阻抗基础上,采用本技术处理所得,它非常清楚地反映了砂体的空间分布情况,是地质家解释的可靠性资料,经与图2对比,其符合率达91%。见附表1,2。
表1 两口井反演对比结果
砂层数目(厚度大于3m) | 砂层总厚度(m) | |
实际测井 | 12 | 130 |
反演结果 | 10 | 102 |
符合率 | 83.3% | 76.9% |
表2 三口井反演对比结果
图5是在波阻抗基础上采用该技术处理所得到的电阻率剖面,该剖面主要反映油气水情况,与自然电位剖面联合可用于解释油水界面。在砂层中,低值区说明可能为含水区,砂质较纯,高值区说明可能为含油气区。如图中2420米深度处,正好与图2中的油水界面相符。图6是该技术所得到的自然伽玛剖面,反映泥质含量情况,负异常为砂岩。图7是该技术所得到的补偿中子剖面,该信息与孔隙度、密度有关,高值可能对应孔隙度大。总之,该技术可反演得到各种测井曲线参数剖面,从而将一孔之见扩展到一个面上,为解释专家提供了一种有效的油藏描述手段。
砂层数目(厚度大于3m) | 砂层总厚度(m) | |
实际测井 | 12 | 130 |
反演结果 | 11 | 115 |
符合率 | 91.6% | 88.5% |
实施例2,是甘肃酒东地区某油田的两条十字相交测线的处理结果。图8是测线A的地震剖面。图9是该技术所得到的测线A的波阻抗剖面,从图中可以看出,砂体横向连续,预测可知为沿河道走向。图10是测线B的地震剖面,剖面信噪比较低,连续性差。图11是该技术所得到的测线B的波阻抗剖面。该剖面砂体横向不连续,且表现为透镜状,因此可预测为垂直河道走向,正好与测线A相对应,从而证明了该预测结果的正确性。
Claims (1)
1、一种地震-测井多信息储层数据采集处理方法,其特征在于:(1)将多个传感器和弹药埋于地层内并引爆弹药或用重锤冲击地层,由埋入地层内的传感器采集到相关数据;(2)将由传感器采集到的数据的地震、测井、地质资料数据输入计算机;(3)由计算机计算地震剖面上空间距离为ΔX的数值差的期望值E|BH(X+ΔX)-BH(X)|2;分数Brown运动BH(t)是一非平稳的具有均值为零的Gauss随机函数,其定义如下:
式中BH(t)是通常的Brown运动,0<H<1,当
时,FBM就是Brown运动,Pentland把FBM扩展到高维情况,给出了分数布郎随机场定义:设:X,ΔX∈R2,0≤H<1,F(y)是均值为零的高斯随机函数,Pr(·)表示概率密度,‖·‖表示范数,若随机场BH(X)满足:
则,BH(X)为分数布郎随机场(FBR),H参数反映了地震剖面的粗糙度,由H参数值可得地震剖面的分形维数D为
D=DT+1-H (3)
式中,DT为地震剖面的拓扑维数;BH(X)具有如下性质:E|BH(X+ΔX)-BH(X)|2=E|BH(X+1)-BH(X)|2·‖ΔX‖2H (4)(4)确定尺度极限参数|ΔX|min,|ΔX|max;确定方法为:画出分维图,即:
logE|BH(X+ΔX)-BH(X)|2相对log|ΔX|的曲线,由图中可见在一段范围内保持为直线,该范围的上下限即为|ΔX|min,|ΔX|max;(5)计算参数H和地震数据正态分布的标准差σ,根据公式(4),可以得到关系:
logE|BH(X+ΔX)-BH(X)|2-2Hlog|ΔX|=logб2 (5)
式中:б2=E|BH(X+1)-BH(X)|2,计算出参数H和б;
这种类似二维分形内插是一种递归中点位移的过程,其递推公式按如下方式进行,对点(i,j),假定当i,j均为奇数时,它的灰度值BH已经确定,则对i,j均为偶数时:
而当i,j中有且仅有一个偶数时,
其中G是Gauss随机变量,服从N(0,1)分布,‖ΔX‖是样本的间距;上述步骤可以重复迭代直到所需的空间分辨率为止;
(6)得到各种测井曲线剖面。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 96118790 CN1073705C (zh) | 1996-10-29 | 1996-10-29 | 地震-测井多信息储层数据采集处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 96118790 CN1073705C (zh) | 1996-10-29 | 1996-10-29 | 地震-测井多信息储层数据采集处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1164651A CN1164651A (zh) | 1997-11-12 |
CN1073705C true CN1073705C (zh) | 2001-10-24 |
Family
ID=5125352
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 96118790 Expired - Fee Related CN1073705C (zh) | 1996-10-29 | 1996-10-29 | 地震-测井多信息储层数据采集处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN1073705C (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100429527C (zh) * | 2005-03-25 | 2008-10-29 | 中国石油天然气集团公司 | 三维地震资料处理质量监控技术 |
CN101329405B (zh) * | 2007-06-20 | 2011-02-09 | 中国石油天然气集团公司 | 一种简单的多参数地震反演方法 |
CN101853522B (zh) * | 2010-04-30 | 2012-03-21 | 中国矿业大学(北京) | 三维地质建模系统及其建模方法 |
CN101915088B (zh) * | 2010-06-25 | 2012-10-03 | 中国石油天然气股份有限公司 | 一种油气运移路径生成方法与装置 |
CN102174888B (zh) * | 2011-03-03 | 2012-05-23 | 康志勇 | 一种基于储集岩岩性参数的地层数据处理方法 |
CN103033846B (zh) * | 2011-10-10 | 2016-06-08 | 中国石油化工股份有限公司 | 地质相控制的地震反演系统和地震反演方法 |
CN103375166A (zh) * | 2012-04-26 | 2013-10-30 | 中国石油集团长城钻探工程有限公司 | 一种确定地层中多种矿物组分含量的方法 |
CN102914797B (zh) * | 2012-10-16 | 2015-04-08 | 中国石油天然气股份有限公司 | 一种获得地层各向异性系数的方法及装置 |
CN103605160A (zh) * | 2013-12-09 | 2014-02-26 | 中国科学院武汉岩土力学研究所 | 深埋隧道开挖过程微震事件空间分形维数计算方法 |
CN104933282B (zh) * | 2014-03-21 | 2017-12-12 | 中国石油化工股份有限公司 | 多层圈闭含油气概率的获取方法和系统 |
CN104345345B (zh) * | 2014-10-30 | 2017-09-01 | 中国石油天然气集团公司 | 一种页岩储层总有机碳toc含量预测方法 |
CN107329171A (zh) * | 2017-06-07 | 2017-11-07 | 中国石油天然气股份有限公司 | 深度域储层地震反演方法及装置 |
-
1996
- 1996-10-29 CN CN 96118790 patent/CN1073705C/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN1164651A (zh) | 1997-11-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101158724B (zh) | 基于偶极小波的储层厚度预测方法 | |
Johnson et al. | Hydrostratigraphic interpretation using indicator geostatistics | |
US6498989B1 (en) | Method for predicting dynamic parameters of fluids in a subterranean reservoir | |
CN1040364C (zh) | 采用统计校准技术导出地质特性的地震记录曲线分析方法 | |
Solano et al. | Cumulative-gas-production distribution on the nikanassin tight gas formation, alberta and British Columbia, Canada | |
CN106368691A (zh) | 基于岩石物理地震信息三维异常孔隙压力预测方法 | |
CN1073705C (zh) | 地震-测井多信息储层数据采集处理方法 | |
Harris | The role of geology in reservoir simulation studies | |
CN113779811B (zh) | 一种构造变形对现今地应力方位扰动的定量分析方法 | |
Hesthammer et al. | The use of dipmeter data to constrain the structural geology of the Gullfaks Field, northern North Sea | |
Lorenz et al. | Subsurface fracture spacing: Comparison of inferences from slant/horizontal and vertical cores | |
CN104698494A (zh) | 一种计算异常地层压力的方法 | |
Taheri et al. | Modeling of the shale volume in the hendijan oil field using seismic attributes and artificial neural networks | |
Prankada et al. | Analysis of wellbore stability by pore pressure prediction using seismic velocity | |
Ozkaya | Geometry of layer-bound fractures based on fracture density and aperture-depth plots from resistivity image logs of deviated wells | |
Winkler et al. | Scale dependent hydraulic investigations of faulted crystalline rocks—Examples from the Eastern Alps, Austria | |
CN103278852B (zh) | 利用地震数据体波形结构特征模型预测油气的方法 | |
Rotimi et al. | Reservoir characterization and modeling of lateral heterogeneity using multivariate analysis | |
CN113376692B (zh) | 致密砂岩气水平井压裂改造方案优化方法及装置 | |
Peijs-van Hilten et al. | Heterogeneity modeling and geopseudo upscaling applied to waterflood performance prediction of an incised valley reservoir: Countess YY Pool, southern Alberta, Canada | |
Worthington et al. | Optimizing the value of reservoir simulation through quality-assured initialization | |
CN110795513B (zh) | 河流相源储异位型致密油气甜点区分布的预测方法 | |
Bruno et al. | Geomechanical analysis of pressure limits for gas storage reservoirs | |
Zene et al. | Volumetric estimation and OOIP calculation of the Ronier4 block of Ronier oilfield in the Bongor basin, Chad | |
Hashemi et al. | Enhancing Reservoir Zonation through Triple Porosity System: A Case Study |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C06 | Publication | ||
PB01 | Publication | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C19 | Lapse of patent right due to non-payment of the annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |