地震-测井多信息储层数据采集处理方法技术

本发明专利技术属于预测地层结构中地下地质体的储层信息的方法，是最优化算法和ＷＦＢＲ方法的结合。将地震信息和测井信息进行分解、提取、合成、重建等手段，反演出多种质量高，实用性强的储层信息，为油气预测找到了一条新的途径。而且分辨率高，计算精度好，能实现多井闭合，可实现三维反演，现重要的是建立各种测井曲线剖面，用于划分地层岩性、油水界面等。使测井的一孔之见扩展到一个面上，从而为下一步的工作，如提取砂泥岩百分含量、孔隙度、渗透率、饱和度等空间分布提供了依据。（*该技术在2016年保护过期，可自由使用*）

【技术实现步骤摘要】
【国外来华专利技术】
本专利技术涉及一种预测地层结构中地下地质体的储层信息的方法，特别是一种地震—测井多信息储层参数反演方法(SEIMPAR)。目前叠后地震反演是一种最常见的反演方法，我国进行了大量的研究工作，并在储层预测等工作上见到了一些成效。早期主要是普通的递归法，要受频带的限制，因此分辨率较低；80年代中、末期出现了稀疏脉冲法，适合于反射系数变化大，有一些突出反射面的地区；模型法则需要严格的约束，以克服多解性问题。在国内，较早引进的方法有Velog、Seislog，近年来引进和研制的slim、Delog，这些方法在反演过程中没有充分利用测井资料的高频信息，因此分辨率较低。所以叠后地震道反演技术存在的问题有纵、横向辨率低；反演结果的多解性严重；横向外推精度不高；多井处理难于闭合；三维反演实现困难；多参数综合反演方法尚未成熟。本专利技术的目的是提供一种专门针对目前叠后地震道反演技术存在的问题，采用“局部最优化算法”和WFBR方法，充分利用地震、测井和地质概念进行综合反演的地震—测井多信息储层参数反演方法。本专利技术是最优化算法和WFBR方法的结合，首次提出地震波场是地下地质体微观波动的结果总和。地震波场是对地下地质体特征的总体反映，测井是对地下地质体特征的局部反映，地震和测井都包含地下地质体的储层信息。其特点是利用地震、测井、地质资料，采用信息优化和WFBR方法，反演出各种测井曲线剖面，然后将多种信息综合，得到储层参数剖面。有别于传统反演方法。传统方法是利用地震、测井、地质资料以褶积模型为基础，采用优化相关外推方法反演出波阻抗剖面，然后利用经验公式或量板，计算出储层参数剖面。本专利技术的实现过程是采用滑动小时窗，首先计算地震剖面上空间距离为△X的数值差的期望值E|BH(X+△X)-BH(X)|2；然后确定尺度极限参数|△X|min，|△X|max；最后计算参数H和地震数据正态分布的标准差σ。分数Brown运动BH(t)是一非平稳的具有均值为零的Gauss随机函数，其定义如下{BH(t)=1Γ(H+12)·Φ(t,s)BH(0)=0]]>Φ(t,s)={∫-∞0dB(s)+∫0t(t-s)H-12dB(s)}---(1)]]>式中BH(t)是通常的Brown运动，0＜H＜1，当 时，FBM就是Brown运动，Pentland把FBM扩展到高维情况，给出了分数布郎随机场定义，设X，△X∈R2，0≤H＜1，F(y)是均值为零的高斯随机函数，Pr(·)表示概率密度，||·||表示范数，若随机场BH(X)满足Pr=F(y)---(2)]]>则，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均为偶数时BH(i,j)=14{BH(i-1,j-1)+BH(i+1,j-1)+BH(i+1,j+1)]]>+BH(i-1,j+1)}+1-22H-2·||ΔX||·H·σ·G---(6)]]>而当i，j中有且仅有一个偶数时，BH(i,j)=14{BH(i,j-1)+BH(i-1,j)+BH(i+1,j)+BH(i,j+1)}]]>+2-H21-22H-2·||ΔX||·H·σ·G---(7)]]>其中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.计算地震剖面上空间距离为△X的数值差的期望值E|BH(X+△X)-BH(X)|2；3.确定尺度极限参数，计算|△X|min、|△X|max；画出分维图，即logE|BH(X+△X)-BH(X)|2相对log|△X|的曲线；4.计算参数H和地震数据正态分布的标准差σ；5.得到各种测井曲线剖面有波阻抗剖面，自然电位剖面，电阻率剖面，自然伽玛剖面，补偿中子剖面。实施例1该方法首次在新疆吐鲁番—哈密地区某油田实际资料上应用，取得了可喜的成果。图1是该油田一条地震测线，由于受野外采集条件的限制，该资料在信噪比和本文档来自技高网...

【技术保护点】
一种地震－测井多信息储层参数反演方法包括局部最优化算法和ＷＦ－ＢＲ方法，用于从地震、测井资料中提取储层多种信息，其特征在于采用滑动小时窗，按下列步骤进行：ａ．计算地震剖面上空间距离为ΔＸ的数值差的期望值Ｅ｜Ｂ↓［Ｈ］（Ｘ＋ΔＸ）－Ｂ↓［ Ｈ］（Ｘ）｜↑［２］；ｂ．确定尺度极限参数｜ΔＸ｜ｍｉｎ，｜ΔＸ｜ｍａｘ；ｃ．计算参数Ｈ和地震数据正态分布的标准差σ。

【技术特征摘要】
【国外来华专利技术】

【专利技术属性】
技术研发人员：梁秀文，撒利明，晁吉俊，张志让，高建虎，张金玲，雍学善，石兰亭，许云泽，
申请(专利权)人：中国石油天然气总公司西北地质研究所，
类型：发明
国别省市：62[中国|甘肃]