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

权利要求书

1、一种地震-测井多信息储层数据采集处理方法，其特征在于：(1)将多个传感器和弹药埋于地层内并引爆弹药或用重锤冲击地层，由埋入地层内的传感器采集到相关数据；(2)将由传感器采集到的数据的地震、测井、地质资料数据输入计算机；(3)由计算机计算地震剖面上空间距离为ΔX的数值差的期望值E｜B_H(X+ΔX)-B_H(X)｜²；分数Brown运动B_H(t)是一非平稳的具有均值为零的Gauss随机函数，其定义如下： Φ ( t , s ) = { ∫ - ∞ 0 [ ( t - s ) H - 1 2 - ( - s ) H - 1 2 ] dB ( s ) + ∫ 0 t ( t - s ) H - 1 2 dB ( s ) } - - - ( 1 ) ]]>

式中B_H(t)是通常的Brown运动，0＜H＜1，当 H = 1 2 ]]>时，FBM就是Brown运动，Pentland把FBM扩展到高维情况，给出了分数布郎随机场定义：设：X,ΔX∈R²,0≤H＜1,F(y)是均值为零的高斯随机函数，P_r(·)表示概率密度，‖·‖表示范数，若随机场B_H(X)满足： P r [ B H ( X + ΔX ) - B H ( X ) | | Δ&KHgr; | | H < y ] = F ( y ) - - - ( 2 ) ]]>

则，B_H(X)为分数布郎随机场(FBR)，H参数反映了地震剖面的粗糙度，由H参数值可得地震剖面的分形维数D为

D=D_T+1-H    (3)

式中，D_T为地震剖面的拓扑维数；B_H(X)具有如下性质：E｜B_H(X+ΔX)-B_H(X)｜²=E｜B_H(X+1)-B_H(X)｜²·‖ΔX‖^2H    (4)(4)确定尺度极限参数｜ΔX｜min,｜ΔX｜max；确定方法为：画出分维图，即：

logE｜B_H(X+ΔX)-B_H(X)｜²相对log｜ΔX｜的曲线，由图中可见在一段范围内保持为直线，该范围的上下限即为｜ΔX｜min,｜ΔX｜max；(5)计算参数H和地震数据正态分布的标准差σ，根据公式(4)，可以得到关系：

logE｜B_H(X+ΔX)-B_H(X)｜²-2Hlog｜ΔX｜=logб²    (5)

式中：б²=E｜B_H(X+1)-B_H(X)｜²，计算出参数H和б；

这种类似二维分形内插是一种递归中点位移的过程，其递推公式按如下方式进行，对点(i,j)，假定当i,j均为奇数时，它的灰度值B_H已经确定，则对i,j均为偶数时： B H ( i , j ) = 1 4 { B H ( i - 1 , j - 1 ) + B H ( i + 1 , j - 1 ) + B H ( i + 1 , j + 1 ) ]]> + B H ( i - 1 , j + 1 ) } + 1 - 2 2 H - 2 · | | Δ&KHgr; | | · H · σ · G - - - ( 6 ) ]]>

而当i,j中有且仅有一个偶数时， B H = ( i , j ) = 1 4 { B H ( i , j - 1 ) + B H ( i - 1 , j ) + B H ( i + 1 , j ) + B H ( i , j + 1 ) } ]]> + 2 - H 2 1 - 2 2 H - 2 · | | Δ&KHgr; | | · H · σ · G - - - ( 7 ) ]]>

其中G是Gauss随机变量，服从N(0,1)分布，‖ΔX‖是样本的间距；上述步骤可以重复迭代直到所需的空间分辨率为止；

(6)得到各种测井曲线剖面。