[发明专利]地震－测井多信息储层数据采集处理方法

摘要

本发明涉及一种预测地层结构中地下地质体的储存信息的处理方法。其特征在于（１）将多个传感器和弹药埋于地层内并引爆弹药或用重锤冲击地层，由埋入地层内的传感器采集到相关数据，（２）将由传感器采集到的数据的地震、测井、地质资料数据输入，（３）计算地震剖面上空间距离为△Ｘ的数值差的期望值Ｅ丨Ｂ＃-[Ｈ]（Ｘ＋△Ｘ）－Ｂ＃-[Ｈ]（Ｘ）丨＃+[２]；（４）确定尺度极限参数，计算丨△Ｘ丨＃-[min]、丨△Ｘ丨＃-[max]；画出分维图，即logＥ丨Ｂ＃-[Ｈ]（Ｘ＋△Ｘ）－Ｂ＃-[Ｈ]（Ｘ）丨＃+[２]相对log丨△Ｘ丨的曲线；（５）计算参数Ｈ和地震数据正态分布的标准差δ；（６）得到各种测井曲线剖面。$＃!

主权项

1、一种地震-测井多信息储层数据采集处理方法，其特征在于：(1)将多个传感器和弹药埋于地层内并引爆弹药或用重锤冲击地层，由埋入地层内的传感器采集到相关数据；(2)将由传感器采集到的数据的地震、测井、地质资料数据输入计算机；(3)由计算机计算地震剖面上空间距离为ΔX的数值差的期望值E｜BH(X+ΔX)-BH(X)｜2；分数Brown运动BH(t)是一非平稳的具有均值为零的Gauss随机函数，其定义如下：Φ(t,s)={∫-∞0[(t-s)H-12-(-s)H-12]dB(s)+∫0t(t-s)H-12dB(s)}---(1)]]>式中BH(t)是通常的Brown运动，0＜H＜1，当H=12]]>时，FBM就是Brown运动，Pentland把FBM扩展到高维情况，给出了分数布郎随机场定义：设：X,ΔX∈R2,0≤H＜1,F(y)是均值为零的高斯随机函数，Pr(·)表示概率密度，‖·‖表示范数，若随机场BH(X)满足：Pr[BH(X+ΔX)-BH(X)||Δ&KHgr;||H<y]=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)确定尺度极限参数｜Δ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均为偶数时：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·||Δ&KHgr;||·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·||Δ&KHgr;||·H·σ·G---(7)]]>其中G是Gauss随机变量，服从N(0,1)分布，‖ΔX‖是样本的间距；上述步骤可以重复迭代直到所需的空间分辨率为止；(6)得到各种测井曲线剖面。